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specification  of  processing  functions  for  nodel  representation. 
Example  fanctions  used  in  this  worX  are  computation,  data  base 
management,  communications  and  terminal/access  interface  as 
would  be  applicable  to  a military  command  and  control  computer 
system.  The  models  may  be  used  for  trade-off  analysis  of 
critical  performance  criteria.  Additionally,  a basic  method  for 
the  synthesis/determination  (design)  of  processor  "computer 
power"  to  ensure  a load  balanced  network  is  ascertained.  The 
three  models  and  their  respective  computer  network  evaluations 
for  trade-off  analysis  are  pursued  with  their  inherent 
advantages  and  disadvantages.  Additionally,  these  models  are 
compared  with  a classical  event  driven  simulation  in  order  to 
demonstrate  their  internal  consistency  and  to  provide  a basis 
for  verification.  All  three  models  accommodate  job  routes 
(types)  with  one  model  providing  exact  solutions  for  each  job 
type.  These  models  accommodate  large  numbers  of  jobs  so  that  an 
effective  total  computer  network  may  be  modelled.  The  sole 
restriction  is  the  computatiocal  complexity.  In  order  to 
implement  the  models,  several  computer  program  modules  used  for 
various  solution  stages  of  the  models  are  provided.  These 
nodules  include  job  typing  or  classes,  solution  of  normalizing 
constants,  and  reducticn  to  an  equivalent  network  of  two  nodes. 
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NCTATICNAL  CCNVENTICNS 

In  order  to  provide  clarity  and  correspondence  with  the 
coaputer  program  modules  developed  in  support  of  this  work  and 
to  facilitate  the  preparation  of  this  dissertation,  a particular 
set  of  recognized  nctational  conventions  has  been  adopted. 
Since  the  chosen  conventions  are  those  normally  employed  in 
scientific  computer  programming  languages,  such  notation  should 
be  no  impediment  to  the  experienced,  scientific  reader  and  will 
be  an  asset  to  the  experienced  computer  user  and  professional. 

Of  particular  importance  is  the  method  for  denoting 
subscripts,  superscripts,  and  functions.  Subscripts  are  given 
within  parentheses  such  as  F(i,j)  fcr  p , and  functions 
similarly  denoted  so  that  Tr(N)  indicates  that  Tr  is  a function 
of  N.  Efforts  have  been  made  to  avoid  confusion  between 
parenthetical  subscripts  and  functions  by  explicitly  pointing 
out  the  function  names  as  they  are  developed.  Where 

superscripts  might  normally  be  used  as  in  p , the  superscript 
indicator  k has  been  appended  to  the  function  letter  p to  yield 
pk.  As  stated  above,  such  a nctational  anomaly  is  one  familiar 
to  the  computer  user.  Furthermore,  where  the  order  of 
operations  in  an  eguation  is  not  explicit,  the  and  the 
have  been  used  as  the  symbols  for  multiplication  and 
exponentiation,  respectively.  The  standard  operator  symbols  for 
addition  (•♦■),  subtraction  (-)  , and  division  (/)  ate  also  used. 


L 


X 


Precedence  of  operators  is  as  in  prograaaing  languages 


Onotdered  sets  are  denoted  fcy 


braces  "{»•  and  " and 


ordered  sets  by  the  deliaitecs  "{!"  and  Vector  and  natrix 


names  are  underlined  as  in  P.  Elements  c£  vectors  and  matrices 
are  enclosed  within  square  brackets  as  in  [23  45  43]. 
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INTPCEOCTIOM 

This  research  deals  with  the  establ ishment  of  an  analytic 
model  of  a projected  computer  network  employing  the  techniques 
of  Network  Queuing  Theory.  Such  a model  is  referred  to  as  an 
Exponential  Queuing  Network  (EQU)  Model.  The  projected  computer 
network  is  one  based  on  identified  types  of  processing 
functions.  The  functions  identified  are  at  the  discretion  of 
the  analyst  or  designer.  However,  in  this  work  the  functions  of 
computation,  data  base  management,  communications,  and 
terminal/access  interface  have  been  selected.  The  analytic 
models  which  have  been  established  can  assist  in  the  formulation 
of  various  network  structures  employing  defined  processing 
functions  through  subsequent  trade-off  analysis  of  critical 
performance  criteria.  Three  methods  of  evaluating  the  model  are 
presented  along  with  their  advantages  and  disadvantages. 
Additionally,  a basic  method  for  the  synthesis/determination 
(design)  of  processor  "computing  power"  to  ensure  a load- 
balanced  computer  network  is  ascertained  through  the  use  of  the 
model.  Consequently,  the  models  provide  a quantitative  means  or 
tool  for  both  a basic  design  and  a thorough  evaluation,  or 
trade-off  analysis,  of  various  computer  network  organizations. 

The  four  types  of  processing  functions  selected  have  been 
based  on  a projected  concept  of  processing  functions  that  can  he 
identified  within  a geographical  command  and  control  system 
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suited  to  the  ailitary  environment.  These  functions  are  seen 
from  the  user  viewpoint  and  consist  of; 

Front-end  Interface  Processing  - That  processing 
associated  with  driving  and  interfacing  terminals 
similar  to  the  Terminal  Interface  Processors  (TIP)  of 
the  AHPA  Network  [Roberts  73]  and  access  control  for 
authorized  users  and  levels  of  usage. 

2.  Ba§e  Hanagement  Processing  - That  processing 
concerned  with  the  update,  retrieval,  and 
creation/maintenance  of  files  comprising  the  data  base 
as  would  be  associated  with  the  data  base  machine 
concept  {[ACB  77],  [Canaday  7«]). 

3.  Communications  Processing  - That  processing  required  to 
act  as  a node  within  a network  or  the  interconnection 
point  for  the  passing  of  data  traffic  to  other 
functional  processors  similar  to  the  concept  of  the 
IMPS  (Interfact  nessage  Processors)  of  the  ARPA  Network 
[Roberts  73]. 

PliPCAPal  Sl  Processing  - That  processing 

normally  associated  with  the  ability  to  do  high  speed 
arithmetic,  "number  crunching"  operations. 

Geographical  separation  or  co-location  and  combinations  thereof 
for  the  processing  functions  could  be  determined  by  operational 
constraints  such  that  any  of  the  following  may  arise: 

1.  All  functions  nay  be  in  separate  locations. 

2.  Where  functions  are  co-located,  the  transmission  medium 

\ 

\ 

i 

I 

i 


Introduction 


1 


I 

I 


I 

i 

1 

[ 

I 


1 


may  be  high  speed  parallel  transfers  such  as  those 
within  a classical  mul ti/parallel  processing 
envircnnent. 

3.  When  all  four  functions  are  physically  co-located,  this 
organization  may  be  viewed  as  a user  f unction- oriented 
multiprocessing  system.  Cne  might  view  this  in  the 
sense  of  a loosely  coupled  system  |LCS)  or  a local 
network. 

4.  Numbers  and  capabilities  cf  processing  functions  would 
vary.  Selection  of  numbers  and  capabilities  is 
dependent  to  a large  degree  on  operational  requirements 
and  other  factors  such  as: 

a.  Fault  tolerance. 

b.  Pile  allocation. 

c.  Existing  communication  facilities. 

d.  Existing  organizational  units,  e.g.  , headquarters, 
etc. 

e.  Hierarchical  structure  of  the  supported 
organization,  e.g,,  chain  of  command  or  updating 
requirements. 

5.  Evaluation  is  based  on  the  total  network  organization 
to  include  the  relative  performance  of  each  functional 
processor  and  not  just  the  communication  subnetwork  in 
support  of  the  computer  network. 

The  applications  for  the  model  as  developed  fall  into  the 
categories  of  analysis  and  of  design.  Three  different  methods 
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of  aodel  evaluation  for  trade-off  analysis  are  pursued  with 
their  inherent  advantages  and  disadvantages.  The  design 
techu.-'ue  is  applied  with  the  nost  basic  of  the  three  nodel 
evale<<r  rs.  Additionally,  these  nodels  are  coapared  with  a | 

classical  event  driven  siaulation  in  order  to  deaonstrate  their 
internal  consistency  and  to  provide  a basis  for  verification. 
Analysis  and  design  are  demonstrated  such  that  an  existing 
system  nay  be  parameterized  in  terns  of  functions  to  natch  the 
model  and  subsequently  evaluated  in  terns  of  the  najot 
performance  criterion  of  response  time,  \iith  certain  guidelines 
in  nind,  the  architect  nay  also  fornulate  a topology  and 
ascertain  the  necessary  processing  power  (design  process) 
required  to  achieve  a specified  response  tine  with  a specified 
number  of  users. 

An  additional  result  of  this  wort  is  the  provision  of 
computer  program  nodules  which  are  used  for  various  solution 
stages  in  the  inplementation  of  the  models.  The  program  nodules 
include  the  provision  for  job  typing  or  classes,  solution  of 
normalizing  constants,  and  reduction  to  an  equivalent  network  of 
two  nodes.  They  may  also  be  implemented  on  an  interactive  basis 
using  the  IBB  TSO. 
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Chapter  I 5 

Chapter  I:  EACKGPCOND  CEVELCPHENT 

Oevelcpaent  of  Parallel  Processing  Computer  Systems 

We  may  commence  the  developmental  review  by  ascertaining 
the  characteristics  of  a multiprocessing  or  parallel  processing 
system.  The  characteristics  may  vary  frcm  designer  to  designer 
but  in  general  are  as  enumerated  in  Snslcw  [197tt]; 

1.  Two  or  more  processors. 

2.  Operation  of  the  processors  in  a coordinated  manner  at 
the  hardware/sof twar e level. 

3.  Sharing  cf  common  main  memory. 

4.  Operation  under  the  control  cf  a single  integrated 
operating  system  which  might  even  be  viewed  as 
operation  under  a single  control  structure  such  as 
would  be  the  case  for  a computer  network  system 
organ ization. 

With  these  four  characteristics  ia  mind  it  is  natural  to 
examine  the  general  rationale  for  multiprocessing  or  parallel 
processing.  There  are  several  motivating  factors  which  have 
brought  about  this  type  of  processing.  Seme  of  these  factors 
are; 

1.  Sharing  of  resources. 

2.  Increasing  the  computaticnal  oi  processing  capability. 

3.  Providing  concurrence  cf  operations  to  achieve 
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increased  speed. 

4.  Recognizing  the  speed  limitations  of  the  physical 
hardware  prompting  the  need  for  other  computational 
means. 

5.  Providing  a fault  tolerant  capability. 

Looking  at  each  of  the  points  enumerated,  we  find  that  the 
motivating  factors  have  indeed  teen  significant.  The  sharing  of 
resources  has  been  important  due  to  their  expense,  and  the 
consequent  need  to  minimize  these’  costs.  It  has,  therefore, 
been  impractical  to  allcw  any  one  resource  to  remain  idle  when 
another  process  can  utilize  it.  The  resources  that  are  referred 
to  include  the  entire  spectrum  of  "classical"  computer  resources 
such  as  processors,  memory,  devices,  aid  information  [Hadnick  £ 
Donovan  74]. 

The  gain  in  computation  and  processing  capability  by 
permitting  an  overlap  of  operations  and  of  processes,  and  in 
particular  of  servicing  many  users,  is  readily  exemplified  by 
the  Hultics  system  ([Crganick  72],  [Sekino  72]).  This  is 
certainly  not  the  only  example  of  such  a gain.  Nearly  every 
major  computer  manufacturer  has  produced  either  or  both  general 
or  special  purpose  systems  in  such  a category  [ Enslow  74]. 

Not  only  has  the  sharing  of  resources  had  a significant 
impact  on  the  development  of  parallel  processing  systems,  but 
the  recognition  of  the  inherent  speed  limitations  of  the 
physical  hardware  devices  is  significant.  Switching  speeds  of 
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devices  are  now  operating  at  the  level  of  10  **  -10  sec  and  a 
comaonly  accepted  liaitation  due  to  the  heat  dissipation  is  at 
the  level  of  10  **  -11  sec.  Furthermore,  even  if  switching 
speeds  were  reduced  to  zero,  propagation  delays  would  become 
dominant.  Yet,  propagation  delays  are  a function  of  device 
circuit  densities  which  are  fast  approaching  the  limits  of 
optical  resolution  [Stone  75],  consequently,  we  cannot  expect 
to  find  much  reduction  nor  will  we  be  able  to  gain  any 
significant  increase  in  speed  from  the  physical  device.  The 
need,  then,  is  to  find  other  imaginative  or  innovative  means. 

Certainly  the  replication  of  similar  units,  or  even  the 
distribution  of  processes  among  dissimilar  units  for  execution, 
leads  to  a certain  amount  of  fault  tolerance.  Typical  straight 
forward  approaches  are  achieved  by  tie  multiple  and  identical 
central  processing  units  of  such  systems  as  the  Honeywell  60CC 
and  the  Univac  1100  systems  [ Enslow  74].  Mote  sophisticated 
systems  for  fault  tolerance  ate  achieved  ty  various  schemes  for 
the  networking  cf  computet  systems  ([Cosell  75],  [Doll  73]). 

The  structure  of  the  operating  system  or  of  a control 
system  to  support  the  different  architectural  approaches  tc 
multi/parallel  processing  systems  has,  quite  naturally,  evolved 
and  adapted  to  meet  the  varied  demands  cf  the  systems.  Perhaps 
this  point  is  best  illustrated  by  the  lEt!  approach  first  through 
the  LCS  (loosely  coupled  system)  and  then  the  DCS  (direct 
coupled  system)  [Freeman  76].  In  the  LCS  it  is  questionable 
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whether  actual  aultiprocessing , in  the  strictest  sense,  is 
occurring,  in  view  of  the  fact  that  the  criterion  of  operation 
under  the  control  of  a single  integrated  operating  systea  are 
not  in  fact  aet.  However,  if  we  note  that  the  two  or  more 
processors  must  operate  in  a controlled  fashion  relative  to  each 
other,  similar  to  the  network  concept,  then  it  is  feasible  to 
consider  such  a systea  in  the  category  of  parallel  processing 
and  hence  in  a broad  category  of  aulti/parallel  processing. 

Furthermore,  the  characteristics  of  the  operating  system 
itself  have  taken  on  various  forms.  The  initial  approaches  of 
tying  together  cooperating  modules  to  accomplish  the  various 
functions  (as  in  OS/360  [nadnick  & Donovan  74])  have  been 
supplanted  by  more  elegant  approaches.  Some  of  these  approaches 
such  as  the  RC-4000  [Hansen  73]  or  Hoare's  monitor  concept 
[ Hoare  74]  have  provided  a more  positive  means  of  control  and 
assurance  of  cooperative  interaction  of  system  functions. 
Perhaps  the  culaination  of  such  complex  operating  systems  is 
represented  by  the  Rultics  system  [Ctganick  72].  This  system 
represents  the  current  epoch  for  a multi-user,  multiprocessor 
operating  systea. 


We  must  note,  however,  that  throughout  the  development  of 
parallel  processing  systems  there  is  always  the  recognition  and 
identif icaticn  of  the  specific  functions  which  must  he 
accomplished.  For  instance,  we  have  already  acknowledged  the 
’’classical"  functions  of  the  operating  system  for  such  systems  - 
processor  management,  memory  management,  device  management,  and 
information  management  [ ffadnick  S Conovan  74],  These  functions 
are  further  subdivided  into  different  mechanisms  to  implement 
various  policies  which  have  teen  valuable  to  a particular 
computer  system  and  set  of  users.  Eut  even  more  than  this, 
identifications  of  other  speciali2ed  functions  associated  with 
the  parallel  environment  have  surfaced.  Two  such  functions  are 
associated  with  the  requirement  for  the  computer  system  to 
communicate  with  users  or  other  computers  and  the  need  to  handle 
vast  amounts  of  information. 

It  has  been  well  established  that  the  concept  of  "front- 
ending" a main  computer  system  by  a computer  to  accomplish  the 
specialized  function  of  communications  to  include  protocols  and 
message  preparation  can  most  effectively  be  performed  by  a 
functionalized  computer  [Wood  73].  in  effect  it  would  appear 
that  this  approach  very  closely  resembles  the  IBfl  concept  of 
LCS.  However,  because  of  the  requirement  to  communicate  and  the 
need  for  the  transmission  of  various  amounts  of  data  to  include 
the  support  of  many  users  sending  small  amounts  of  data,  it  has 
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teen  determined  to  be  more  productive  tc  employ  a functionalized 
computer.  In  fact  this  approach  is  taken  with  many  of  the 
current  systems  to  include  the  Honeywell  Hultics  600C  and  the 
CDC  600C/7000  series  [ Enslow  74], 

The  other  function  that  has  received  recognition  in  the 
recent  past  is  the  concept  of  a "back-end"  computer  or 
essentially  a data-base  machine  ([Canaday  74],  [Baum  6 Nsiac 
76],  [ farill  and  Stern  75]).  This  concept  is  analogous  to  that 
of  the  front-end  computer  development.  In  the  former  it  was 
recognized  that  a large  processing  load  was  being  generated 
solely  to  support  the  communications  traffic  and  it  was  more 
efficient  for  a specialized  processor  to  perform  these  tasks. 
Likewise,  the  increased  traffic  or  processing  loads  encountered 
solely  to  access,  maintain  and  retrieve  information  from  data 
bases  has  been  identified  as  being  a candidate  for  specialized 
processing  ([Aschim  73],  [ Baryanski  75]).  In  essence,  we  are 
beginning  to  see  a move  towards  identification  of  specific 
functions  to  be  accomplished  within  the  parallel  processing 
system  and  the  assignment  of  these  tasks  to  specific  processors 
within  the  system.  In  other  words,  it  is  more  effective  for 
these  specialized  functions  to  be  accomplished  within  a 
dedicated  processor  with  an  architecture  tailored  to  the 
reguirement  rather  than  by  a system  of  multiple,  general  purpose 
processors  which  are  essentially  homogeneous. 
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Networking  of  Computer  Systems 

The  factors  of  functionalized  development  and  the 
recognized  limitation  for  a totally  homogeneous  parallel 
processing  system  to  be  able  to  meet  the  increasing  reguirements 
[Flynn  72]  have  supported  the  development  of  the  networking 
concept.  With  networking  we  are  no  longer  concerned  that  the 
entire  network  operate  under  the  control  of  a single  integrated 
operating  system.  We  are  concerned,  however,  with  the 
cooperative  sharing  of  resources  amcng  individual  computer 
systems.  In  essence,  "a  computer  network  can  be  defined  as  an 
interconnected  group  of  host  computers  that  automatically 
communicate  with  one  another  and  that  can  share  such  resources 
as  programs,  data  bases,  memory  space,  and  the  long-haul  links 
with  each  other.  In  a true  network  situation,  therefore,  each 
constituent  computer  system  can  either  operate  in  a local  mode 
under  its  own  operating  system  or  participate  in  network 
activity  under  the  direction  of  a higher-level  network 
supervisor  program,  or  it  can  do  both"  [Coll  73].  It  can  be 
observed  that  this  does  not  necessarily  imply  that  there  is,  or 
is  not,  a geographical  separation  of  computers.  The  commonly 
accepted  thoughts  of  networking  were  such  that  we  generally  took 
the  systems  to  be  physically  separated.  However,  this  may  not, 
nor  need  not,  be  the  case.  Furthermore,  we  may  note  that  in  the 
network  operaticn,  control  is  now  under  a higher-level  network 
supervisor  program.  The  supervisor  program  may  not  need  tc 
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concern  itself  with  the 
for  specific  purposes 
system. 


detailed  control  of 
as  is  the  case  with 


special  resources 
the  local  operating 


At  this  point  it  is  interesting  to  note  some  of  the 
motivating  factors  for  networking  and  realize  the  similarity  to 
the  motivating  factors  for  multi/paralle 1 processing: 

1.  Resource  sharing  - particularly  for  files  and  data. 

2.  Greater  computational  capability. 

3.  Load  sharing. 

4.  Increased  efficiency  and  use  of  resources. 

5.  Fault  tolerance. 

Some  of  these  factors,  if  not  identical,  are  quite  similar  to 
those  postulated  earlier  for  the  raulti/parallel  processing 
systems.  In  fact,  at  this  point,  we  can  state  that  the  concept 
of  networking  is  merely  another  way  of  looking  at  the 
multi/parallel  processing  concept  hut  perhaps  on  a larger,  macro 
scale.  Being  cognizant  of  this  situation,  we  are  now  ready  to 
examine  those  architectural  changes  that  are  occurring  in  these 
systems. 
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Architectural/Crgani^ational  Ccnsiderations 

■^he  introduction  of  low  cost  ccb futational  lodules  has 
brought  forth  auch  new  or  revised  thinking  in  the  area  of  system 
architecture  and  organization  and  hence  in  its  functioning.  In 
fact,  it  is  difficult  to  separate  what  might  be  classified  as  a 
multi/parallel  processing  system,  or  even  a network  of  several 
such  systems  and  whether  or  not  they  are  locally  interconnected 
or  geographically  separated  throughout  the  world.  This 
situation,  however,  is  no  longer  the  primary  factor  as  it  is  not 
germaine.  Hhat  is  essential  is  the  identification  of  the 
functions,  the  assignment  of  the  fancticns  to  processing  units 
and  the  necessary  control  and  interaction.  Recently  there  has 
been  a plethora  of  concepts  and  proposed  architectures  embracing 
the  notion  of  low  cost  computational  modules  [IEEE  Symp  77]. 

To  note  several  of  the  more  established  examples  we  could 
begin  by  examining  the  structure  of  the  Carnegie-Mellon  C.mmp 
project  ([Cohen  S Jefferson  75],  [levin  et  al  75],  [Mulf  et  al 
75]).  This  architecture  is  essentially  a crossbar  structure 
with  15  processor  units  and  16  memory  units.  Each  of  the 
processors  may  be  assigned  any  of  the  identified  tasks  within 
the  system.  There  is  no  attempt,  other  than  within  the 
protection  features  themselves,  to  assign  any  specialized 
function  to  a specific  processor.  Consequently,  this  approach 
appears  to  be  no  different  than  what  has  been  encountered 
earlier.  Perhaps  the  large  number  of  available  processors  to  do 
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tasks  and  the  lessening  requirewent  for  a processor  to  switch 
from  one  task  to  another  once  initiated  and,  therefore,  a 
resulting  lowering  of  overhead  would  be  an  exception.  However, 
as  the  authors  readily  point  out,  the  operating  system  itself. 
Hydra,  is  as  complex  and  as  comprehensive  as  any  encountered  to 
date.  Another  point  to  note  i.«  that  this  system  is  a local 
system  which  participates  in  the  AFPA  Network. 


Another  approach  is  taken  wit 
Hult  iprocessor  ([Eressler  et  al],  [C 
architectural  concepts,  the  process 
interconnected  via  multiple  tusses 
Although  the  operating  system  that  is  em 
consistent  treatment  of  all  processo 
parts  of  the  software,  which  are  known  t 
resources,  are  provided  as  multiple  co 
associated  with  a particular  processor, 
provide  for  a certain  degree  of  function 
not  identified  as  such.  The  main  advant 
others  seems  to  be  the  fault  tolerance 
centralized  switch  to  fail.  Although  it 
processes  are  functionalized,  there 
equal  treatment  that  implies  a sophistl 
for  the  control  and  management  of  system 


h the  BBN  Pluribus 
osell  75  ])  . In  terms  of 
ors  and  memory  are 
as  the  name  implies, 
ployed  provides  for  a 
rs  as  equal  units,  some 
c make  heavy  demands  on 
pies  in  a private  memory 
In  effect  this  seems  to 
alization  although  it  is 
age  of  this  system  over 
, as  there  is  not  even  a 
would  appear  that  some 
is  still  a generality  of 
rated  operating  system 
functions. 
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The  macro-module  or  multiple  function  unit  architecture  is 
which  makes  a more  refined  step  towards  functionalization. 
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At  least  the  recognition  is  made  of  the  possible  utilization  and 
advantages  of  identifying  and  iapleaenting  particular  functions. 
Such  is  the  case  of  the  multiple  function-unit  processors  as 
described  by  Eovra  [197ft]  and  Crnstein  [19ft7]  and  implemented 
in  such  machines  as  the  CDC  ftftOO  and  the  IBil  360/91  [ Enslow  7u]. 
These  approaches  go  to  the  other  end  of  the  spectrum  since  they 
assign  the  computational  functions  to  specific  units  or  modules 
and  hence  imply  very  little  change  in  the  operating  system  for  a 
large  scale  machine.  Instead  they  are  concerned  primarily  with 
the  hardware  means  cf  implementation.  We  cannot  overlook, 
however,  the  associated  problem  of  identification  of  the 
parallel  computations  and  the  necessary  control  criteria. 

Lest  one  feel  that  these  are  the  only  structural  approaches 
that  may  be  taken  we  cculd  examine  a host  of  others  such  as  ting 
structures  [ Parber  7*5],  hierarchical  structures  [ Aschenhurst 
75],  and  virtual  channel  structures  [Fraser  75].  Even  some  of 
the  architects  of  C.mmp  are  studying  systems  composed  of  large 
numbers  of  microcomputers  and  clusters  of  microcomputers  for 
multiprocessor  system  structures  ([Jones  et  al  75],  [Fuller 
76  ])  . 

The  most  interesting  and  possible  ma^or  differing  structure 
is  that  proposed  by  Arden  [1975].  Employing  the  concept  of 
modularization  and  functionalization,  "it  should  be  possible  to 
construct  a system  that  is  at  least  as  powerful,  is  highly 
modular  in  construction,  yet  permits  the  use  of  much  simpler 
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operating  system  structures"  [Arden  75].  The  limitation  due  to 
cost,  complexity,  and  reliability  of  merely  constructing  bigger 
and  faster  central  processors  has  been  recognized.  It  appears 
that  in  the  conceptual  approach  presented  by  Arden,  the 
operating  system  is  simplified  by  the  transfer  of  many  control 
functions  to  hardware.  That  is  to  say  that  "by  identifying  the 
state  of  a process  with  a processor,  a number  of  process  control 
tables  and  task  dispatcher  routes  are  obliterated"  [Arden  75]. 
It  appears  that  the  control  structure  is  in  fact  simplified  and 
more  easily  associated  with  a processor.  This  is,  in  itself,  an 
important  trend  and  implies  that  the  more  functionalization  that 
can  be  achieved  perhaps  the  simpler  the  control  structure  can 
be.  Perhaps  one  can  even  arrive  at  a pcint  whereby  seme  of  the 
classical  schemata  may  be  apnlied  [Karp  5 Riller  69]  in  order 
that  correct  control  is  assured  and  models  for  deterministic 
evaluation  may*  be  developed. 

It  is  interesting  to  observe  that  the  trend  described  above 
may  result  in  a simplified  means  cf  implementing  a total 
computer  system,  which  includes  the  operating  system,  whether  it 
is  applied  to  a single  local  computer  processing  system  or  to  a 
computer  network  of  processing  systems.  In  effect  it  is  left  to 
the  definition  of  the  processing  function  being  performed.  It 
could  be  on  a very  low  level  such  as  at  the  computation  unit  or 
at  the  higher-level  cf  a front-end  or  back-end  (data-base 
machine)  computer  function,  either  local  or  networked.  As  a 
consequence,  we  have  now  arrived  at  a pcint  whereby  we  have  both 
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a reorientation  of  the  system  architecture  and  organization  and 
the  associated  operating  system  from  that  typified  by  Multics  to 
that  presented  by  Arden  [1975]  which  will  result  in  more 
emphasis  on  the  classical  forms  of  schemata  for  both  a means  of 
controlling  various  functions  and  as  a means  of  evaluation  of 
the  functional  organization  of  such  a system. 
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Otilization  of  Functionalized  Digital  Processors 


Application  of  the  above  trends  to  the  developnent  of 
Coaaand  and  Control  Systeas  used  by  the  National  Coaaand 
Authorities  appears  to  be  a particularly  fertile  area.  Already 
developaent  and  exper iaentation  are  underway  for  a preliainary 
network  of  coaaand  systeas  [CCTC  76].  Cne  of  the  aajor  probleas 
encountered  is  that  of  evaluating  a system  of  network  structures 
to  satisfy  reguireaents  short  of  building  and  evaluating  the 
systea  or  even  a long,  involved  and  detailed  siaulation  of  the 


systea.  If,  however. 

such 

architectural 

approaches 

as 

previously  presented 

are 

examined  hand 

in 

hand 

with 

functionalization,  it  is 

quite 

conceivable 

that  a 

means 

of 

evaluating  the  structure  of  the  proposed  systea  in  terms  of 
quantitative  measures  can  be  ascertained  on  a comparative  basis 
in  respect  to  another  structure.  Through  the  use  of  analytical 
modelling  as  a means  of  evaluation  we  can  gain  the  advantages  of 
flexibility,  relatively  low  cost,  responsiveness  and  significant 
information  capacity  ([Coletta  et  al  74],  [Klein  cock  75]). 

When  we  examine  the  idea  of  functionalization  within  the 
parallel  system  what  we  are  really  locking  for  is  the  trade-off 
of  hardware  costs  versus  complexity  costs.  Even  more  subtle  is 
the  fact  that  complexity  impinges  on  reliability  and  hence  on 
fault  tolerance.  Presumably  the  processing  hardware  costs  are 
constantly  being  driven  downwards  so  that  functionalization  is 
more  and  more  attractive.  As  a result,  we  must  reevaluate  some 
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of  the  factors  that  originally  motivated  the  classical  move 
towards  the  use  of  multi/parallel  processing  under  control  of  a 
complex  operating  system.  A primary  factor  for  examination  is 
that  of  resource  sharing.  In  actuality  the  principal  resource 
that  must  be  shared  frcm  an  operational  point  of  view  in  a 
military  command  and  control  situation  is  that  of  information. 
Next  a high  level  of  fault  tolerance  must  be  achieved, 
yunctionalization  within  a worldwide  network  seems  to  be  the 
most  obvious  answer  whereby  we  can  achieve  the  following  desired 
results : 

1.  Informaticn  resource  sharing 

2.  Fault  tolerance 

3.  Greater  capability 

u.  Load  sharing 

5.  Greater  efficiency 
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Network  Functions  and  Applications 

The  functionalization  within  a network  lends  itself  to  a 
natural  network  model  definition  and  ensuing  construct.  At 
first  analysis,  it  appears  that  four  user  oriented  processing 
functions  may  be  considered: 

1.  Front-end  Interface  Processing  - That  processing 
associated  with  driving  and  interfacing  terminals 
similar  to  the  TIPs  of  the  AHP*  Network  [Roberts  73] 
but  without  their  conmunica ticn  network  functions  and 
access  control  for  authorized  users  and  levels  of 
usage. 

2.  Data  Base  Hanagement  Processing  - That  processing 
concerned  with  the  update,  retrieval,  and 
creation/aaintenance  of  files  comprising  the  data  base 
as  would  be  associated  with  the  data  base  machine 
concept  ([ACB  77],  [Canaday  74]), 

3.  Communications  processing  - That  processing  required  to 
act  as  a node  within  a network  or  the  interconnection 
point  for  the  passing  of  data  traffic  to  other 
functional  processors  similar  to  the  concept  of  the 
INPs  of  the  ABPA  Network  [Roberts  73]. 

4.  Principal  Element  of  Processing  - That  processing 
normally  associated  with  the  ability  to  do  high  speed 
arithmetic,  "number  crunching"  operations. 
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The  interconnecticn  of  these  user  processing  functions  need 
not  he  dependent  on  any  particular  geographical/ locality 
constraints.  Geographical  separation  or  co-location  and 
corahinaticns  thereof  for  the  four  processing  functions  could  be 
determined  by  operational  constraints  such  that  any  of  the 
following  may  arise: 

1.  All  functions  may  be  in  separate  locaticns. 

2.  Where  functions  are  co-located,  the  transmission  medium 
may  be  high  speed  parallel  transfers  such  as  exist 
within  a classical  multi/parallel  processing 
environment. 

3.  When  all  four  functions  are  physically  co-located,  then 
this  organization  may  be  viewed  as  a user  function 
oriented  multiprocessing  system.  One  might  view  this 
in  the  sense  of  a loosely  ccupled  system  (LCS)  or  a 
local  network. 

4.  Numbers  and  capabilities  cf  processing  elements  would 
vary.  Selection  of  numbers  and  capabilities  is 
dependent  to  a large  degree  on  operational  requirements 
and  other  factors  such  as; 

a)  Fault  tolerance. 

b)  Pile  allocation. 

c)  Sxisting  communications  facilities. 

d)  Existing  organizational  units,  e.g.  headquarters. 

e)  Hierarchical  structure  of  the  supported 
organization,  e.g.  chain  of  command  or  updating 
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regairenents. 

Although  specific  functions  have  been  identified  in  this 
basic  network  discussion,  it  is  obvious  that  a particular 
architect  night  desire  to  identify  other  types  of  processing 
functions  to  describe  an  entirely  different  network.  It  is 

desirable,  therefore,  to  develop  a model  which  provides  for  ' 

! 

generality  in  the  definition  of  the  network  topology  and  ; 
processing  functions.  What  then  beccmes  important  is  the  | 
network  system  processing  parameters,  topology  and  relative  i 
performance.  I 

Conseguently , such  a generality  will  provide  the  desired 
objective  of  using  the  model  to  perform  comparative  tradeoff 
analysis  of  competitive  design  approaches  and  will  additionally 
provide  a basic  tool  for  design  and  synthesis  of  a network  to 
meet  specific  processing  reguirements. 


i 


I 
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Chapter  II;  EXPONENTIAL  QDEDING  NETWOSK  (HCN)  MODELS 


In  this  chapter  we  will  look  at  modelling  of  computer 
systems  as  they  pertain  to  computer  networks.  In  particular  we 
shall  examine  the  use  of  network  queuing  theory  as  a basis  for 
an  analytic  model  of  a computet  network  composed  of 
functionalized  processors.  At  this  point  one  should  keep  in 
mind  Frank's  [1975]  observation  that: 


Modelling,  analysis  and  design  techniques  have 
made  substantial  progress,  but  the  best  results  are 
still  obtained  by  the  artistic  blend  of  theory, 
experience,  and  pragmatism. 


In  other  words,  it  should  not  he  expected  that  the  development 
of  the  models  about  to  be  presented  will  solve  all  anticipated 
problems.  It  will,  however,  be  shown  that  the  models  provide  a 
very  useful  tool  to  support  the  artistic  blend  needed  in 
computer  network  design  and  analysis. 
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Model  Objectives 

First,  we  should  note  that  there  ate  various  objectives  for 
Bodelling  which  will  influence  any  deternination  of  a nodel. 
Basically,  one  may  model  an  existing  computer  network  or  system 
to  ascertain  or  obtain  a projected  behavior  under  different 
workloads  or  configurations.  On  the  ether  hand,  one  might 
desire  to  model  a projected  system  not  knowing  all  of  its 
detailed  and  specific  components.  it  is  this  latter  use  of 
modelling  for  which  our  objective  is  set  - to  model  a projected 
system  not  knowing  all  of  its  detailed  and  specific  components. 
Consequently,  it  is  important  to  he  able  to  study  the  projected 
network  system's  behavior  in  order  to  gain  insight  and 
conceptual  clarity  into  the  overall  network  performance  rather 
than  perform  a detailed  examination  of  its  operations.  In  other 
words,  a global  approach  to  the  network  is  being  made.  As  a 
result  we  expect  to  be  able  to  accomplish  three  things: 

1.  Begin  with  a given,  projected  configuration  and 

topology  and  then  ascertain  whether  the  projected 
network  performs  according  to  a desired  set  of  criteria 
(trade-off  analysis) . 

2.  Begin  with  a desired  set  of  criteria  and  then  design  a 
network  to  meet  these  criteria  or  performance 
specifications  (design).  Obviously,  a practical  mix  of 
iterations  between  the  twe  capabilities  must  be 
employed  to  arrive  at  an  acceptable  network  solution. 
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3.  Make  the  inodel  natheniatically  tractable  in  order  tc 
petfori  the  ccmputaticns  rapidly  and  inexpensively. 
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' Sinulation  and  Analytic  Models 

I 

Throughout  this  presentation  it  beccmes  apparent  that  the 
intention  has  been  to  work  with  exponential  queuing  theory.  Why 
not  use  conventional  simulation?  Or  have  we  prejudiced  the 
choice  towards  the  analytic  model  and  exponential  queuing 
theory?  The  most  common  types  cf  convertional  simulation  deal 
either  with  event-oriented  models  to  represent  the  actual 
operations  of  a system  as  they  occur  event  by  event  or  through 
the  use  of  empirically  derived  data  represented  in  a 
deterministic,  analytic  form.  As  a consequence,  the  inherent 
criterion  for  conventional  simulation  is  the  desire  to  obtain 
precise  knowledge.  To  obtain  precise  knowledge,  however, 
implies  that  detailed  information  concerning  a computer  system 
to  be  modelled  is  known.  Such  is  not  always  the  case.  This 
deficiency,  however,  can  be  tolerated  if  we  are  willing  to  limit 
our  precision  of  the  model  as  we  have  already  stated.  The  type 
model  which  appears  to  work  well  in  this  situation  is  an 
analytic  model  based  on  stochastic  processes  and  exponential 
queuing  theory.  Such  a model  is  often  termed  an  Exponential 
Queuing  Network  (EQN)  model.  It  is  very  well  established 
([Basket  and  Huntx  73],  [Bhandiwad  and  Williams  7«],  [Browne  et 
al  75],  [Buzen  75],  [ Buzen  76a],  [Fuller  75],  [Giammo  76], 

I [Moore  71])  that  the  use  of  exponential  queuing  networks 

j satisfies  the  following  criteria: 

1.  Is  simple  and  inexpensive  to  use. 

I 


f 
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2.  Say  ieal  with  situations  where  only  partial  workload 
knowledge  exists. 

3.  Say  be  used  over  a wide  range  of  parameters. 

u.  Is  particularly  well  suited  for  the  isolation  of 
parameters. 

5.  Accepts  small  changes  in  input  values  and  does  not  lead 
to  large  changes  in  predicted  performance. 

6.  Is  several  orders  of  magnitude  more  cost  effective  than 
conventional  simulation  models. 

7.  Can  provide  a range  of  quick  answers  useful  in  planning 
evaluations. 

8.  Can  provide  an  insight  into  the  network  in  terms  of 
gross  computer  performance  measures  such  as  thruput  and 
response  time.  In  other  tords,  gualitative  or 
comparative  quantitative  and  approximate  answers  are 
provided . 

9.  Predicts  with  some  certainty  the  network  performance. 

19.  Permits  the  follow-on  development  of  arbitrary  networks 

of  the  models  under  general  conditions. 

It  is  this  set  of  EQS  criteria  which  meshes  with  the  desired 
model  objectives  and  is  therefore  selected.  In  particular,  it 
has  been  said  that  "as  a mathematical  approach,  probability 
models  are  generally  more  realistic  than  deterministic  models 
because  they  can  represent  the  irregular  and  unpredictable 
demand  by  computer  users."  [Coffman  and  Cenning  73] 
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Application  of  ^QN  rcdels 


The  validity  and  usefullness  of  the  EQN  acdels  has  been 
well  established  in  practice  and  in  publications.  However,  the 
models  employed  to  date  seem  to  he  concentrated  either  in  the 
area  of  analyzing  the  communications  aspect  of  a computer 
network  {[Kleinrock  74],  [FCleinrock  76]),  the  evaluation  of  a 
local  computer  system  ([  Buzen  71],  [Rose  76])  or  locally 
interconnected  systems  [Browne  et  al  75],  The  approach  to  be 
taken  here  is  to  examine  the  totality  of  the  computer  network 
from  the  standpoint  of  its  complete  processing/compu tational 
ability.  This  means  that  the  communications  aspect  is  just 
another  subsystem  {or  subnetwork)  of  the  total  computer  network 
and  must  be  parameterized  as  part  of  the  model.  Once  this  is 
accomplished,  there  is  complete  similarity  in  viewing  the  model 
in  its  local  or  geographical  separation.  At  this  point,  it 
becomes  valuable  to  analyze  the  system  from  a hierarchical  point 
of  view  to  obtain  values  of  specified  parameters  just  as  Browne 
et  al  [1'?75]  have  done  for  a specific  system.  This  is  the 
employment  of  the  technique  of  the  micro  to  macro  evolution  or 
of  first  examining  subsystems  of  the  computer  network  to  obtain 
specified  parameters  and  then  using  these  parameters  in  the 
higher  level  model.  In  this  way,  some  details  cf  lower  level 
models  may  be  incorporated  into  the  higher  level  model  ([Basket 
T,  «untz  73],  [Browne  et  al  76],  [Buzen  71],  [Colella  et  al  74]). 
The  approach  taken  in  this  research  is  to  employ 
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functionalization  as  a aeans  of  incorporating  the  lower  level 
parameters  into  the  higher  level  network. 


One  common  criticism  of  the  fCN  models  is  that  they  do  not 
seem  to  relate  to  the  actual  operatior  of  computer  networks  or 
systems.  As  we  shall  see  in  the  next  chapter,  certain 
assumptions  regarding  exponential  distributions  and  independence 
of  servers  and  arrivals  ate  known  to  be  false.  Yet  this  type  of 
model,  as  noted  above,  has  been  shown  to  demonstrate  very  good 
results.  Op  to  this  point,  most  analysts  have  accepted  the  ECU 
models  as  tools  and  are  extremely  careful  in  the  model 
employment  because  of  the  contradictions  that  exist  between  the 
necessary  assumptions  and  the  known  operational  variables  in  the 
systems.  In  1971  Buzen  [1971]  pointed  out  that  the  EQH  provided 
an  entirely  accurate  representation,  considering  that  only  the 
most  significant  aspects  were  included.  As  a consequence,  the 
model  should  not  be  judged  on  the  goodness  of  fit  of  the 
exponential  assumptions  hut  rather  on  the  validity  and  utility 
of  the  insights  which  are  ultimately  gained.  However,  Buzen  has 
recently  demonstrated  ([  1976a],  [ 1976t])  that  perhaps  these 
results  are  not  accidental  after  all.  Through  the  use  of  the 
concept  of  an  "Operational  Hethod"  of  computer  system  analysis 
based  on  a set  of  concepts  that  correspond  naturally  and 
directly  to  observed  properties  of  real  computer  systems,  the 
same  analytical  results  can  be  obtained  as  with  the  EQK . These 
results  are  encouraging  since  they  ace  further  indication,  in 


addition  to  the  excellent  results  obtained  in  practice,  that  the 
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approach  to  fce  taken  with  the  use  of  5CN  Bodelling  is  in  fact  a 
qcoi  one  and  say  he  accepted  even  though  some  of  the  assumptions 
nay  not  he  satisfied. 

For  the  application  of  the  model,  it  will  be  assumed  that 
the  computer  network  to  be  modelled  is  operating  under  peak 
workload  conditions  in  a steady-state  environment  over  a 
particular  interval  of  interest.  This  is  a reasonable 
assumption  since  one  would  expect  that  in  a Command  and  Control 
situation  the  most  critical  possibility  must  be  considered  and 
designed  for  or  evaluated.  Furtheraote,  it  cannot  be 
overemphasized  that  our  interest  lies  iu  gaining  an  insight  into 
performance  of  a £tojected  s^s^em.  One  cannot  expect  to  obtain 
results  that  are  mote  refined  than  tte  parameters  provided  to 
the  model.  Such  is  the  case  for  a projected  system.  It  must  he 
recognized  that  the  input  parameters  will  he  projections  so  we 
should  expect  only  projections  or  insights  as  results. 
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Chapter  TII:  NETWORK  CEFINITICN/STRUCTORE 


General  Approach 


For  discussion  purposes  a queuing  netxork  will  te 
established  to  represent  a computer  network  architecture 
consisting  of  the  four  functions  as  described  in  Chapter  I 
(Front-end  Interface  Processing  (FIP) , Communications  Processing 
(COflP) , Data-Base  Management  Processing  (DBH) , and  Principle 
Element  of  Processing  (PEP)) . The  queuing  network  model  is  not 
necessarily  restricted  to  these  four  processing  functions  but 
may  represent  any  functions  as  defined  by  the  architect.  Only 
the  complexity  of  computation  will  be  a restriction.  This 
complexity  is  a function  of  the  number  cf  nodes  (or  processors)  , 
number  of  jobs  and  number  of  job  types.  As  will  be  seen,  this 
complexity  is  not  a limiting  factor  for  practical  models  even 
for  large  networks. 

Procedures  are  developed  in  order  that  the  following  may  be 
accompl ished : 

1.  Entry  to  the  processing  network  through  the  FIP  by 
terminal  users. 

2.  Specify  various  routes  for  different  job  types  to 
represent  workload  criteria  in  terms  of  functionalized 
processor  requirements  and  assimilate  them  into  the 
network  representation. 
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3.  Determine  a response  time  expression  used  for 

evaluation  of  the  network  capacity. 

4.  Select  processor  service  rates  as  a design  procedure  to 
obtain  balanced  utilization  of  the  functions. 

With  these  procedures  we  are  able  to  accoanodate  the  two 
problems  of  analysis  and  design.  Three  different  approaches  for 
the  assimilation  or  representation  of  job  types  into  the  model 
will  be  obtained.  As  a result,  three  different  models  will  be 
derived  each  having  advantages  and  disadvantages.  All  three 
will  be  presented  in  this  work. 


I 
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Establishsent  of 

It  is  well  Xnown  that  the  following  are  inportant  in  a 
queuing  aodel  [Fuller  15]: 

1.  Arrival  process 

2.  Service  aechanisD 

3.  Queuing  discipline 

4.  Pouting 

Scheaat ica lly  cne  denotes  these  factors  as  shown  in  Figure 


arrival 
of 

custoners 


Queue 


departure 

of 

custcoers 


Server 


Figure  I1I-1 

A detailed  explanation  of  these  factors  is  contained  in  Fuller 
[1975]  and  Kleinrock  [1975]  and  will  not  be  repeated  here. 
Furtheraore,  an  extensive  discussion  of  these  factors  in  the 
foraulation  of  specific  aodels  for  the  evaluation  of  quantua 
controlled  service  disciplines  plus  an  extensive  tabulation  of 
research  and  survey  papers  in  this  area  are  included  in  Buzen 
[1972].  For  our  purposes  the  diagraaaatic  reference  above  aay 
be  considered  as  a node  in  a queuing  network  with  a specified 
probability  p(i,j)  that  a custoaer  exiting  a server  of  node  i 
will  proceed  to  node  j as  represented  in  Figure  III-2.  (Note: 
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In  Chapter  I a definition  was  given  for  a computer  network  and 
we  are  now  speaking  of  a queuing  network  which  is  a model  with 
more  than  one  node  and  will  be  used  to  represent  the  computer 
network . ) 


The  analysis  of  such  a network  of  nodes  (cr  queues)  has 
been  established  over  a rumber  of  years,  culminating  in  specific 
algorithms  for  solution.  Burke's  Theorem  ([Burke  56], 
[Kleinrcck  7*5])  established  that  for  an  l/M/m  queue  in  steady 
state  condition  with  a Poisson  input  process  (see  Kleinrock 
ri'>75]  for  discussion  of  system  categorization  such  as  .T/;i/a, 
which  represents  a Poisson  input  process,  an  exponential  service 
distribution,  and  a servers)  the  output  is  also  Poisson  and  is 
independent  of  the  other  processes  in  the  system.  Jackson 
[1956]  further  established  that  for  general  networks  of  nodes, 
even  with  feedback,  each  node  in  the  network  behaves  as  though 
it  were  an  S/Jl/a  system  with  Poisson  input.  This  is  so  even 
though  the  total  input  to  a node  is,  in  fact,  not  a Poisson 
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process. 

Subsequently,  Gordon  and  Newell  [ 1967  ] examined  a 
particular  subcase  of  Jackson  [1957],  They  have  considered 
queuing  networks  with  completely  general  interconnecting  paths 
as  did  Jackson  but  with  the  exception  that  only  closed  systems 
were  examined.  A closed  system  is  cne  where  there  is  no 
possibility  of  outside  arrivals  or  departures  of  customers  from 
the  system.  The  closed  queuing  network  systems  are  shown  to  be 
stochastically  equivalent  to  open  systems  (those  of  Jackson)  in 
which  the  number  of  customers  is  fixed  to  a specified  value,  say 
N. 

The  general  solution  for  these  networks  is  in  the  form  of 
the  probability  of  finding  a specified  number  of  customers  at 
each  of  the  nodes  (queue/server  pair)  in  the  queuing  network. 
Eased  on  this  solution,  other  performance  measures  of  the 


network  may  be 

obtained 

to 

include 

the 

response  time. 

utilization  and 

thruput. 

The 

concept 

of 

local 

balance  was 

introduced  by  Chandy  [ 1972] 

to 

provide 

an 

easier 

method  of 

solution.  local  balance  equates  the  rate  at  which  customers 
enter  and  leave  a single  queue  of  a network,  whereas  global 
balance  equates  the  rate  at  which  a Harkov  process  enters  a 
state  to  the  rate  at  which  the  process  leaves  the  state.  The 
work  of  Buxen  [1973]  and  Chandy's  study  of  the  local  balance 
technique  are  responsible  for  an  efficient  method  of  computing 
the  performance  measures  of  an  EQN. 
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The  aatheaatical  tasis  for  the  classical  queuing  theory  is 
quite  well  enumerated  in  the  references  cited  in  the  previous 
paragraphs.  However,  as  a suminary,  it  should  be  noted  that  the 
following  assumptions  have  been  made  in  order  to  provide  a 
mathematically  tractable  solution  [Buzen  S Denning  77]: 

1.  Stationary  stochastic  processes. 

2.  Stochastic  independence  of  jobs. 

3.  Job  transitions  follow  a Barkov  chain. 

u.  Stochastic  equilibrium  has  been  achieved. 

S.  Service  time  at  each  node  follows  an  exponential 
distribution. 

b . Ergcdicity. 

Although  this  set  of  assumptions  may  not  agree  with  the  computer 
architect's  concept  of  the  actual  computer  system's  operation, 
nor  are  the  assumptions  generally  and  readily  understood  by  the 
computer  system  manager  or  the  decision  maker  for  system 
selection,  the  fact  remains  that  the  use  of  queuing  models 
produces  excellent  results  for  performance  analysis. 
Consequently,  Suzen  R Denning  [1977]  have  proposed  a different 
approach  to  queuing  networks  based  on  the  concept  of 
"operational  analysis"  [ Euzen  76].  The  tasic  premises  employed 
are; 

1.  Finite  observation  periods.  No  assumption  of  steady 
state  is  made. 

2.  Work  conservation.  The  number  of  entries  to  a given 
state  must  equal  the  number  of  exits  from  the  state. 
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3.  Houcgene ity.  The  output  rate  at  any  node  is  dependent 
only  on  the  nunher  of  jots  at  that  node. 

This  operational  analysis  approach  yields  the  sane  solution 
forras  as  does  the  classical  queuing  theory  ([Buzen  75],  [ Buzen 
76a]).  It  would  appear,  however,  that  there  is  validity  to  the 
claia  that  this  latter  approach  provides  a higher  correlation  to 
actual  computer  system  operation.  Hence,  it  helps  to  make  the 
modelling  techniques  employing  queuing  theory  more  readily 
acceptable. 

A further  extension  to  the  theory  which  permits  service 
time  distributions  other  than  exponential  has  been  made  by 
Chandy  [ 1972].  He  shows  that  when  a processor  sharing  (PS)  or 
last-coae-f irst-served  preemptive-resume  (ICISPE)  discipline  is 
assumed,  the  only  releva 't  parameter  in  a service  distribution 
with  a rational  laplace  transform  is  the  mean  value. 
Furthermore,  a queue  with  either  of  the  two  conditions  stated 
behaves  the  same  way  as  a queue  with  a first-ccme-f irst-served 
discipline  and  exponential  service  distribution  with  the  same 
mean.  As  a consequence,  we  may  consider  distributions  other 
than  exponential  when  we  ace  interested  in  only  the  mean  values. 
This  fact  allows  more  flexibility  in  network  representation  in 
the  model.  Cn  the  other  hand,  it  also  demonstrates  the 
robustness  of  the  queuing  modelling  technique  since  the  solution 
forms  have  not  changed.  The  modelling  done  in  this  work  is 
restricted  to  the  set  of  service  disciplines  and  distributions 
just  discussed. 
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ASQ  (Arithmetic  Solaticn  to  Queues)  as  a Tool 

The  solution  methods  for  queuing  networlcs  provide  closed 
form  soluticns  but  require  considerable  computation  even  with 
the  techniques  employed  for  the  reduction  of  the  computations. 
As  a result,  specialized  computet  prcctams  have  been  developed 
for  their  implementation  such  as  that  of  Moore  [ 1971],  the  ASQ 
program  ([chandy  72],  [Asplund  76])  and  the  QNET4  program 
[Peiser  76],  Each  of  these  programs  has  specific  limitations, 
""he  Moore  [1971]  program  models  the  specific  system  then 
employed  at  the  tJniversity  of  Michigan  and  accommodates  one  job 
type  only.  The  Chandy  ASQ  program  has  been  extended 
rinformation  Research  Associates  75]  to  accommodate  up  to  three 
iob  types  and  a total  of  20  jobs  and  to  operate  in  a general 
purpose  mode.  It  is  implemented  in  Fortran.  The  Duke  version 
of  ASQ  [Asplund  76],  also  in  Fortran,  is  a general  purpose 
program  model  to  accommodate  up  to  99  jots  of  a single  type. 
Finally,  the  QNETU  program,  written  in  API,  is  also  of  a general 
purpose  nature  with  a limitation  of  two  job  types.  Each  of  the 
programs  mentioned  (except  Moore's)  is  proprietary  and  thus  not 
available  for  general  use. 

One  of  the  general  program  implementations  using  the 
principles  and  efficient  computational  methods  discussed  in  the 
previous  section  has  been  developed  by  Chandy  et  al  [1972]  and 
further  modified  to  operate  in  an  interactive  mode  by  Asplund 
[1976].  For  the  basic  model  development  in  Chapter  IV,  the 
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latter  interactive  iaplemen tation  at  Duke  University  known  as 
ASQ  (Arithmetic  Solution  to  Queues)  will  be  employed  as  a 
computational  tool. 


Duke  ASQ  provides  a number  of  perfc 
queue  length  at  a node,  utilization,  thr 
the  network  being  in  a specified  state. 
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ptors.  Then  a procedure 
type  will  be  developed. 


Two  further  model  representations  will  be  developed  which 
do  not  rely  entirely  on  the  Duke  ASQ  and  the  use  of  single  job 
types.  Pach  of  these  models  is  developed  with  a set  of  programs 
which  lend  themselves  to  use  on  a general  purpose  basis.  They 
have  the  ability  to  accommodate  multiple  job  types  and  to 
compute  response  time  directly  for  interactive  computer  network 
systems. 
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Network  Cescripticn 
Topology  and  Parameters 

A construct  for  network  definition  is  presented  and  will  be 
used  throughout  for  model  determination  and  evaluation-  The 
four  functions  of  the  selected  network  can  he  viewed  to  he 
topologically  interconnected  as  in  Figure  III-3. 


PIP  (or  F)  denotes  Front-end  Interface  Processing 
CCflP  (or  C)  denotes  CCHmunicalions  Processing 
DBM  (or  D)  denotes  Data  Ease  Management  Processing 
PEP  (or  P)  denotes  Principal  Element  of  Processing 

Figure  III-3 

Each  of  the  nodes  represents  a gueue/server  pair  with  parameters 
for  the  type  of  processing  function  involved.  In  order  to  have 
an  operational  concept  for  this  representation  of  a computer 
network  as  a queuing  network,  a projected  and  hopefully  logical 
interpretation  of  processing  functions  was  presented  in  Chapter 
I.  At  this  time  it  should  be  recalled  that  the  system  architect 


or  designer  must  go  through  similar  analysis  for  the  functional 
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representation  of  a conputer  network  to 
types  of  processing  functions  selected 
on  an  anticipated  trend  towards  the  fc 
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and  will  be  connected 


represented  as  a single 
to  the  FIP  of  the  Conputer 
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Processing  Network.  The  TEP.*?  node  will  consist  of  N servers 
representing  N users.  N is  then  indicative  of  the  number  of 
users  active  in  the  entire  system  (the  terminals  plus  the 
computer  Processing  Network,  CFN)  and  will  be  a major  factor  in 
the  design/evaluation  process. 

The  network  is  formally  represented  as  a quintuple 

QN  ::=  <a,p,s,D,R>  (III-I) 

which  is  specified  by: 

1)  A finite  set  H ::=  {n  (T) , n ( F)  ,n  (C)  , n (D)  ,n  ( P) , . . . } of 
nodes  in  the  network  such  that 


n(T) 

::=  TERB 

n (D) 

::=  CBB 

n(F) 

::=  FTP 

n (P) 

: :=  PEP 

n(C) 

::=  CCBP 

2)  A matrix  P ::=  [p{irj)]  fcr  all  n(i),n(j)  in  H,  where 
p(i,j)  is  the  probability  that  a job  leaving  node  n(i) 
will  proceed  directly  to  node  n(j). 

3)  An  ordered  set  S {|s(i)|}  for  all  n (i)  in  H,  where 

s(i)  = <u(i),a(i)>  such  that  u(i)  ::=  service  rate  of 

each  server  and  a(i)  is  the  number  of  servers. 

u>  An  ordered  set  T5  : :=  (i)  1}  tor  all  n (i)  in  n,  where 

d(i)  in  (FCFS,PS,LCFSPP)  representing  the  type  of 
queuing  discipline  used: 

FCFS  - First  Come  First  Served 
PS  - Processor  Sharing 

LCFSPR  - Last  Come  First  Serve  Pre-emptive  Resume 


1 


I 
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^ 5)  A finite  set  H : where  r(k)  is  an  ordered 

j pair  coaposed  of 

a)  An  ordered  tuple  <n  (T) ,n (F)  . ,n  (T) > traversing  n 

nodes  which  we  shall  call  a job  route  with  each 
route  representing  a job  type.  Another  view  of  the 
route  is  as  a circuit  in  the  network  graph 
initiating  at  n (T)  and  terninating  at  n (T)  with  no 
interaediate  occurrences  of  n (T) . In  our  case,  we 
will  consider  four  routes,  1 < k < 4.  In  theory 
there  is  just  the  practical  conputational  limit  on 
k. 

b)  The  frequency  of  occurrence  (normalized  to  1)  , 
q (k)  , of  route  r (k)  . 


i 
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Network  Cescripticn 

Job  Type  (Route)  Sub-Network  and  Terninal  Representation 


At  this  stage  there  are  two  further  items  that  require 
clarification  before  actual  i npleaentation  of  the  queuing 
network.  The  first  of  these  points  is  the  representation  of  the 
k (in  this  instance  k=4)  permissible  ]cb  types  or  routes  into 
the  single  network  and  the  second  is  the  situation  of  the  TERK 
node. 


In  order  to  represent  the  various  job  routes,  we  form  a k- 
th  subnetwork  for  each  distinct  r (k)  as  a quadruple.  (Note  that 
the  k corresponds  to  the  k-th  route  and  is  appended  as  a 
subscript  for  the  subnetwork  and  as  a superscript  for  the 
elements  of  the  quintuple.) 

SN(k)  ;:=  <Mk,Pk,Sk,Dk>  (III-2) 

which  is  specified  by: 

1)  A finite  set  nk  ::=  {n(i))  of  nodes  in  the  subnetwork. 
Note  that  the  following  necessary  condition  must  hold 
in  order  that  topological  integrity  be  maintained: 

Hk  ^ (the  first  element  cf  r(k))  ^ H 

(111-3) 

2)  A matrix  Fk  : :=  [pk(i,j)  ] for  n (i)  and  n ( j)  in  Nk, 


where  pk(i,j)  is  the  probability  that  a job  leaving 
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node  n (i)  will  proceed  directly  to  node  n(j)  for  the 
k-th  route. 

3)  An  ordered  set  Sk  ::=  ns(i)|}  for  all  n ( i)  in  Hk, 
where  s (i)  = <u{i),a(i)>  such  that  u (i)  service 

rate  of  each  server  and  a (i)  is  the  number  of  servers. 
Note  that  the  following  necessary  condition  holds: 


Sk  3 s 


(III-U) 


Furthermore,  the  service  rate  at  a processor  s(i)  is 
not  a function  of  the  route,  hut  is  constant  over  all 
routes.  We  shall  see  in  later  model  development  the 
relaxation  of  this  condition. 

U)  An  ordered  set  Ck  ::=  {|d(i)|}  for  all  n{i)  in  Hk, 
where  d(i)  is  in  {FCFS,PS,LCFSFE) . Again  note  that  the 
following  necessary  condition  holds: 


Dk  3 D 


(III-5) 


The  TERM  is  structured  such  that  s (T)  ::=  <u(T),N).  By  our 
terminology  1 is  the  mean  think  time  of  a user  at  a terminal  or, 
in  other  words,  the  mean  time  between  completion  of  a user's  job 
and  his  subsequent  initiation  of  another  job.  As  a result,  with 
i terminals  (0  < i < N)  in  the  thinking  mode,  the  effective 
service  rate  for  the  TFRH  node  is 


u{T)  = f(i)  = i / Z 


(III-6) 


Note  that  the  terminal  node  can  be  treated  as  an  Infinite 
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Service  (IS)  station  and,  therefore,  no  queuing  occurs  at  the 
node.  This  structure  for  TEEH  then  yields  the  desired  result 
and  is  the  saae  approach  as  that  established  in  Kleinrock  [ 1967, 
1976]  and  Suntz  and  Song  [1974], 


J 
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Chapter  IV:  BASIC  nCDBL  FCPMOLATICS  AND  IflPLEMENTATION 


Basic  nodel  Formulation 


Previous  to  this  chapter  a general  description  of  queuing 
theory  as  a modelling  technique  and  a means  of  providing  a 
network  description  has  been  presented.  Ne  are  now  ready  to 
establish  the  first  of  the  three  models  to  be  presented  which  we 
will  call  the  Basic  nodel  (BH) . The  model  approach  is  to  find  a 
means  of  consolidating  or  assimilating  the  k job  types  into  one. 
Consolidating  job  types  into  one  can  reduce  the  computational 
complexity  and  make  it  possible  to  use  the  available  ASQ  program 
[Asplund  76]  at  Duke  University  which  simplifies  the 
implementation  of  the  problem  solution.  Such  a procedure  is 
reported  in  leech  [1976].  The  consolidation  of  job  types  is 
done  by  assuming  that  superposition  of  the  k subnetworks  into 
one  queuing  network  may  be  accomplished. 


The  procedure  for  assimilation  of  the  k routes  and  their 


respective 

k 

subnetworks 

into  one 

queu ing 

network 

is 

accompl ished 

by 

noting  from 

equations 

III-1  and 

II  I- 3 

the 

elements 

s. 

and  D may  be 

formed  by  the 

ordered  me 

rging  of 

Hk 

for  all  k into  H 

, Sk  for  all 

k into  S and 

Dk  for  all 

k into 

D. 

This  merging  implies  the  restriction  that  s |i)  and  d(i)  are 
identical  for  all  routes  regardless  of  their  occurrence  in  a 
particular  route.  P,  must  be  specified  by  ascertaining  the  r (k) 
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7 'c.  what  remains  to  te  determined  in  equation  III 
probability  matrix  P.  First  however,  PJ;  V k must  be 
The  method  for  determining  Pk  is  applicable  not  on 
basic  model  but  will  be  used  for  succeeding  model 
references.  First  it  is  noted  that  the  information 
route  is  in  the  form  of  a deterministic  description, 
obtain  the  subnetwork  transition  probability  matrix 
description.  The  following  technique  can  be  us 
purpose. 

Define  tk{i,j)  to  be  the  number  of  <n(i),n(j) 
pairs  occurring  in  the  first  element  of  r(k).  The 
pairs  define  the  arcs  in  the  route.  Next  ascertain 
of  times  that  a job  traverses  or  visits  a node  vk(i) 

vk  (i)  = tk{i,j) 

j 

Now  the  proportion  of  the  job  visits  that  leave  n (i) 
n{j)  is 

pk  (i,  j)  = tk  (i,  j)  / vk  (i) 

This  procedure  then  completely  defines  the  Pk  V k. 


In  order  to  evaluate  P 

from 

the 

Pk  ¥ k,  the 

elements  p(i,j)  of  £ are 

evaluate 

d based  on  the 

pk(i,j) , the  i-th  and  j-th 

element 

of  Pk,  which 

calculated,  conditioned  on 

the 

occ 

urrence  of  the 

r (k) . The  probability  that 

route 

^ 9 

r (k)  , occurs 

5 3 
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which  is 

(IV-1) 

to 

go  to 

(IV-2) 

individual 
probability 
has  been 
k-th  route, 
is  q (k)  . 
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Therefore, 

p{i,j)  = g(i)  pk(i»j)  ’♦q(V),  {IV-3) 

k 

for  n (i)  ,n  (-j)  in  {n,n k) 

(Recall  that  the  appendage  of  k denotes  a superscript  over  which 

the  susaatioD  occurs.) 

where 

g{i)  = 1/(5^  I(i,k)  ♦ q()c))  (IV-4) 

k 

and 

I(i,)c)  = if  n(i)  is  in  r{k)  (IV-5) 

/O  otherwise 


The  assumption  made  here  is  that  the  q(k),  defined  as  the 
likelihood  of  the  occurrence  of  the  k-th  route,  is  also  the 
likelihood  for  the  occurrence  of  a job  at  the  k-th  route  at  a 
particular  node  (normalized  based  on  the  node  being  in  the  k-th 
route) . This  is  a reasonable  assumption  if  the  amount  of  time 
that  a job  remains  in  the  computer  processing  network  is  much 
greater  than  the  think  tine.  That  is  tc  say,  Tr  (N)  >>  Z.  Put 
in  other  words,  this  means  that  the  CPN  is  loaded  to  the  maximum 
and  nearly  no  jobs  are  thinking.  Such  an  approach  should 
provide  the  worst  case  situation  in  regards  to  response  tine  or 
an  upper  bound.  The  point  to  be  ascertained  is  what 
inaccuracies  are  incurred  and  to  what  extent.  Such  a comparison 
for  the  selected  network  will  be  presented  in  Chapter  V. 
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Response  Time  Expression 


As  discussed  earlier,  we  must  determine  a means  of 
evaluating  the  mean  response  tine,  Tr  (N)  , as  a function  of  the 
number  of  users  in  the  total  system.  The  response  time  is 
defined  as  the  time  that  a job  is  initiated  into  the  network 
until  the  job  is  returned  to  the  terminal  from  the  network. 
Eortunately,  an  extremely  useful  relationship  is  available  to 
us.  Little's  Result  [ Kleinrock  76]  states  that: 


J = L * T 


(IV-6) 


where  J ::=  number  of  jobs  in  the  computer  network  system 

L ::=  mean  arrival  rate  of  jobs  tc  the  computer  network 
system 

T ::=  mean  time  in  the  computer  network  system 
It  is  the  mean  time  in  the  system  that  we  have  defined  as  the 
response  time.  if  we  know  the  number  of  jobs  in  the  terminal 
node,  say  Q (T) , then  J = N - C(7).  Now  recall  the  general 
configuration  of  the  network  in  Figure  IV-1  (same  as  Figure 
III-U)  : 


Compu  ter 

Processing  Network 


(TFRB 


Q(T)  jobs 


J jobs 


Figure  lV-1 
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Hy  Little's  result,  equation  IV-6  becomes 

L'  = Q(T)  / Z {IV-7) 

which  is  the  input  to  the  terminal  node  coming  from  the  rest  of 
the  network  and 

L = J / Tr  (N)  (IV-8) 

which  is  the  output  from  the  terminal  node  going  to  the  rest  of 
the  network.  Also  since  equation  III-2  holds  such  that  no  queue 
is  allowed  to  form,  then  in  the  steady  state  condition 

L = L*  (17-9) 

Consequently,  combining  IV-7  and  IV-8  into  IV-9: 

Q(T)  / Z = J / Tr  (N) 

or  re-arranged 

Tr(N)  = (N  - Q(T))  Z / N (T)  (IV-1C) 

= [S  / Q(T)  - 1]  Z 

which  is  the  desired  result.  In  the  latter  expression  it  may  be 
noted  that  Tr  (N)  is  inversely  related  to  Q (T) , the  queue  length 
or  number  of  jobs  at  TE^n,  which  is  the  only  parameter  being 
varied  at  this  point.  This  condition  is  supported  by  intuitive 
reasoning.  The  fewer  the  number  of  jobs  at  the  terminal  (or  in 
a think  state)  the  greater  is  the  number  of  jobs  being 
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processed,  since  this  is  a closed  network.  With  a larger  namber 
of  jobs  in  the  processing  state  or  loading  the  Computer 
Processing  Network,  the  longer  it  can  be  expected  for  the  jobs 
to  be  processed. 


1 
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Parameter  Cefiniticn 


With  the  augmentation  of  the  Terminal  node  (TERK)  as  in 
Figure  III-4  the  network  previously  shown  in  Figure  IV-1  is 
shown  in  Figure  11-2. 


FTP 

(or 

denotes 

CCflP 

(or 

C) 

denotes 

DEM 

(or 

D) 

denotes 

PEP 

(or 

P) 

denotes 

denotes  Front-end  Interface  Processing 
denotes  CCMmunica tions  Processing 
denotes  Data  Base  Management  Processing 
Principal  Element  cf  Processing 


Figure  IV-2 


This  network  according  tc  our  structure  is  represented  as 


QN  ;;=  <5,P,S,D,R>  (IV-11) 

with  the  following  element  representations: 

1.  n ::=  {n  (T)  ,n  (F)  ,n  (C)  ,n  (D)  ,n  (P)  } {IV-12) 


where 

n(T)  = 

TERH 

n (C) 

::=  DBH 

n(F)  ::  = 

PIP 

n (P) 

: :=  PEP 

n(C)  ::  = 

COHP 

2.  The  matrix  P is  unknown  at  this  point  and  must  be 
calculated  based  on  the  Pk  as  described  earlier  in  this 
chapter. 

s • • = 

•J  » m 


3 


{|s  (T)  ,s  (F)  ,s(C)  ,s  (C)  ,s  (P)  n 


(IV-13) 
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which 

does  not 

provide  distinctions 

between  routes  for 

this 

model. 

The  number  of  servers  at  a node  and  the 

service 

rates 

of  these  servers 

are  specified  as 

folio 

ws: 

s(T) 

• • 

<1/2 

, s> 

(IV-14) 

s(F) 

: : = 

<1.5 

, 1> 

(IV-15) 

s(C) 

: : = 

<1.0 

, 1> 

(IV-16) 

s(D) 

: : = 

<0.  2 

, 1> 

(IV-17) 

S(P) 

: : = 

<0.2 

, 1> 

(IV-18) 

The  u (i)  have  been  initially  selected  as  though  these 
are  known  mean  service  rates  fcr  existing  or  projected 
processors  to  perform  each  of  the  functions, 
a.  D ;;=  {IPS,  FCFS,  FCFS,  PS,  PS|}  (IV-19) 

The  d (i)  are  selected  based  on  the  same  rationale  as 
the  u (i) . That  is,  the  selection  is  based  on  existing 
or  projected  processors  to  perform  each  of  the 
functions.  Again,  distinction  between  routes  is  not 
possible  for  the  basic  model. 

5.  P {r{1),  r|2),  r(3),  r (4)  } (IV-20) 

Pour  routes  have  been  selected  as  being  representative 
of  a possible  ordering  of  processing  functions  to  be 
performed  for  the  four  different  job  types: 
a.  r(1)  ::=  (<n  (T)  , n (?)  , n (T)  > ; 0.6)  {IV-21) 

This  route  may  be  visualived  as  one  job  type 
requiring  query  generation  processing  or  menu 
selection,  or  even  security  access.  In  other 
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words,  the  short  interactive  requests  that  would  be 
serviced  within  a "local"  geological  area  to  the 
FTP.  Me  expect  a large  percentage  of  the  total 
number  of  requests  to  be  of  this  type. 
Consequently,  q(1)  is  expected  to  be  rather  large. 


b.  r{2) 

::=  (< 

n(T) 

r n{F),  n(C), 

(IV- 

22) 

n(D) 

, n(C),  n(F), 

n(T)> 

0 

.2) 

This 

route 

is  representative 

of  a ( 

da 

ta 

base 

upd 

ate 

or 

retr ie 

val 

requirement 

with 

n 

o 

large 

scale 

computation 

involved  as  no  de 

imands 

a 

re 

made 

to 

the 

PEP. 

c.  r (3) 

{< 

n(T) 

, n (F)  , n (C), 

n(D), 

{IV- 

23) 

n(C) 

, n (P)  , n (C) , 

n (F)  , 

n 

(T)  > ; 0 

.1) 

This 

route 

is 

representative 

of 

a 

retrieval 

or 

update  followed  by  a subsequent  requirement  to 
perform  considerable  processing  prior  to  job  return 
to  the  user. 

d.  r(tt)  (<n(T),  n (F)  , n (C)  , (IV-2«) 

n(P)  , n |C)  , n(F),  n (T)  > ; 0.1) 

This  route  represents  a strictly  computational  job, 
i.e.  not  requiring  access  to  the  EBd. 

It  is  important  to  realize  at  this  point  that  the 
definition  of  routings  and  associated  probabilities  is  part  of 
the  parameter  selection  and  to  a great  extent  is  derived  from 
the  designers*  or  analysts'  artistic  blend  of  theory. 


: i 


:■  i 


p, 


I 


experience,  and  pragmatism  discussed  in  Chapter  I 


Values 
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selected  in  this  iirpleaentation  ate  representative  of 
configuration  and  are  in  no  way  intended  to  justify  a 
system . 


a possible 
particular 
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I 

I 


ASQ  lapletDentaticn  for  Trade-Cff  Analysis 


The  parameters  defined  for  the  implementation  of  the  Basic 
:iodel  yield: 


r 
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ASQ  Calculation  for  Trade-Cff  Analysis 


Once  the  paraaeter  definition  and  deteraination  have  been 
completed,  it  is  relatively  straight  fcrvard  to  introduce  these 
values  into  Duke  ASQ  following  the  ASQ  Users'  Manual  [Asplund 
76],  For  the  preceeding  values  the  ASQ  Network  Description  is 
at  Figure  IV-3.  It  should  be  noted  that  the  u (T)  is  calculated 
in  accordance  with  equation  III-2  and  is  based  on  Z = 15 
seconds.  Values  for  1 < N < 75  are  provided.  The  Z is  obtained 
based  on  the  reasoning  that  r(1)  jobs,  ty  virtue  of  their 
description,  should  have  a mean  think  time  equal  to  five  seconds 
as  where  the  other  jobs,  r (k)  , 2 < k < 4,  should  have  a mean 
think  time  equal  to  thirty  seconds.  Weighting  the  think  times 
in  accordance  with  their  respective  probabilities  of  occurrence, 
we  have  produced  a virtual  mean  think  time  as  follows: 


Z = ^ zk  q(k)  (IV-30) 
K 

where  zk  is  the  think  time  for  route  k.  Note  that  the 
assumption  of  maximum  loading  or  that  all  jobs  are  in  the  CPN  is 
also  made  here  as  was  done  for  determining  q (k) . 

At  Table  IV-1  are  the  computed  results  for  the  Basic  Model 
as  defined.  Shown  is  the  queue  length  at  n(T),  Q(T),  and  the 
respective  response  time  function  Tr(N)  in  accordance  with 
equation  (IV-10).  Additionally,  the  thruput  and  utilization  of 
each  functional  processing  node  may  be  easily  computed  using  ASQ 
(Tables  IV-2  and  IV-3,  respectively).  Curves  for  the  Tr(N)  and 
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utilization  for  the  processor  with  the  largest  values  (DBfJ)  are 
also  shown  in  Figure  lV-4.  This  data  thus  provides  inforaation 
regarding  the  overall  network  perfcraance  and  is  but  one 
iteration  of  an  evaluation  process.  It  is  evident  that 
additional  iterations  aay  be  perforaed  varying  the  different 
parameters  to  obtain  a relative  comparison  of  the  network 
performance  under  various  conditions. 

A noteworthy  point  at  this  juncture  is  the  utilization  for 
DBfl  of  Table  17-3.  For  this  particular  selection  of  the  u (i)  , 
the  DBJ1  acts  as  the  bottleneck  as  defined  in  Huntz  and  Hong 
[1974].  The  next  bottleneck  occurs  with  the  FIP  followed  by  the 
PEP  and  then  the  CCMP.  Cr,  in  ether  words,  we  can  increase  the 
u (F)  , u (P)  , and  u (C)  and  will  not  improve  the  response  time  or 
thruput  as  long  as  the  DEH  remains  the  bottleneck.  Jobs  just 
simply  cannot  proceed  through  the  network  any  faster.  This 
observation,  consequently,  further  illuminates  the  desire  for  a 
design/synthesis  procedure. 
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TFArr-CFF  AN  5LySTS 


ASO  VERSION  7 RELEASE  3 
DUKE  UNIVERSITY  JUNE  1976 

THE  NETUORK  DESCRIRTION  RESIDES  K ! <FILENAME>  ? 
TRADEOFF 

ENTER  NUMBER  OF  JOB  TYPES  (1  OR  MORE) 

1 

ENTER  NUMBER  OF  NODES  + BRANCHING  POINTS 
5 

ENTER  NUMBER  OF  NODES 
S 

MILL  YOU  HAVE  EXTERNAL  POISSON  SOURCES? 

NO 

ENTER  EDGES 


FIP 

0.20000 

COMP 

FIP 

0.80000 

TERM 

COMP 

0.45834 

FIP 

COMP 

0.33333 

DBM 

COMP 

0.20834 

PEP 

TERM 

1.00000 

FIP 

DBM 

1.00000 

COMP 

PEP 

1.00000 

COMP 

ENTER  RATES  FOR  NODES  WITH  FCFS  DISCIPLINE 

FIP  1.50 

COMP  1.00 

* 

ENTER  RATES  FOR  NODES  UITH  PS  OR  LCFSPR  DISCIPLINE 
DBM  0.20 

PEP  0.20 

TERM 


0.06667 

0.13333 

0.20000 

0.26667 

0.33333 

0.40000 

0.46667 

0 . 53333 

0.60000 

0 . 66667 

0.73333 

0.80000 

0.86667 

0.93333 

1.00000 

1.06667 

1 . 13333 

1.20000 

1.26667 

1.33333 

1.40000 

1.46667 

1.53333 

1.60000 

1.66667 

1.73333 

1.80000 

1.86667 

1.93333 

2.00000 

2.06667 

2.13333 

2.20000 

2.26667 

2.33333 

2.40000 

2.46667 

2.53333 

2.60000 

2.66667 

2.73333 

2.80000 

2.86667 

2.93333 

3.00000 

3.06667 

7.13333 

3.20000 

3.26667 

3.33333 

3.40000 

3.46667 

3.53333 

3.60000 

3 . 66667 

3.73333 

3.80000 

3.86667 

3.93333 

4.00000 

4.06667 

4.13333 

4.20000 

4.26667 

4.33333 

4.40000 

4.46667 

4.53333 

4.60000  • 

4.66667 

4.73333 

4.80000 

4.86667 

4.93333 

5.00000 

* 

ENTER  DEGREE(S)  OF  MULTIPROGRAMMING 
1 

THE  TEMPLATE  IS: 

FIP  COMP  DBM  PEP  TERM 

riiiire  IV-T 
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FSSPONSE  lines 

TRADE-OFF  ANALYSIS 


N 

Q(T) 

Tr  (N) 
|sec) 

1 

0. 84'^04 

2.86 

2 

1. 66889 

2.98 

3 

2. 48554 

3.  10 

4 

3. 28890 

3.24 

5 

4. 07782 

3.  39 

6 

4.85102 

3.55 

7 

5.60716 

3.73 

8 

6.  34480 

3.91 

9 

7.06241 

4.12 

10 

7.75838 

4.33 

1 1 

8. 43108 

4.57 

12 

9.07884 

4.83 

13 

9.69999 

5.  10 

14 

10.29290 

5.54 

15 

10.85605 

5.73 

16 

11.38808 

6.07 

17 

11.88782 

6.45 

18 

12. 35437 

6.85 

19 

12.78717 

7.29 

20 

13. 18600 

7.75 

21 

13.55106 

8.25 

22 

13. 88299 

8.77 

23 

14. 18279 

9.33 

24 

14. 45189 

9.91 

25 

14.69198 

10.52 

30 

15.53124 

13.97 

35 

15.96338 

17.89 

40 

16. 18793 

22.06 

45 

16.31142 

26.38 

50 

16. 38304 

•0.78 

55 

16. 42639 

25.  22 

60 

16. 45313 

39.70 

65 

16. 47008 

44.20 

70 

16. 48073 

48.71 

75 

16. 48756 

53.  23 

Table  IV-1 
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N 

FIP 

COMP 

DBM 

PEP' 

N 

.00000 

0.07001 

0.03055 

0.01018 

0.00636 

1.00000 

.00000 

0.13907 

0.06069 

0.02023 

0.01264 

2.00000 

.00000 

0.20713 

0.09038 

0.03013 

0.01883 

3.00000 

.00000 

0.27408 

0.11960 

0.03987 

0.02492 

4.00000 

. 00000 

0.33982 

0.14828 

0.04943 

0.03089 

5.00000 

.00000 

0.40425 

0.17640 

0.05880 

0 . 03675 

6.00000 

.00000 

0.46726 

0.20390 

0.06797 

0.04248 

7.00000 

.00000 

0.52873 

0.23072 

0.07691 

0.04807 

8.00000 

.00000 

0.58853 

0.25682 

0.08560 

0.05350 

9.00000 

. 00000 

0.64653 

0.28212 

0.09404 

0 . 05873 

10.00000 

.00000 

0.70259 

0.30659 

0.10219 

0.06387 

11.00000 

.00000 

0.75657 

0.33014  - 

0.11005 

0.06878 

12.00000 

.00000 

0.80833 

0.35273 

0.11757 

0.07349 

13.00000 

.00000 

0.85774 

0.37429 

0.12476 

0.07798 

14.00000 

, 00000 

0.90467 

0.39477 

0.13159 

0.03225  ! 

15.00000 

.00000 

0.94901 

0.41411 

0. 13804 

0.08628 

16.00000 

.00000 

0.99065 

0.43229 

0.14409 

0.09006 

17.00000 

.00000 

1 .02953 

0.44925 

0.14975 

0.09360 

18.00000 

.00000 

1.06560 

0.46499 

0.15499 

0.09688 

19.00000 

.00000 

1.09883 

0.47949 

0.15983 

0.09990 

20.00000 

.00000 

1.12925 

0.49277 

0.16425 

0.10266 

21.00000 

.00000 

1.15691 

0.50484 

0.16828 

0.10518 

22.00000 

.00000 

1.18190 

0.51574 

0.17191 

0.10745 

23.00000 

.00000 

1.20432 

0.52552 

0.17517 

0.10949 

24.00000 

.00000 

1.22433 

0.53425 

0.17808 

Q.11I3I 

25.00000 

C h 'T  D t-.  e r IV 


6.00000 

7.00000 

8.00000 

9.00000 

10.00000 


;io. 00000 

117.00000 

1 18.00000 

119.00000 

120.00000 


*30.00000 

35.00000 

Uo. 00000 

;45. 00000 
>50.00000 


t]  T ^ T r' 


^^^T’'CT3 


FIP 

0.04667 
0.09272 
0.13008 
0.18272 
0.22655 
0.26950 
0.31151 
0.35249 
0.39236 
0.43102 
0.46839 
0.50438 
0.53889 
0.57183 
0.6031 1 
0.63267 
0.66043 
0.68635 
0.71040 
0.73255 
0.75284 
0.77128 
0.78793 
0.80288 
Qt  81622 
0.86285 
0.88685 
0.89933 
0.90619 
0.91017 
0.91258 
0.91407 
0.91500 
0.91560 
0.91598 


COMP 

DBM 

PEP 

0.03055 

0.05091 

0.03182 

0.06069 

0.10114 

0.06322 

0.09038 

0.15064 

0.09415 

0.11960 

0.19933 

0.12459 

0.14829 

0.24714 

0.15447 

0.17640 

0.29400 

0.18376 

0.20390 

0.33983 

0.21240 

0.23072 

0.38453 

0.24034 

0.25682 

0.42802 

0.26753 

0.28213 

0.47020 

0.29389 

0.30659 

0.51097 

0.31937 

0.33014 

0.55023 

0.34391 

0.35273 

0.58788 

0.36744 

0.37429 

0.62381 

0.33990 

0.39477 

0.65794 

0.41123 

0.41412 

0.69019 

0.45139 

0.43229 

0.72047 

0.45031 

0.44925 

0.74875 

0.46799 

0.46499 

0.77498 

0.48438 

0.47949 

0- 79.915 

0.49949 

0.49277 

0.82127 

0.51332 

0.50484 

0.84139 

0.52589 

0.51574 

0.85956 

0.53725 

0.52553 

0.87587 

0.54744 

0.53426 

0.39042 

0.55654 

0.56478 

0.94128 

0.58833 

0.58049 

0.96747 

0.60470 

0 • 58866 

0.98109 

0.61321 

0.59315 

0.98857 

0.61788 

0.5957a 

0.99291 

0.62060 

0.59733 

0.99554 

0.62223 

0.59830 

0.99716 

0.62325 

0.59891 

0.99818 

0.62389 

0.59931 

0.99883 

0.62430 

0.59956 

0.99925 

0.62456 

Table  IV-t 


• N 
1 1 1.00000 
2.00000 

3.00000 

4.00000 

5.00000 

6.0000C 

7.00000 

8.00000 

9,00000 

10.00000 
■■  11.00000 
12.00000 

13.00000 

14.00000 

15.00000 
"■  16.00000 

17.00000  ! 

18.00000 

19.00000 

|20. 00000 

21.00000 
22.00000 

23.00000 

24.00000 
25 1 00000 

30.00000 

35.00000 

40.00000 

45.00000 

150.00000 

55.00000 

60.00000 

65.00000 

70.00000 

75.00000 


I Chapter  IV  7C 

^ Balanced  Processor  Design 

I 

I 

Now  suppose  we  would  like  to  have  the  nost  efficient 

* processor  utilization,  which  implies  an  optiiua  thrupot  rate,  or 

at  least  know  what  it  should  be  for  a specified  network.  We  do 
know  that  such  a network  should  te  a balanced  one  ([Koore  71], 
[Baekeura  and  Boyd  71]).  This  does  not  imply  that  what  results 
is  the  most  efficient  network  in  terms  of  other  parameters.  It 
does  indicate  where  efficiency  may  he  achieved.  Therefore,  we 
would  like  to  be  able  to  select  service  rates  to  ensure  a 
balanced  system  [ Buntz  and  Wong  7U].  in  terms  of  our  structure 
QN  = <B,  F,  S,  D,  R>  we  need  to  determine  S.  However,  let  os 

consider  the  case  of  a (i)  = 1 V i since  we  know  [Kleinrock  76] 

that  for  any  s(i)  = <u(i),a(i)>  we  can  obtain  an  equivalent 

s(i,eq)  = <u(i,eq),1>.  Therefore,  we  will  assume  a(i)  = 1 in 

this  discussion.  Furthermore,  we  can  now  proceed  in  the 

development  of  the  u(i,eq)  V n{i)  so  that  all  nodes  approach 

saturation  uniformly  and  no  one  noce  becomes  a bottleneck. 
Hence  we  will  have  a balanced  system.  Eased  on  the  work  of 
Buntz  and  Wong  [1974]  we  have  the  following  asymptotic 

properties  for  Tr  (N)  as  N becomes  large: 

Tr(N)  = N {v{s)  / u(s))  - Z (IV-31) 

Here  we  let 

I v{i)  = L(i)  / L(T)  2 < i < IM  (IV-32) 


X 
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I 


where  v (i)  is  the  average  number  of  visits  to  node  n(i)  per  user 
request  as  before  and  I (i)  is  the  mean  arrival  rate  of  jobs  at 
node  i.  T denotes  the  TEPfl  node  as  before.  The  subscript  s 
denotes  the  saturation  node  n (s) . This  asymptotic  condition  is 
shown  in  Figure  IV-5.  The  line  for  Tr(1)  is  also  shown  and  is 
given  by 

« 

Tr  (1)  = S (IV-33) 

k 


Figure  IV-5 

As  N becomes  large  then  from  equation  IV-31  the  slope  m of  the 
linearized  response  time  is  expressed  as 

m = V (s)  / u (s)  (IV-34) 

Since  saturation  is  being  approached  at  the  same  rate  for  all 
nodes,  in  a balanced  system  one  has 

v(i)  / u(i)  = m ¥ n(i),  i / T (IV-35) 

Osing  IV-35  in  IV-31  and  TV-33  we  have 


W 
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Tr  (N)  = N B - Z (IV-36) 

Tr(1)  = (im  - 1)  B (IV-37) 

equations  IV-36  and  IV-37  may  be  equated  in  order  to  solve  for 
N*,  the  point  of  intersection  of  the  twe  asymptotes  given  by 

N*  = {{\n\  - 1)  m ♦ Z)  / a (IV-38a) 

such  that 

0 = Z / (N*  - I B I ♦ 1)  (IV-38b) 

He  can  observe  that  IV-35  nay  be  rewritten  in  the  form 

u(i)  = v(i)  /mV  n(i),  i / T (IV-39) 

Thus,  if  we  are  given  the  acceptable  level  cf  response  time, 
Tr(l),  equation  IV-37  can  be  used  to  evaluate  m.  Alternatively, 
if  we  are  given  the  number  of  terminals  that  must  be  supported 
without  saturating  the  coaputer  network,  then  equation  IV-38t 
can  be  used  to  evaluate  a.  Consequently,  if  we  should  know 

V (i) , then  equation  IV-39  may  be  employed  to  provide  a means  of 

deriving  the  desired  service  rates,  u(i).  But  recall  that  IV-32 
does  establish  the  v (i) . He  aerely  need  to  solve  for  the  L(i) 
and  let  L (T)  = 1 and  the  desired  u (i)  are  determined.  To  solve 
for  L(i)  we  may  employ  the  balance  equations  of  Gordon  and 
Vcwell  [1967]  for  the  entire  network  such  that: 

v(i)  = v(j)  p(i,j)  n{i),n(j)  in  H 


I 


(IV-UC) 
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The  L is  not  uniquely  defined  but  by  fixing  the  L (T)  =1  for  all 
n (i)  where  i / T the  L can  then  be  uniquely  determined. 
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Synthesis  Procedures  fcr  Load  Balancing 


The  conpleaentary  capability  for  utilization  of  the  Basic 
Model,  as  previously  described,  is  that  of  a synthesis  procedure 
to  achieve  the  most  "efficient  processor  utilization."  We  have 
said  that  this  should  be  done  by  selecting  processor  service 
rates  so  that  all  processing  nodes  approach  saturation  uniforaly 
vith  no  one  node  becoaaing  a bottleneck. 

Eaploying  the  procedure  just  described,  we  first  proceed  tc 
obtain  the  v (i)  by  use  of  equation  IV-UO  with  v(T)  = 1: 

[ Vv  (?)  ,v  (C)  ,v  (D)  ,v  (P)  ] = [ 1 v (F)  V (C)  v(C)  v(P)  ] P 

(IV-41) 


where  P is  given  in  equation  TV-29.  At  A 
prograa  which  solves  the  v(i)  given 
program  uses  a routine  for  the  solution 
unknowns  of  Gallie  and  Paam  [1976]  and 
and  Sigmon  [1977],  Osing  the  Appendix  I 
following: 


ppendix  I is  a PL/I 
the  P as  input.  This 
of  n equations  in  N 
was  provided  by  Ballard 
prograa  we  obtain  the 


v(F)  = 1.2500 

V (C)  = 0.54546 

v(D)  = 0.181813 

v(P)  = 0.1  13638 


Assuaing  a saturation  of  N*  = 24  is  tolerable,  then  by  IV-38b  we 


have 
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a = ^ / (N*  - Ml  ♦ 1)  = 15  / (2U  - 5 ♦ 1) 

and  by  equation  IV-37  it  is  seen  that 

Tr(l)  = (ini  - 1)  = a (0.75)  = 3.0  seconds 

This  is  the  horizontal  asymptote  of  Figure  IV-5 
assumed  to  be  an  acceptable  value.  Consequently 


u (F) 

= 1.25  / 0.75  = 1. 

6 6 667 

(a) 

U(C) 

= 0.54546  / C.75  = 

0.72728 

(b) 

u (D) 

= 0. 181813  / 0.75  = 

0.  2424  1 

(c) 

u (P) 

= 0.113638  / 0.75 

= 0.15  152 

(d) 

= 0.75  (IV-U2) 

(IV-43) 

which  may  be 
the  u(i)  are: 

(IV-4U) 


The  network  description  fcr  the  balanced  implementation 
governed  by  equations  IV-11  and  IV-41  is  at  Figure  IV-a.  Table 
TV-4  lists  the  response  times  which  is  followed  by  the  thruput 
and  utilization  tables,  IV-5  and  IV-6,  respectively. 
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ASQ  VERSION  7 RELEASE  3 
tiUKE  UNIVERSITY  JUNE  1976 

THE  NETWORK  DESCRIPTION  RESIDES  » ! <F1LENA«E>  ? 
BALANCED 

ENTER  NUMBER  OF  JOB  TYPES  (1  OR  MORE) 

1 

ENTER  NUMBER  OF  NODES  + BRANCHING  POINTS 
5 

ENTER  NUMBER  OF  NODES 


WILL  YOU  HAVE  EXTERNAL  POISSON  SOURCES? 
NO 

ENTER  EDGES 


FIP 

0.20000 

COMP 

FIP 

0.80000 

TERM 

COMP 

0.45834 

FIP 

COMP 

0.33333 

DBM 

COMP 

0.20834 

PEP 

TERM 

1.00000 

FIP 

DBM 

1.00000 

COMP 

PEP 

1.00000 

COMP 

ENTER  RATES  FOR  NODES  WITH  FCFS  DISCIPLINE 

FIP  1.66667 

COMP  0.72728 

« 

ENTER  RATES  FOR  NODES  WITH  PS  OR  LCFSPR  DISCIPLINE 
DBM  0.24241 

PEP  0.1S152 

TERM 


0.06667 

0.13333 

0.20000 

0.26667 

0.33333 

0.40000 

0.46667 

0.53333 

0.60000 

0.66667 

0.73333 

O.COOOO 

0.86667 

0.93333 

1.00000 

1.06667 

1.13333 

1.20000 

1.26667 

1.33333 

1.40000 

1.46667 

1.53333 

1 .60000 

1.66667 

1.73333 

1.80000 

1.86667 

1.93333 

2.00000 

2.06667 

2.13333 

2.20000 

2.26667 

2.33333 

2.40000 

2.46667 

2.53333 

2.60000 

2.66667 

2.73333 

2.80000 

2.86667 

2.93333 

3.00000 

3 . 06667 

3. 13333 

3.20000 

3.26667 

3.33333 

3.40000 

3.46667 

3.53333 

3.60a00 

3.66667 

3.73333 

3.80000 

3.86667 

3.93333 

4.02000 

4.06667 

4.13333 

4.20000 

4.26667 

4.33333 

4.40000 

4.46667 

4.53333 

4.60000 

4 . 66667 

4.73333 

4.80000 

4.86667 

4.93333 

5.00000 

* 

ENTER  DEGREECS)  OF  MULTIPROGRAMMING 
1 


THE  TEMPLATE  IS: 

FIP  COMP  DBM  PEP  TERM 


Yiqnrp  IV- F 
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RESPONSE 

TIMES 

BALflNCEC 

SYSTEM 

N 

Q(T) 

Tr  (N) 
(sec) 

1 

C. 83332 

3.00 

2 

1. 65517 

3.  12 

3 

2. 46459 

3.  26 

4 

3.26058 

3.  40 

5 

4. 04212 

3.55 

6 

4.80809 

3.72 

7 

5. 55732 

3.89 

8 

6. 28863 

4.08 

9 

7. 00077 

4.28 

10 

7.69251 

4.50 

11 

8. 36261 

4.73 

12 

9.00986 

4.98 

13 

9.63312 

5.24 

1 4 

10. 23131 

5.33 

15 

10. 80349 

5.83 

16 

11 . 34888 

6.  15 

17 

11. 86685 

6.49 

18 

12.35701 

6.85 

19 

12. 81916 

7.23 

20 

13.  25339 

7.64 

21 

13.  65998 

8.06 

22 

14. 03950 

8.51 

23 

14. 39270 

8.97 

24 

14.72056 

9.46 

25 

15.02423 

9.96 

30 

16. 22632 

12.73 

35 

17.02826 

15.83 

40 

17. 57477 

19.14 

45 

17. 96167 

22.58 

50 

18. 24664 

26.  10 

55 

18. 46381 

29.68 

60 

18. 63420 

33.30 

65 

18. 77121 

36.94 

70 

18.88358 

40.60 

75 

18. 97751 

44.28 

Table  lv-4 
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t]  T?  T?  r T 

FALANCrr  SYSTE* 


N PIP  COM.P  ' PEP • DBM N 


1.00000 

2.00000 

3.00000 

4.00000 

5.00000 

0.06945  0.03030  0.01010  0.00631 
0.13793  0.06019  0.02006  0.01254 
0.20538  0.08962  0.02987  0.01867 
0.27172  0.11857  0.03952  0.02470 
0.33684  0.14699  0.04900  0.03062 

1.00000 

2.00000 

3.00000 

4.00000 
5.00000_ 

6.00000 

7.00000 

8.00000 
9.00000 

10.00000 

0 . 40067  0.17484  0.65028  0.03643 
0.46311  0.20209  0.06736  0.04210 
0.52405  0.22868  0.07623  0.04764 
0.58340  0.25457  0.08486  0.05304 
0.64104  0.27973  0.09324  0.05828 

6.00000 

7.00000 

8.00000 
9.00000 

10.00000 

1 1 .00000 
12.00000 

13.00000 

14.00000 

15.00000 

0.69688  0.30410  0.10136  0.0o336 
0.75082  0.32763  0.10921  0.06826 
0.80276  0.35030  0.11676  0.07298 
0.85261  0.37205  0.12401  0.07751 
0.90029  0.39286  0.13095  0.08185 

1 1 .00000 

12.00000 

13.00000 

14.00000 

15.00000 

16.00000 

17.00000 

18.00000 

19.00000 

20.00000 

6.94374  0.41269  0.13756  0.08593 
0.98890  0.43152  0.14384  0.03990 
1.02975  0.44935  0.14978  0.09362 
1.06826  0.46615  0.15538  0.09712 
1.10445  0.48194  0.16065  0.10041 

16.00000 

17.00000 

18.00000 

19.00000 

20.00000 

'21.00000 

22.00000 

23.00000 

24.00000 

25.00000 

1 .13835  0.49673  0.16537  0.10347 
1.16996  ' 0.51053  0.17017  0.10636 
1.19939  0.52337  0.17446  0.10904 
1.22671  0.53529  0.17843  0.11152 
1.25202  0.54634  0.13211  0.11382 

21.00000 

22.00000 

23.00000 

24.00000 

25.00000 
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117  ATT?':  I 


F-'iLA’iCZr  ' I 


N 

FiP 

COMP 

DBM 

PEP 

N 

1.00000 

0.04167 

0.04167 

0.04167 

0.04167 

1.00000 

2.00000 

O.OS276 

0.08276 

0.08276 

0.03276 

2.00000 

3.00000 

0.12523 

0.12323 

0.12324 

0.12323 

3.00000 

4.00000 

0 . 16303 

0.16303 

0.16304 

0.16303 

4.00000 

r,.  00000 

0.2021 1 

0.20211  ■ 

0.20212 

0.2071 1 

5.00000 

6.00000 

0.24040 

0.24040 

0.24042 

0.24041 

6.00000 

7.00000 

0.27787 

0.27787 

0.27788 

0.27787 

7.00000 

8.00000 

0.31443 

0.31443 

0.31445 

0.31443 

8,00000 

9.00000 

0.35004 

0.35004 

0.35006 

0.35004 

9,00000 

10.00000  1 

0.38463 

0.38463 

0.33465 

0.38463 

10.00000 

11.00000 

0.41813 

0.41813 

0.41815 

0.41813 

11,00000 

12.00000 

0.45049 

0.45049 

0.45052 

0.45050 

12.00000 

13.00000 

0.48165 

0.48165 

0.48168 

0.48166 

13.00000 

14.00000 

0.51156 

0.51156 

0.51160 

0.51157 

14.00000 

IG. 00000 

0.54017 

0.54017 

0.54021 

0.54018 

15.00000 

16.00000  ; 

0.56744 

0.56744 

0.56748 

0 • S6743 

16.00000 

17.00000  1 

0.59334 

0.59334 

0.59338 

0.59335 

17.00000 

18.C0C00  ! 

0.617S5 

0.61735 

0.61789 

0.61786 

18.00000 

19.00000  } 

0.64096 

0.64096 

0.64099 

■0 . 64096 

19.00000 

20.00000  ! 

0.66267 

0.66267 

0.66271 

0.66267 

20.00000 

21.00000 

0 . 68300 

0 .6C3uO 

0 . 6830-* 

0 1 o S 3 0 0 

21 .00000 

22.00000 

0.70197 

0.70197 

0.70201 

0.70198 

22.00000 

23.00000 

0.71963 

0.71963 

0.71967 

0.71964 

23.00000 

24.00000 

0.73603 

0.73602 

0.73607 

0.73603 

24.00000 

25.00000 

0.75121 

0.75121 

0.75125 

0.75121 

^Z^Ll-COO^O 

30 . OOOyO 

0.811^2 

0.31132 

0.S1136 

0.81132 

30 . 00000 

35.00000 

0.85141 

0.85141 

0.85146 

0.85142 

35.00000 

40.00000 

0 * S7G7  4 

0.37874, 

0.87879 

0.87875 

40.00000 

45.00000 

0.89809 

0.89309 

0.39814 

0.89809 

45.00000 

50.00000 

0.91233 

0.91233 

0.91238 

0.91234 

50.00000 

55.00000 

0.92319 

0.92319 

0.9232<» 

0.92320  1 

55.00000 

60.00000 

C. 93171 

0.93171 

0.93177 

0.93 172 

60.00000 

65.00000 

0.93856 

0.93B56 

0.93862 

0.93357 

65.00000 

70.00000 

0.94419 

0.94418 

0.94424 

0.94419 

70.00000 

75.00000 

0.94885 

0.94887 

0.94893 

0.94888 

75.00000 

pr-f' J A • ' ’ 

" ' ; 

'■>  ' • 

r 

L/  i f i ^ ‘ 

* 

'V  -f' 
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As  is  expected,  Tr{1)  = 3.0  which  is  in  accord  with 
equation  IV-43.  Purthernore,  Table  IV-6  does  verify  that  the 
selection  of  the  u{i)  for  all  n (i) , i ^ T,  does  in  fact  provide 
a totally  balanced  systea  as  the  utilisation  factor  is  the  sane 
for  each  processor  for  all  N.  Figure  IV-7  shows  the  asymptotic 
values  and  the  Tr(S)  from  the  model  for  various  value  of  m.  It 
is  obvious  that  as  m is  reduced  the  horivontal  asymptote  drops 
and  the  large  N linearization  asymptote,  equation  IV-36,  is 
reduced  in  slope.  Consequently,  a corresponding  improvement  is 
gained  for  all  N and  the  N*  is  increased  as  well.  However,  as 
shown  in  Figure  IV-8,  a corresponding  decrease  in  utilization  is 
encountered  as  it  should  be.  'What  is  mcst  interesting,  however, 
is  that  for  a corresponding  N*  and  m,  0.375  < m < 1.33,  there  is 
a difference  of  less  that  10  percent  in  the  balanced 
utilization.  Furthermore,  this  point  is  occurring  near  the  knee 
of  the  curve,  indicating  that  N*  may  be  thought  of  as  the 
optimum  load. 
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Chapter  V:  GFSS  SI’IULSTICN 


Discussion  and  Description 


It  is  recognized  that  in  the  develcpaent  of  any  computer 
system  model,  validation  against  the  actual  system  being 
modelled  is  highly  desirable.  However,  we  are  developing  a 
model  for  general  purpose  use  in  analysis  and  design  of  a 
computer  network  within  the  constraint  of  functionalization.  We 
may  not  have  the  actual  system  at  hand  for  validation. 
Consequently,  we  will  verify  the  analytic  model  developed  for 
internal  consistency  and  to  ensure  that  it  is  performing  as  we 
have  anticipated.  Then,  through  the  designer/architect's 
implementation  and  analysis  of  the  computer  network  functions,  a 
realistic  match  to  the  model  must  he  achieved.  With  this 
approach,  we  can  feel  confident  of  having  a good  model  to 
represent  the  desired  computer  network. 

The  procedure  for  verification  of  the  analytic  model  is  to 
implement  the  analytic  model  as  an  event  driven  simulation  model 
and  compare  the  results  based  on  cur  primary  performance 
criterion  of  mean  response  time.  The  simulation  system  selected 
is  the  General  Purpose  Simulation  System  (GPSS)  [ ISS  71],  This 
system  was  selected  because  of  its  well  established  use  for 
event  driven  simulation  of  computer  systems  ([Colella  et  al  74], 
[Schriher  74]),  its  ready  availability,  the  convenience  of  use 
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at  the  higher,  laacro  level  and  the  ability  to  readily  control 
desired  parameters. 

The  representation  and  implementation  in  GPSS  of  the 
networic  as  determined  in  Chapter  III  is  included  in  Appendix  E. 
In  this  GPSS  simulation,  exact  deterministic  routings  as 
described  for  r (k) , 1 < k < 4,  are  implemented,  provisions  for 
other  service  time  distributions  are  made,  and  the  capability 
for  a much  mote  detailed  repent  generation  is  available. 
Although  exponential  service  distributions  are  used  for  exact 
comparison  to  the  EQN  model,  other  distributions  can  be  used  in 
GPSS.  Additionally,  ve  will  be  able  to  readily  note  the  effects 
of  our  assumptions  used  in  assimilating  the  k routes  into  one. 
A sample  of  a typical  GPSS  output  is  at  Figure  V-1.  It  is  easy 
to  see  in  this  figure  some  of  the  additional  information 
provided  by  the  simulation  compared  to  that  provided  by  the 
Basic  Model  calculations  of  Chapter  IV.  First,  general 
information  concerning  the  simulation  run  is  tabulated.  At  this 
point,  it  should  be  noted  that  an  initialixation  of  the 
simulation  run  is  required  prior  to  data  collection  [ Schriber 
74].  This  initialization  is  necessary  to  force  the  simulation 
into  a steady  state  condition.  Other  data  is  the  mean  think 
times  zk  for  each  route  with  a k correspondence  of  {|1,2,3,U|} 
::=  {|A,B,C,C|].  The  mean  service  times  for  each  function  and 
the  type  of  distribution,  in  this  case  exponential,  are  listed. 
This  is  followed  by  information  concerning  both  the  queue  and 
processor  statistics.  Next  a matrix  showing  the  absolute  count 
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on  row  one  and  the  relative  frequency  {normalized  to  1000)  on 
row  two  of  the  k job  types  (columns)  entering  the  CFN  is 
displayed.  This  relative  frequency  is  the  job  lead  distribution 
in  the  CPN.  The  mean  composite  response  time  for  all  job  types 
and  each  job  type's  mean  response  time  (PESRk)  and  their 
respective  standard  deviations  are  enumerated  (in  milliseconds) . 
This  data  is  gathered  from  tables  maintained  by  GPSS  which  may 
also  be  individually  listed,  if  desired,  as  well  as  a 
probability  density  function  plot  of  the  response  times. 
Lastly,  a matrix,  processor  entry  count  by  routes,  showing  the 
number  of  visits  of  each  job  type  (columns  1-U)  to  each  node 
(rows  1-5)  , is  tabulated.  Note  that  these  values  compare  to  the 
vk(i)  with  an  i correspondence  of  { 1 1 ,2,3, 4,5  |}  ::  = 
{ I T, F, C , D, P I } . The  fifth  column  is  the  sum  of  all  jobs  entering 
a functional  processor. 
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GFSS  Sirini  A"‘ICN  OF 
CASIC  ^'CTEL,  N = T- 


RonsiR  OF  joes  is  io 

IRITIAUZAriCR  HOODS  ADD  6.00 
CATA  COLLRCTICD  HOODS  ARE  2.00 
THINK  TINES  ARE  • A:  ISOOO 
B:  1S000 

C;  15000 
D:  15000 

SEETICE  TIRES  (NSEC)  ARE  - 


TBE  JOB  COORT  BT  ROOTS 
FOLIHORD  HATRIX  RTCRT 


BCR/COIDHB 

1 

2 


2T1B 

729 


rxP: 

EXP 

600 

coflp: 

EIP 

1375 

DBn: 

EXP 

4125 

PEP: 

EXP 

6600 

559 

150 


r-T 


r-r 


U' 


iTi  i ‘ 


15  e 


i i ■ 


7i 


198 

53 


251 

67 


COlfPOSITI 

BEAR  R TIRE 

IS  4317.410 

TAELE  BESB1 

ERTRIES 

IE  TABLE 

REAR  AR60RERT 

2718 

891.017 

TABLE  RBSD2 

EBTEIfS 

IE  TABLE 

REAR  ARGDRERT 

559 

10747.160 

TABLE  RESR3 

ENTBIES 

IN  TABLE 

HEAR  ARGOBENT 

198 

21322.492 

TABLE  BBSR4 

EWTRIBS 

IE  TABLE 

REAR  ARGORERT 

251 

13686.769 

QOEOE  STATISTICS 

QO  EOE 

AT  ED  AGE 

AT  ERA6E 

CONTENTS 

TIHE/TBARS 

TEBflQ 

7.765 

14970.199 

PIPQ 

.582 

885.714 

COBPQ 

.603 

1962.601 

DBflQ 

.543 

5166.402 

PEPO 

.504 

8078.222 

PROCESSOE 

STATISTICS 

PXCILITT 

AVERAGE 

NORBER  ATERAGB 

0TIII2ATI0R 

ENTRIES  TIRE/TBAR 

ripp 

.394 

4734.  599.330 

CCBP? 

.437 

2215  1421.429 

cBnp 

.425 

757  4048.403 

PEPP 

.412 

450  6595.597 

PROCESSOR 

ENT  CODNT  BE 

ROOTS 

FOLLRORC 

BATRIX  ERT 

STARDABD  OETIATIOB 

899.000 


STARDABD  DETIATIOB 

5360.000 


STANDARD  DETIATIOB 
9RR0.000 


STARDABD  DETIATIOB 

8544.000 


RCR/COIOBB 


1 

2 

3 

4 

5 

2718 

559 

198 

251 

3726 

2718 

1118 

396 

502 

4734 

0 

1118 

594 

503 

2215 

0 

559 

198 

0 

559 

0 

0 

19B 

251 

449 

•iqtltp  V-1 
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Conpariscn  to  Basic  Jlodel  Results 


A ccaparison  of  the  Basic  Kodel  (BB)  results  with  those 
froffl  the  GPSS  simulation  (S)  is  at  Table  V-1  which  follows.  Due 
to  the  relatively  large  amount  cf  computer  time  required  for  the 
GPSS  simulation,  only  major  job  points  have  been  simulated.  The 
following  initialization  and  run  times  for  the  GPSS  simulation 
were  used: 


N 

Initialization 

(hours) 

Tine 

Pun  Time 
(hours) 

10 

6 

2 

20 

12 

4 

30 

18 

6 

♦ UO 

6 

2 

* Due  to  the  lengthy  computations  involved,  this  sample 
is  run  for  an  indication  only. 

^ 10, 

20,  and  30,  a copy 

of  the 

output  is 

at  Figures  7-1, 

id  7-3 

respectively. 

Comparison  of 
Basic  nodel  and 

GPSS 

Q(T) 

Tr  (H) 

^lerror 

N 

GPSS  EB 

GPSS 

EH 

10 

7.77  7.69 

4.32 

4.50 

4.2 

20 

14.19  13.25 

6.13 

7.64 

24.6 

30 

18.73  16.23 

9.C2 

12.73 

41.1 

*tt0 

20.66  17.57 

14.04 

19.14 

36.3 

Table  7-1 
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GFSS  ST’^nLS.TTCN  O? 
"ASTC  r.CD3L,  N = 2C 


lonazf  cf  jots  is  to 

ISlTUtLIZATICN  HOURS  ARE  12.00 
DATA  CCIISCTIOR  HOOBS  ARE  «.00 
THINK  IIHES  ARE  - A: 


1SOOO 
B:  15000 

C;  15000 
D:  15000 

SERVICE  TIRES  (HSEC)  ABE  - 


n-rri 


/ 


A f A A ^ 


TBE  JCB  COONT  BT  ROOTE 
EtJLLBORC  HATPIX  STCHT 


BCH/CCIOHR 


10055 

767 


CCBPOSITI  HEAR  R TIRE  IS 
TABLE  RESR1 
ENTRIES  IS  TABLE 
10*55 

TABLE  RESR2 
ENTRIES  IN  TABLE 
1759 

TABLE  BESR3 
ENTRIES  IN  TABLE 
601 

TABLE  RESR* 

ENTRIES  IN  TABLE 
815 


TIP: 

EIP  600 

^ i <N  At'  • il  < 

coap: 

EXP  1375 

Daa: 

EXP  *125 

PEPS 

EXP  6600 

2 

3 

1759 

601 

815 

129 

46 

59 

6133. 

80* 

aEAE 

ARGORENT 

STAIDARO  DEVIATXOI 

1533.592 

1 *07.000 

HSAtf 

ARGORENT 

SIAEDARD  OEflATIOl 

17737.531 

8096.000 

REAM 

ARGORENT 

r 

OBVIATION 

32938.  9*5 

12832.000 

REAR 

ARGORENT 

STA 

"I  DEVIATION 

20335.535 

10*16.000 

^ rryj 

- I 


j 

i! 


1 


QOEOB  STATIS71CS 


QO  BUS 

AVERAGE 

average 

CONTENTS 

TISE/TRANS 

TEBBC 

1*.19* 

1*978.070 

PIPQ 

1.906 

15*7.6*0 

coapc 

1.**2 

2986.607 

DB8Q 

1.*22 

8676.1*0 

PEPQ 

1.13* 

11533.371 

PI0CESS08 

STATISTICS 

FAC  II.  ITT 

AVERAGE 

ROBBER 

ATFRAGB 

OTIITZATIOR 

ENTRIES 

TIfle/TSA8 

riPF 

.688 

16808  • 

589.957 

CCHPf 

.675 

6*53 

1399.170 

DEBT 

.693 

2359 

9230.937 

PIPE 

.6*9 

1*16 

6601.636 

PROCESSOR 

ENT  COONT  BI 

ROOTE 

POLLRORO 

HATRIX  ENT 

BOV/COIOBI 

1 

2 

3 

* 

5 

1 

10*55 

1759 

601 

815 

13630 

2 

10*56 

3519 

1203 

1630 

168  08 

3 

0 

3518 

1802 

1630 

6950 

9 

0 

1760 

601 

0 

1760 

5 

0 

0 

600 

816 

1*16 

1 

Plan  r = 

V-? 

I 


1 


F 

t 

i 
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ST !"nL '■v  0? 


BAS^c  ICRTI, 


KDRBER  OF  JOES  IS  10 
IltlTlALIIXTICII  HOODS  ABE  IS. 00 
DATA  COLIECTIOD  HOODS  ABE  6.00 
THINK  TIBBS  ABE  - A:  15000 
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Proa  the  results  in  Table  V-1,  it  is  apparent  that  the 
Basic  nodel  may  not  predict  response  tiaes  as  accurately  as  one 
would  desire.  This  is  true  even  when  one  considers  that  for  EQN 
Bodelling,  we  should  be  satisfied  with  response  tiae  accuracies 
with  a tolerance  of  up  to  30*  [ Bu2en  and  Denning  77].  Saall 
values  of  N are  producing  gocd  results  and  tend  toward  the 
unacceptable  range  near  B = 20.  The  figures  do,  however, 
provide  insights  as  anticipated,  Pecall  that  the  design 
criteria  for  the  balanced  systea  was  with  a saturation  of  N*  = 
24.  In  other  words,  as  long  as  we  are  below  saturation,  we 
obtain  usable  results  but,  above  saturation,  the  assuaptions 
aade  in  the  developaent  of  the  Basic  Model  show  their  effect 
draaatically . The  principal  assuapticn  aade  was  in  the 
assiailation  of  the  k job  routes  into  cne  in  order  to  utili2e  a 
quick  and  rapid  coaputation  facility  such  as  the  single  job  type 
ASQ.  As  explained  in  Chapter  IV,  we  have  assuaed  a job  load 
distribution  in  the  CPS  as  (.6,  .2,  .1,  ,1).  In  fact,  this  is 
not  true  as  can  be  verified  by  the  frequency  distribution  in  the 
job  count  by  route  aatrix  of  the  GPSS  output  which  shows  the 
following  lead  distributions: 


N 

Job  load  Distributiens 

19 

(.729,  .150, 

.053,  .067) 

20 

(.767,  .129, 

.044,  .059) 

30 

(.793,  .113, 

.038,  .054) 

Due  to  their  faster  response,  we  expect  to  find  a large  nuaber 
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of  route  1 jobs  and  a fewer  noaher  of  jets  of  routes  2,  3,  and  « 
in  the  CPN.  Since  route  1 jots  require  only  the  FIP  processing 
function,  whereas  the  others  require  two  units  of  processing 
froa  the  FIP  and  at  least  two  other  processing  functions,  it  is 
natural  to  find  this  skewed  job  load  distribution  froa  that 
assnaed.  Furthernore,  to  a large  degree,  this  frequency 
distribution  will  depend  on  the  network  structure  and  the 
processing  functions. 

The  utilization  of  the  processors  reaains  relatively 
balanced  as  intended  throughout  the  range  of  N as  shown  in  Table 

v-2. 

GPSS  Otilizaticn  Suaiary 
of  a Balanced  Systea 


R BN 

(all  i 1) 


GPSS 


PIP 

CCHP 

CBn 

PEP 

10 

.385 

.394 

.437 

.425 

.412 

20 

.663 

.688 

.675 

.693 

.64  9 

30 

.811 

.898 

.781 

.787 

.726 

ttO 

.879 

.983 

.813 

.813 

.804 

Table  V-2 

This  is  expected  since  it  is  realized  that  utilization  figures 
are  affected  less  by  the  job  load  distribution  than  ace  response 
tiaes  [ Buzen  and  Denning  77],  Again  the  effect  of  the  loading 
of  the  FIP  by  route  1 jobs  is  shown  by  a utilization  slightly 
greater  than  the  balanced  values  of  the  Basic  Hodel. 
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Suaaary  and  dotivation  for  Farther  aefinement 

The  Basic  Hodel  (EH)  is  shown  to  provide  good  results  when 
used  with  a saall  number  of  jobs,  that  is,  within  the  saturation 
design  criteria.  However,  near  the  upper  liait  of  the 
saturation  the  response  time  error  is  greater  than  desired.  As 
is  expected,  the  utilixations  provide  ouch  better  results 
throughout  the  range.  A significant  advantage  of  this  aodel  is 
its  sioplicity  of  use  when  a general  puipose  progran  such  as  ASC 
is  available.  Both  trade*off  analysis  and  balanced  design  aay 
be  accoaplished  very  quickly.  Conseguen tly,  it  is  expected  that 
the  Sasic  Hodel  aay  be  effectively  utilixed  to  give  insights  but 
we  are  motivated  to  search  for  yet  a more  accurate  model, 
realizing  that  a trade-off  in  terms  cf  simplicity  both  in  use 
and  in  computation  may  be  required.  He  should  also  note  that  in 
the  Basic  Hodel  the  service  rates  and  service  disciplines  at  the 
functionalized  processors  mast  be  identical  for  all  routes. 
This  restriction  is  not  felt  to  be  significant  since  a 
particular  route  may  require  multiple  visits  at  a functionalized 
processor.  The  multiplicity  of  visits  may  account  for  differing 
processing  requirements  by  route. 
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Chapter  VI:  POOR  JOB  TYPES  MODEL  USING  INTEGER  ITERATION 


General  Approach 


The  greatest  disadvantage  of  the  Basic  Model  presented  in 
Chapter  IV  is  the  assnaption  required  to  assinilate  the  k job 
types  (routes)  into  one.  This  assimilation  was  done  to  simplify 
the  computation  procedure  and  to  take  advantage  of  the  Duke  ASQ 
general  purpose  program.  Another  approach  would  be  to  use  one 
of  the  general  purpose  programs  with  multi-job  capabilities  such 
as  the  ASQ  from  Information  Research  Associates  [1975].  This 
possibility  was  investigated  and  some  trial  runs  were  made  using 
a version  available  at  Federal  Simulations  in  the  Department  of 
Defense  executing  on  a Honeywell  (GE)  635.  However,  due  to  the 
constraints  of  providing  only  three  job  types  and  a total  of  20 
jobs,  this  program  was  not  felt  to  be  entirely  satisfactory  for 
our  purposes.  Furthermore,  the  amount  of  CPU  execute  time 
required  was  relatively  high,  making  the  program  expensive  to 
use.  Since  the  exact  algorithms  used  for  solution  were  not 
immediately  available,  it  is  assumed  that  additional  time  was 
required  because  of  the  general  purpose  nature  of  the  program. 
The  other  simulation  program,  QNET4,  has  not  been  released  for 
use  outside  IBM.  It  too,  however,  restricts  the  user  to  three 
job  types  instead  of  our  desired  four. 
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Consequently,  another  approach,  which  is  pursued  in  this 
chapter,  is  to  develop  a special  purpose  progran  enploying  the 
solution  foras  of  Bu2en  [1973]  and  flilliaas  and  Bhandivad 
[1976].  Not  only  are  we  able  tc  tailor  the  progcaa  to  model  our 
computer  network,  but  we  can  ensure  that  some  of  the  most 
computationally  efficient  techniques  and  procedures  are 
employed. 
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Solution  Forms  with  Normalization  Constant 

The  general  solution  form  for  a single  job  type  in  a closed 
FQN  model  ([Gordon  and  Newell  67],  [Williams  and  Bhandiwad  76]) 
is  based  on  the  state  of  the  network  at  any  time,  denoted  by  a 
job  vector,  n = [ N ( 1)  , N (2)  , . . . , N (n)  ] , where  N(i)  is  the  number 
of  jobs  at  node  i (both  those  in  service  and  those  waiting  in 
the  queue).  Since  each  job  is  either  in  service  or  waiting  in 
cne  of  the  queues,  we  must  have 

2 N(i)  = N (VI-1) 

i 

The  solution  form  for  the  probability  of  being  in  a particular 
state,  often  referred  to  as  the  product  form  solution,  is 

Pr[N(1)  ,N{2)  ,...,N  (|H|)  ] = 

(1  / G{N))  n X(i)  **[N(i)]  (VI-2) 

i 

where 


X(i)  v(i)  / u{i)  (VI-3) 

Note  that  the  v(i)  represents  the  relative  visits  to  node  i as 
solved  by  the  balance  equations  (IV-ttO) . Consequently  the  X*s 
are  often  referred  to  as  the  relative  utilizations.  The 
normalizing  constant,  G,  is  determined  by  the  requirement  that 
the  probabilities  summed  over  all  states  is  unity.  Therefore 

E n [X  (i)  **  (N  (D)  ] 

St{N,|f1|)  i 


G(N) 


(VI-U) 
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where  the  sua  is  over  all  n in  the  state  space  of  N jobs  and  | fl  | 
nodes,  denoted  as  St(N,|Hl).  It  is  then  shown  that  the 
utilization  of  the  i-th  node  with  N jobs  in  the  network  is 

0(i,N)  = Pr[N{i)  >1]  = X(i)  ♦ G(N-1)  / G(N)  (VI-5) 

From  the  utilization,  the  thruput  say  be  coaputed  as  follows. 

T (i,N)  = 0 (i,N)  • u(i)  (VI-6a) 

= [X(i)  / 2(i)  ] * G(n-1)  / G(N)  {VI-7) 

where  z (i)  is  the  aean  service  time  at  node  i (equal  by 
definition  to  1 / u (i)  ) . 
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Solution  Forms  for  Bultiple  Job  Types 

At  this  juncture  we  should  recall  that  the  solution  forms 
in  the  previous  paragraph  are  for  the  situation  of  only  one  jot 
type.  Basket  et  al  [1975]  have  extended  the  work  of  Jackson 
[1957]  and  Gordon  and  Newell  [1967]  to  orovide  solution  forms 
for  multiple  job  types.  Euten  [ 1971,  197!!]  developed  the  idea 
of  using  the  normalizing  constant  to  characterize  network 
properties.  He  also  provided  the  recursive  specification  of  the 
normalizing  constant  as  a ccmputational  technigue. 
Subsequently,  Williams  and  Ehandiwad  [1976]  neatly  present  the 
solution  forms  for  two  job  types  requiring  an  extension  to  k jot 

types  in  our  case.  We  can  designate  Nk(n)  as  the  number  of  jobs 

of  the  k-th  route  or  job  type  in  the  complete  network  n such 
that  the  sum  over  k of  Nk(n)  equals  N,  the  total  number  of  jobs 

in  the  network.  Likewise,  Nk(i)  is  the  number  of  jobs  of  the 

k-th  route  at  the  i-th  node  in  the  network.  Then  for  two  jot 
types  or  routes,  k = 1,  2,  the  solution  form  is 

Pr[Nl(1),N2(1)  ;...;N1(i)  ,N2(i);...;IJ1{|f1|),N2{|H))  ] 

= [1/(52]  n FECD2  (i)  {VI-8a) 

i 

where 


PRCD2(i)  ;;=  ([  Hi  (i)  ♦ N2(i)]!  / [N1(i)!  ?'2  (i)  ! ]) 

* X1(i)  [Nl(i)]  * X2(i)  ••  [N2{i)  ] 

(VI-8b) 


where  G2  is  again  the  normalizing  constant  but  in  this  case  for 
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two  job  types.  Thus 

02-1.  n PR0C2(i>  (VI-9) 

St[Nl  (n),N2(n),|H|  ] j 

for  all  [ N1  (i)  ,N2  (i)  ] such  that 

Nlc(i)  = N)c{n)  for  h = 1,  2. 
i 

Extending  equation  VI-8  to  k job  types  or  specifically  in 
our  case  where  k=4  we  have 


Pr[  N1  (1)  . ,Nk  (1)  ; . (i)  ,Nk  (i)  ] 

= [ 1 / Gk ] n FRCDk  (i) 


(Vl-lOa) 


where 


C ^ Nk(i)  ]! 

Prodk(i)  = ♦ ** 

^ (Vl-lOb) 

K 


L 

k 

n [iik(i) !] 

K 

and  again  G (k)  is  the  notnalizf'ng  constant  for  k job  types: 

«=  E n PRODk(i)  (VI-11) 

St[N1,N2,...,Nk,| n|  ] ' 


For  the  sake  of  sinplicity  and  clarity  in  the  prograo 
implementation  and  in  further  notation,  the  k superscript, 
indicative  of  the  job  type  or  route,  will  he  specifically 
designated  by  the  job  types  as  A,  E,  C,  and  D.  Likewise  the 
corresponding  relative  utilizations  for  each  job  type  at  the  i- 
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th  node  can  fce  respectively  designatea  as  M (i) , X (i)  , I (i)  and 
Z (i)  where 


XI  (i) 

fcecoaes 

V (i) 

X2  (i) 

becoaes 

X |i) 

X3  (i) 

hecoaes 

Y (i) 

X4  (i) 

hecoaes 

2 (i) 

Othei 

: variables  will 

have 

the 

lieu 

of  the 

superscri 

.pt  appen 

v1  (i) 

becomes 

vA(i) 

z2  (i) 

becomes 

2E  (i) 

Ni 

becomes 

NC 

T4  (i) 

becomes 

TD  (i) 

and  so 

forth . 

Recall  that 

the  vk(i) 

is  the 

type 

(route)  k fo« 

ind  by 

th 

zk(i) 

= 

1 / u(i 

-)  . 

In 

specification  of  Chapter  III,  the  s (i)  = <u  (i)  , a (i) > are 
identical  fcr  all  routes.  Strictly  speaking,  we  have  the 
capability  to  specify  s(i)  for  each  route.  Following  the 
previous  nouencla ture,  we  would  have  sk  (i)  = <uk(i),  ak  (i)  >. 

Wow  for  our  specific  case  the  solution  fora  becones 

Pr[  NA  ( 1)  ,riB  ( 1)  ,NC  (1)  ,ND  (1)  ;.  . . ;NA  (4)  ,NB  (4)  , KC(4)  ,ND  (4)  ] 

= [ 1 / C (NA,Pe,NC,ND)  ] n Prod4  (i)  (VI-12a) 
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where 


Prc<ia(i)  = [ NA  (i)  ♦ NB  (i)  ♦ NC(i)  ♦ ND(i)]!  (Vl-I2b) 

/ [ NA  (i)  ! NE  (i)  ! NC  (i)  ! ND  (i)  ! ] 

♦ « (i)  ♦*  [NA(i)  ] ♦ X (i)  •*  [ NB(i)  ] 

* Y (i)  [ NC  (i)  ] * 2 (i)  *♦  [ HD<i)  ] 


Looking  at  four  job  types,  the  utilization  at  node  i is 


0 (i,NA,NB,NC,ND)  = (71-13) 

[B(i)  G (NA-1,NB,KC,ND)  ♦ 7(i)  G (N  A , NE- 1 , NC,  N D) 

♦ T (i)  G (NA,NB,NC-1  ,ND)  ♦ Z (i)  G (NA,  NB,NC,  ND- 1)  ] 

/ G {NA,VB,NC, NC) 


Similarly  to  equation  Vl-6,  the  thruput  for  each  job  type  is 

TA  (i)  = a (i,NA,NE,NC,SD)  ♦ u (i)  (VI-14a) 

= [X(i)  / 2A  (i)  ] G {NA-1,NE,NC,ND)  / G (NA  , NB  ,NC  , ND) 

(VI-14b) 

This  equation  holds  similarly  for  job  types  B,  C,  and  D. 
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Ccaputation  of  the  Nor aaliza tion  Constant 


From  the 

previous 

sections 

we 

have 

seen  that 

the 

normalization 

constant 

, G,  is 

central 

to 

solutions 

for 

utilization  and 

thruput. 

Therefore, 

we 

need 

to 

examine 

the 

computational 

aethods 

available 

for 

G. 

One 

of  the 

more 

efficient  methods  in  the 

single  job 

type  case 

is 

presented 

by 

Euzen  [1973]  and  then  by  Nilliatos  and  Ehandiwad  [ 1976  ], 

g(nA,  a)  = g(nA,  a-l)  ♦ H (a)  g(nA-1,  a)  {VI-15) 

The  latter  reference  extends  the  single  job  type  algoritha  to 
aultiple  job  classes.  In  the  case  of  tac  job  types,  it  is 

g(nA,  nB,  a)  = g(nA,  nE,  a-l)  (VI-16) 

♦ W (a)  g (nA- 1 , nB , a) 

♦ A (a)  g (nA,  nB- 1 , a) 

Extended  to  font  job  types,  it  is 


g(nA,  nB,  nC, 

nD, 

a)  = 

g (nA, 

n 

B , nC 

, nD, 

B- 

1) 

(VI- 

17) 

♦ i(B) 

g 

(nA-l 

, nl 

E, 

nC 

, nD, 

B) 

♦ A{B) 

g 

(nA, 

nB- 

Ir 

nC 

, nD, 

B) 

♦ Y(*) 

g 

(nA, 

nB, 

nC 

-1 

, nD, 

a) 

♦ Z(B) 

g 

(nA, 

nB, 

nC 

nD- 1 , 

■) 

The  a-th  subscript 

indicates 

the  B- 

th 

node 

in 

th 

e 

network  and 

in 

our  case  has  a aax 

iauB 

value 

of  mi 

• 

The 

nk  is 

such  that  1 < 

nk 

< NA.  The  final  value  of  g,  G,  is  obtained  when  a reaches  its 
aazinaa  value. 

It  is  seen  froa  the  above  equations  that  this  is  a 
recursive  specification  such  that  the  current  value  of  g is 


Chapter  VI 


based  on  previously  ccaputed  values.  It  is  a matter  of 
establishing  the  initial  values  so  that  the  basic  equation  of 
the  solution  fora  is  aet.  For  the  single  job  type,  g(0,a)  = 1 V 
a and  g(n&,0)  =0  V nA  are  the  initialization  requirements. 
Furtheraore,  it  is  noted  that  although  g is  in  the  fora  of  a 
two-dimensional  aatria,  only  the  current  column  vector  need  be 
retained  for  computation  and  not  the  entire  matrix. 
Consequently,  computational  storage  reguireaents  are  reduced 
considerably. 


In  analyzing  the  equations  VI-16  and  VI-17  as  extensions  of 
equation  VI-15,  we  note  that  the  ainimui  dimensions  required  for 
computation  of  G are  NA,  NE,  NC,  and  NE.  The  m value  for  the 
nodes  is  maintained  only  as  the  current  computation  of  g.  In 
order  to  visualize  more  clearly  how  the  computational  procedure 
is  extended  to  four  or  more  job  types.  Figure  VI-1  illustrates., 
the  case  of  two  job  types.  Here  we  have  let  a and  b be  the 
running  variables  for  NA  and  NB  respectively  in  lieu  of  the  nA 
and  nB.  In  this  case  the  full  g would  be  a three-dimensional 
matrix  or  cube  and  the  initial  values  vould  lie  adjacent  to  the 
two  faces  of  the  cube  as  shown.  These  initial  values  represent 
the  case  where  one  of  the  job  types  has  a zero  value,  e.g. 
g(0,b,a)  or  g(b,m)  and  g(a,0,a)  or  g(a,n).  But  to  obtain  the 
g(a,m)  and  g(b,m)  we  need  to  employ  the  procedure  for  the  single 
job  type  case. 
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Schema  for  G CoKfutaticn 
for  Twc  Goh  ypes 


M 


Z 


M 


z: 


z 


g(a,m) 

^ G(NA,1m|) 


1 


FiGiir® 
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As  explained  in  the  single  job  type  case,  it  is  not  necessary  to 
compute  the  entice  t vo-dinensicnal  iiatcix  represented  in  Figure 
VI- 1 as  j '=ace  to  the  cute  but  one  colunn  vector  for  each  n 
value  vould  suffice.  likewise  only  one  "a  slice"  of  the  cube  is 
required  for  the  current  coaputation  of  g.  In  order  to  siaplify 
the  coaputational  algorithm,  however,  we  will  compute  the  entire 
initial  value  aatcices  although  the  a diaension  of  the  final  g 
coaputation  will  be  eliminated. 

Using  the  procedure  just  described,  a G Computation  Tree 
(GCT)  is  developed: 


Computation  Tree  for  Two  Job  Types 


g 


g(a,a) 


(a,b,n) 

\ 

g 


(b,m) 


Figure  VI-2 


For 

more 

than  two 

jobs,  hypercubes  and  hypecplanes 

result.  A 

GCT 

is 

shown  at 

Figure  Vl-3  for  three  jobs  and  at 

Figure  VI-4 

for 

four 

jobs. 

our  case  of  interest.  The 

following 

correspondence  of 

running  variables  is  used  throughout: 

n A 

becomes 

a 

nB 

becoaes 

b 

nc 

becomes 

c 

nD 

becoaes 

d 
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This  conpatational  procedure  is  inplemented  by  the  progran 
listed  at  Appendix  C. 
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G Conputaticn  Tree 
for  Three  Job  Types 


g{a,b,c,iD) 


Figure  Vl-3 


1 

I 

I 
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G Conputaticn  Tree 
for  Four  Job  Types 


Figure  71-4 
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In  order  to  provide  additional  assurance  and  a greater 
degree  of  verification,  several  test  cases  of  the  G 
Noraalizaticn  Coaputation  were  considered. 

1.  Single  job  type  coapatations  were  ccapared  against 

Buzen's  example  [1973]  and  against  the  Duke  ASC 
employing  the  five  network  nodes  as  described  in 
Chapter  III. 

2.  Two  job  type  conputations  were  compared  against  a small 
two-node  network  computed  by  hand  and  also  run  using 
the  Information  Research  Associates  [1975]  ASQ  at 
Federal  Simulations. 

3.  Three  job  type  conputations  were  run  at  Federal 

Simulations  within  the  constraints  of  the  ASC 
implemented  there.  The  model  run  was  similar  to  that 
of  Chapter  III,  Only  r («)  was  eliminated  and  g(3)  was 
increased  from  C.1  to  0.2. 

All  of  the  above  cases  yielded  exact  comparisons  for  the 
utilization  and  thruput.  It  should  be  noted  that  a direct 
comparison  of  G values  was  not  made  since  they  are  a function  of 
the  relative  utilizations  which  in  turn  are  dependent  on  the 
solution  to  the  balance  eguations  (IV-OO)  which  dc  not  have  a 
unique  solution. 
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Model  Fornulation  and  Implenentaticn 

By  using  equations  VI-13  and  Vl-lh  and  the  computation  of 
the  normalization  constant,  we  can  compute  the  utilization  and 
the  thruput  at  a particular  node.  Hhat  we  are  interested  in 
obtaining,  however,  is  the  response  time  as  a function  of  the 
number  of  users  Tr  (N) . Hecall  from  Figure  Vl-5  (same  as  Figure 
Iv-1)  that  the  response  time  is  the  time  a job  spends  in  the 
CPN. 


Figure  Vl-5 


It  is  recognized  that  we  can  apply  Norton's  equivalent 
[Chandy,  Herzog  and  Woo,  75a  and  75b]  to  the  Computer  Processing 
Network.  The  CPN  is  then  replaced  by  its  short-circuited 
thruput  equivalent.  The  short  circuit  thruput  for  each  job  type 
is  computed  directly  from  equation  VI- 1«  as  long  as  the  vk(i) 
Vi  > 1 V k are  computed  based  on  vk(1)  = 1 v k.  Now  we  can 
apply  a job  load  vector  n = (NA,NB,HC,ND)  to  the  CPN  and 


W 


I 


I 

I 

I 
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determine  the  equivalent  short-circuit  thruput.  In  other  words, 
with  the  CPN  loaded  by  n,  the  equivalent  thruput  of  the  CPN  for 
each  job  route,  Tk  (CPN) , is  obtained.  In  fact,  we  may  compute 
Tk  (CPN)  as  a function  of  the  number  of  each  job  type  in  the  CPN 
(Nk  (CPN) ) varied  over  its  range  of  values  to  fora  distinct  job 
vectors  for  Tk(CPN).  As  a result,  there  is  a Tk(CPN)  for  each 
value  of  Nk  (CPN)  . 

We  now  further  recognize  that  we  can  obtain  a set  of 
thruputs  by  varying  each  job  type  over  the  range  of  its  values 
in  the  CPN  while  holding  the  number  of  ether  job  types  fixed  at 
the  maximum  values  of  the  job  load  vector.  This  will  yield  a 
minimum  thruput  for  each  value  cf  each  job  type  since  the  CPN  is 
maximally  loaded  with  the  remaining  job  types.  As  a 
consequence,  we  will  obtain  an  upper  bcund  on  response  tine.  To 
obtain  this  response  time  we  have  used  the  Duke  ASQ  for  the 
solutions  as  a matter  of  convenience.  We  represent  the 
equivalent  network  as  two  nodes,  the  TIEH  and  the  CPN,  foi:  each 
of  the  k job  types. 


Figure  VI-6 


I 


The  TENN 
equation 


is  as  represented 
ITI-2  as  u (T)  = 


previously  in  Chapter  III  and  by 
N / 2 where  N is  now  Nk  and  Z would 
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he  zk.  The  effective  service  rates  of  the  CPN  is  the  set  of 
thruputs  obtained  for  each  job  type  such  that 


uk  (CPN)  = {Tk(nk)}  1 < nk  < Nk  (VI-18) 

(The  Tk  is  for  the  CPN  node  of  Figure  VI-6  so  the  CPN  subscript 
is  implied  for  the  remainder  of  this  chapter.)  From  ASQ  we 
solve  for  the  number  of  jobs  of  the  respective  type  in  the  CPN 
(the  queue  length  of  the  CPN  node)  and  the  response  time.  The 
number  of  jobs  in  the  CPN  is  a new  approximation  to  the  job  load 
for  the  respective  job  type.  Doing  this  for  the  k job  types,  we 
now  have  a new  job  load  vector  n*.  With  this  new  job  load 
vector  we  can  again  compute  the  thruputs  and  obtain  a new 
loading  of  the  CPN.  We  are  restricted  to  integer  values, 
however,  since  the  computational  procedure  for  G is  restricted 
to  integer  job  values.  We  will  converge  to  an  approximate  value 
as  we  obtain  better  estimates  of  the  integer  values  for  the  jot 
load  vector. 


We  are  aware 

that 

there  will 

be  an  exact 

so lution 

for  the 

job  load  vector 

if 

we  were  to 

acconodate 

other  than 

integer 

values  [Baskett  et  al  75].  Starting  with  an  assumption  of 
maximum  load,  we  would  obtain  successively  better  approximations 
until  the  exact  value  is  reached.  In  our  case  we  will  obtain 
the  best  integer  load  vector  near  the  exact  values. 

The  algorithm  for  the  Iterative  Hodel  (IM)  is  formulated  as 
follows: 


r 


I 


1 
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1.  Determine  the  relative  atili2ations  Xk(i)  based  on 

vk  (T)  = 1 V k. 

2.  Assume  a maximum  job  load  vector  n = (N 1,N2, . . . ,Nk)  . 

3.  Solve  for  the  CPU  thruput  vector  for  each  job  type 
[T!:(nk)  ],  assuming  that  the  other  job  types  maximally 
load  the  CPN.  (The  program  at  Appendix  III  is  used.) 

tt.  Determine  a new  job  load  value  by  using  ASQ  as  in 
Figure  Vl-6  and  the  uk(CFN)  as  in  equation  VI-18. 

5.  Approximate  each  of  the  job  load  values  by  the  nearest 
positive  integer. 

6.  Compare  previous  job  lead  values. 

7.  If  there  is  no  change,  then  the  test  integer  load  value 
has  been  determined;  else  set  a new  load  vector  and 
repeat  starting  at  step  3. 

8.  Determine  the  [Trk(Nk)  ] from  the  TERH  queue  length  and 
equation  IV-1C. 

A composite  response  time  must  be  formulated  based  on  the  k 
response  times.  This  is  done  ty  weighting  the  response  times  by 
the  job  load  distribution  or  the  proportion  of  the  k job  types 
flowing  through  the  CPN  which  is  based  on  the  ratio  of  job  type 
thruputs. 

CTr  (N1  ,N2,. . . ,Nk)  = Tr(M)  = p (k)  Trk{Nk)  (71-19) 

k 

where 


1 

i 
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Tk  (Nk) 

p(K)  = 7 y (VI-20) 

2 Tk(Nk) 
k 

Once  the  G values  have  been  computed  for  the  value  of  |H| 
nodes,  the  data  may  be  retained  and  used  to  provide  G values  for 
any  job  load  vector  where  the  individual  job  type  values  are 
less  than  or  equal  to  those  initially  computed  for  G. 
Consequently,  for  the  second  and  subsequent  passes  of  step  3, 
the  G values  need  not  be  recomputed.  An  interactive  program 
(TPOT)  at  Appendix  D provides  such  a capability.  Also,  if 
memory  size  constraints  are  not  inhibiting,  then  it  is  possible 
to  compute  G values  for  the  maximum  anticipated  load  vector  and 
retain  these  values  in  a data  file,  thus  eliminating  the  need 
for  further  G computations. 
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Verification  Using  GESS 

I 

The  saae  GPSS  results  as  presented  in  Chapter  7 were 
I compared  against  our  integer  iterative  model.  The  Network 

Description  found  in  Chapter  ITT  and  used  for  the  Basic  Model 
and  the  balanced  service  rates  of  Chapter  IV  are  used  in  this 
verification. 


Computation  for  NA=6,  NB=2,  NC=ND=1,  and  N=10 


Iteration 

Number 

Computed  Thruput 
[Tkfnk)  ] 

Proportion 
Types  F(k) 

of  Tr[  N(k)  ] 

Tr  (N) 

1.  318 

0.883 

1.08 

1 

0.099 

0.066 

19.  1 

0.031 

0.021 

32.3 

(6,2,1,  1) 

0.044 

C.03C 

22.7 

Computed 

job 

load  vector 

= [0.402 

1.121  0.683 

0.602] 

1.290 

0.873 

1.00 

3.22 

2 

0.086 

0.058 

16.7 

0.039 

0.027 

25.6 

0.063 

0.043 

15.  9 

Computed 

jcb 

load  vector 

= [0.377 

1.054  0.631 

0.514] 

Table  VI-1 


r 


I 

1 

t 
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Coaputaticn  fcr  NA= 

12,  NE=4,  NC=ND=2, 

and  N=20 

t 

Iteration 

Computed  Thruput 

Proport icn 

of 

Tt[N(k)  ] 

Tr(N) 

Number 

[Tk(n)c)  ] 

Types  p(k) 

1.243 

0.855 

1.91 

1 

0.117 

0.081 

37.  1 

0.040 

0.028 

48.  9 

(■>2, 4, 2, 2) 

0.054 

0.037 

35.  9 

' 

Computed 

job 

load  vector 

= [1.357 

2.715 

1.  530 

1.410  ] 

1.165 

0.823 

1.79 

2 

0.  147 

0.  104 

19.4 

0.053 

0.038 

36.  0 

(1,3,2, 1) 

0.051 

0.036 

30.7 

Compu  ted 

job 

load  vector 

= [1.282 

2.253 

1.  142 

1.348] 

1.210 

C.841 

1.68 

5.44 

3 

0.133 

0.093 

19.7 

0.035 

0.024 

47.  3 

(1,2, 1,1) 

C.060 

0.042 

25.4 

Computed 

job 

load  vector 

= [1.206 

2.  270 

1.519 

1.258] 

1.206 

0.839 

1.69 

5.73 

4 

0.  120 

0.084 

22.5 

0.059 

0.041 

32.  1 

(1,2,2, 1) 

0.051 

C.036 

30. 4 

Computed 

job 

load  vector 

= [1.213 

2.400 

1.  363 

1 .344  ] 

Final  result 

0.840 

1.68 

5.58 

0.088 

21.1 

0.032 

39.7 

0.039 

28.  1 

Table  VI-2 
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In  Table  VI-2  the  exact  value  for  convergence  lies 
somewhere  between  the  values  obtained  in  iterations  three  and 
four.  Conseguently , there  is  an  oscillation  between  the  two 
sets  of  values.  This  situation  is  due  to  the  restriction  of 
using  only  integer  values  for  the  job  lead  vector.  A judicious 
analysis  is  made  to  obtain: 

E = [.840  .008  .032  .039] 

Tr  = [1.68  21.1  39.7  28.1] 

and 


Tr  (N)  = 5.58  seconds 

as  the  final  results  for  twenty  jobs. 


Chapter  VI 


1 17 


Computaticn  fcr  NH»18, 

NB=6,  NC=ND=3,  and  N=30 

Iteration 

Coaputed  Thruput 

Proportion 

of  Tr[N(lc)  ] 

Tr  (N) 

Nunber 

[Tk(nlc)  3 

Types  p(k) 

1.206 

C.724 

3.78 

1 

0.126 

0.  151 

44.  4 

0.046 

0.055 

64.  1 

(18,6,3,3) 

0.059 

0.070 

49.  4 

Conpu  ted 

job 

load  vector  = 

[3.626 

4.484  2.431 

2.301  ] 

1.118 

0.803 

3.61 

2 

0.  166 

C.  119 

27.  4 

0.040 

0.029 

64.  4 

(“,5,2,2) 

0.068 

0.049 

35.0 

Conpu  ted 

job 

load  vector  = 

[ 3.494 

3.876  2.434 

2.101 3 

1.  128 

0.807 

3.51 

9.43 

3 

0.154 

0.  110 

27.  8 

0.043 

0.031 

59.4 

(3,4,2, 2) 

0.042 

0.052 

32.7 

Computed 

job 

load  vector  = 

[3.413 

3.895  2.395 

2.057  3 

Table  VI-3 
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The  comparison  of  this  model  and  GPSS  can  be  made  not  only 


on  the  final  composite  response  time,  Tr(N),  but  also  with  the 
individual  job  type  or  route  response  times,  Tr)c[n)c].  The 
proportionality  of  job  types  in  the  CEN  for  this  model  is  also 
available  and  is  compared  to  that  computed  by  GPSS.  Comparisons 
of  these  values  yield  further  insight  into  the  interaction  or 
interference  of  the  job  types  within  the  computer  network. 


Comparison  of 
Iterative  flcdel  and  GPSS 

N £ Tr  Tr(N) 


Iter 

GPSS 

Iter 

GPSS 

Iter 

GPSS 

10 

.873 

.729 

1.00 

.89  1 

3.  22 

4.32 

.058 

. 150 

16.7 

10.7 

.027 

.053 

25.6 

21.3 

.043 

.067 

15.9 

13.7 

error  in 

Tr(II) 

25  % 

20 

.840 

.767 

1.68 

1.53  5.  58  6.1  3 

.088 

.129 

21 . 1 

17.7 

.032 

. 044 

39.7 

32.9 

.039 

. 059 

28.1 

20.  3 

error  in  Tr(N)  = 9.0  % 

30 

.807 

.793 

3.51 

3.16 

9.43  9.02 

.110 

.113 

27.8 

27.  2 

.031 

. 038 

59.4 

47.  8 

.052 

.054 

32.7 

29.  3 

error  in  Tr(N)  = 4.5  % 


j. 


i 

t 

( 


I 


L 


Table  VI-4 
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A review  of  the  percent  error  at  each  of  the  comparisons 
reveals  two  major  items.  First,  the  model  provides  quite 
comparable  results  to  the  Basic  (lodel  with  sufficient  accuracy 
to  offer  somewhat  better  overall  insights  (see  Table  7-1) . 
Secondly,  as  the  total  number  of  jobs  is  increased,  the 
Iterative  Model  yields  increasingly  tetter  results.  This 
situation  is  the  reverse  of  that  given  by  the  Basic  Model.  The 
fact  that  larger  values  of  N lead  to  better  results  is  expected. 
The  error  introduced  due  to  integer  representation  of  the  jot 
load  vector  tecomes  less  pronounced  as  N is  increased. 


A 
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Sunnary  and  Motivation  for  Further  Refioesent 

The  four  job  type  raodel  using  integer  iteration  (IH) 
provides  a good  aodel  and  meets  out  objective  of  obtaining 
insights  concerning  response  times,  the  major  performance 
criterion.  In  particular,  it  complements  the  Basic  Model  by 
providing  greater  accuracies  with  large  N.  However,  we  do  incur 
a greater  complexity  of  computation. 

A useful  result  of  this  acdel  formulation  has  been  the 
development  of  the  specialized  programs  for  the  multijob  type 
normalization  constant  computation  based  on  the  GCT  (Appendix  C) 
and  the  solution  of  the  short  circuit  thruput  values  (Appendix 
D)  using  the  Norton  equivalent.  Both  of  these  programs  have 
been  developed  to  operate  in  an  interactive  mode  using  the  IBM 
ISO.  Although  the  Duke  ASQ  is  used  for  the  evaluations  in  step 
U of  the  algorithm,  we  have  the  basis  for  our  own  program 
solution  since  the  computation  of  G is  readily  available. 

Although  we  have  obtained  a more  accurate  model  for  larger 
values  of  N,  we  are  motivated  to  develop  a still  more  accurate 
model  throughout  the  entire  range  of  N.  Based  on  the  knowledge 
gained  by  development  of  the  Iterative  Model,  it  appears  that 
such  a new  model  may  be  at  hand. 
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Chapter  VII:  EIACI  HNCFR  MODEL 


General  Approach 


In  the  fornulatioD  of  the  Iterative  Model  (Chapter  VI)  we 
have  reduced  our  network  to  a two  ncde  representation  of  the 
TEHM  and  CPN,  Figure  Vll-1  (same  as  Figure  Vl-6) . 


Figure  Vll-1 

We  have  regarded  the  two  node  representation  as  a closed  system 
of  two  load-dependent  devices.  For  the  moment,  let  us  suppose 
that  we  are  dealing  with  one  joh  class  and  a total  of  N jobs  and 
let  n be  the  number  of  jobs  in  the  CPN.  For  this  specific  case 
we  find  that  Buzen  and  Denning  [1977]  have  provided  a means  of 
obtaining  exact  solutions.  The  solutions  provided  yield  a 
direct  computation  of  thruput  as  a function  of  N,  T(N)  and  queue 
lengths.  In  our  specific  case,  the  queue  length  we  are 
interested  in  is  the  average  number  of  jcbs,  n,  in  the  CPN. 

T(N)  = [ (N  / Z)  h(N-1,  m)  ] / h(N,  B)  (VIl-1) 


n 


N [ 1 - (h  (N-1,  m))  / h (N,m)  ] 


(VII-2) 


All  parameters  are  as  previously  defined  except  for  the  h,  which 
is  the  value  obtained  in  the  computation  of  a special 
normalizing  constant,  H,  which  in  tern  is  related  to  our 
familiar  G normalizing  constant  of  Chapter  VI.  Being  able  to 
compute  h permits  us  to  solve  equations  VII-1  and  VII-2. 
Furthermore,  we  recognize  that  the  mean  response  time  (time 
spent  in  the  CPN)  can  then  be  simply  computed  [Buzen  76a]  as: 

Tr  (N)  = n / T (VII-3) 


The  H relationship  is  established  as 


H = 


N! 

2^ 

^ (N-n)  ! 2 ♦*  n 


(VII-U) 


which  leads  to  an  iterative  solution  technique  analogous  to  the 
G computation.  The  recursive  specification  for  h is 


h{a,  m)  = h(a,  m-1)  ♦ a / 2 B(m)  h(a-l,  m)  (VII-5a) 

for  a > 0 and  m > 0 


with  the  following  initial  conditions: 


h(0,  m)  = 1 

for  ■ > 0 

(Vll-Sb) 

h (a,  0)  = 1 

for  a > 0 

(VII-Sc) 

Here,  the  a is  a running  variable  representing  the  number  of  the 
maximum  N jobs  in  the  two  node  network  representation  of  Figure 
VTI-1.  Be  immediately  recognize  that  the  recursive 


I 


specification  of  Vli-5  is  the  same  as  that  for  g,  equation 
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VI-15,  with  the  exception  of  the  nultipl icative  term,  a / Z,  for 
the  relative  utilization,  H(ni).  It  is  this  multiplicative  term 
which  accciints  for  the  think  time  of  jots  at  TERM. 

We  have  already  established  a means  of  computation  for  the 
G normalization  constant  in  Chapter  IV  and  can  simply  use  the 
same  technique  for  the  solution  of  H.  In  fact,  if  we  take  an 
engineering  approach  tc  the  program,  we  can  by  direct  analogy 
obtain  the  solution  form  of  H for  multiple  jot  classes.  We  have 
already  noted  that  the  sole  difference  between  g and  h is  the 
multiplication  of  the  relative  utilization  by  the  term  a / Z. 
Equations  VI-16  and  VI-17  show  the  extension  of  g from  the 
single  job  type  case  to  that  for  two  and  four  job  types, 
respectively.  Recall  that  a is  the  running  variable  for  the 
number  of  type  A jobs  in  the  network  and  zA  is  the  think  time  of 
job  type  A.  Then  the  multiplicative  term  in  the  H normalization 
computation  for  four  job  types  is: 


W(m) 

becomes 

a 

/ 

zA 

* 

s (a) 

X (m) 

becomes 

b 

/ 

zB 

X{m) 

Y(m) 

b ecomes 

c 

/ 

ZC 

♦ 

T{m) 

Z (1) 

becomes 

d 

/ 

zD 

Z (m) 

Note  that  we  have  the  ability  to  consider  individual  think  times 
for  each  job  route,  not  possible  with  the  Basic  Model.  Thus  h 
for  four  job  types  is 
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h(nA,  cE,  nC,  dD,  m)  - h{nA,  dE,  nC,  nD,  u-l) 


♦ 

a 

/ 

2 A 

♦ 

W (n) 

h (n  A- 

1 , nB,  nC , nD, 

B) 

+ 

b 

/ 

2E 

♦ 

X (ni) 

h (n  A, 

nB-T,  nC , nC, 

B) 

+ 

c 

/ 

2C 

Y (ni) 

h (nA, 

dE,  nC-1 , nD, 

B) 

■f 

d 

/ 

zC 

Z (in) 

h (n  A , 

nB,  nC,  nD-1, 

») 

{VII-6) 
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(lodel  PocBulation  and  I apleiaentaticn 

Due  to  our  previous  work  with  the  iterative  model  of 
Chapter  VI,  the  HNORB  Bodel  formulation  and  implementation  is  a 
natural  result.  Be  have  already  established  the  technique  of 
iterative  computation  of  g for  nultifle  job  classes  which  was 
shown  as  a GCT.  The  same  computational  tree  structure  for  the 
recursive  specification  of  h is  also  applicable  here.  Now, 
having  h for  multiple  job  classes,  we  can  extend  the  equations 
for  thruput,  average  number  of  jobs  in  the  CPN,  and  mean 
response  time  for  multiple  job  types  to  obtain 


T(NA) 

= [NX  / ZA  H(NX-1,  NB,  NC,  ND)  ] 

/ H{NA,  HB,  NC,  ND) 

(VII-7) 

nX  = 

NX  [1  - H(NA-1,  NB,  NC,  ND)  ] 

/ H{NA,  NB,  NC,  ND) 

(VII-B) 

Tr  (NX) 

= nA  / T (NX) 

(VII-9) 

and  so  forth  for  job  types  B,  C,  and  D. 

We  establish  the  two  node  representation  of  the  network  as 
in  Figure  VII-1  and  then  solve  for  the  parameters  of  interest. 
We  may  additionally  obtain  the  proportion  of  jobs  of  a 
particular  type,  p(k),  that  are  in  the  CPN  by  taking  the 
proportion  of  thruputs  as  was  done  in  equation  VI-20: 

Tk  (Nk) 

2 Tk(Nk) 

k 


P(^) 


f k 


(VII-1G) 
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As  was  done  in  the  previous  model  (IJI),  the  results 
obtained  using  the  H normalization  procedure  were  verified 
against  test  runs  on  a more  limited  basis. 

1.  The  single  job  type  example  provided  by  Buzen  and 
Denning  [ 1977  ]. 

2.  A simple  two  node,  two  job  type  case  compiled  by  hand. 

3.  A three  job  type  computation  run  at  Federal  Simulations 
where  the  model  of  Chapter  III  was  modified  by  the 
elimination  of  r(4)  and  g (3)  was  increased  from  0.1  to 

0.  2. 


The  results  obtained  were  in  all  cases  exactly  the  same  using 
this  model  with  the  H normalization. 


For  the  network  of  functionalized  processors  using  the 
balanced  service  rates  as  before,  the  Exact  HNCRf!  Model  gave  the 
following  results: 

Exact  H Model  Results 

£ ll 


10 

20 


.731 
. 1«8 
.053 
.068 

.766 
. 130 
.044 
.060 

.793 

.115 

.038 

.053 


.907 

11.1 

21.6 

13.6 

1.60 

17.6 
33.2 

20.6 

3.  22 
26.8 
47.9 
30.1 


Tr  (N) 
4.38 

6.20 

9.08 


Figure  VII-2 


Chapter  VII 


127 


Although  only  the  job  points  of  10,  20,  and  30  are  shown,  all 
possible  values  for  job  types  nay  be  readily  coaputed  since  all 
H values  are  ascertained.  In  Appendix  E is  a program  to  compute 
the  H values  and  the  parameter  values  of  interest  (thruput, 
number  of  jobs  in  the  CPN,  and  the  mean  response  times) . These 
latter  parameter  values  are  for  NA,  MB,  NC,  and  ND.  The  values 
of  interest  may  be  readily  obtained  for  any  other  network  load 
(nA,  nB,  nc,  nD)  where  nk  < Nk  f k once  the  H values  are 
obtained.  In  Appendix  F there  is  an  interactive  program  using 
the  IBB  TSO  which  takes  the  H values  as  a file  input  from  the 
Exact  RKORM  Model  Computation  and  computes  the  values  of 
interest  for  any  network  load  as  specified  above. 
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Verification  Osin'g  GESS 


Following  the  procedure  previously  established,  we  have 
compared  the  HNCRIl  flodel  (HH)  results  to  the  GPSS  simulation. 
At  Table  VII-2  are  the  comparative  results.  It  is  evident  that 
this  model  representation  is,  in  fact,  an  exact  computation  when 
we  acknowledge  that  the  event  driven  simulation  model  will  yield 
accuracies  of  5-10  percent  of  the  actual  system. 


Comparison  of 
Exact  H flodel  and  GFSS 

N £ Tr  Tr(N) 


HH 

GPSS 

HH 

GPSS 

HH  GPSS 

10 

.731 

. 729 

.907 

.891 

4.  38  4.32 

. 148 

. 150 

11.1 

10.7 

.053 

. 053 

21.6 

21.3 

.068 

.067 

13.6 

13.  7 

error  in  Tr  (R) 

= 4.0  % 

20 

.766 

.767 

1.60 

1.53  6.  20  6.1  3 

.130 

.129 

17.6 

17.7 

.044 

. 044 

33.2 

32.9 

.060 

. 059 

20.6 

20.  3 

error  in  Tr(N)  < 1.0  % 

30 

.793 

.793 

3.22 

3.  1 6 

9.  C8 

9.02 

.115 

.113 

26.8 

27.  2 

.038 

. 038 

47.9 

47.  8 

.053 

. 054 

30.  1 

29.  3 

error  in  Ir(N)  < 1.0  % 


I 


I 


Table  VII-2 
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Susmaty  of  Exact  HSOEH  Model 

ie  have  called  this  the  Exact  HNOBB  Model  because  within 
the  assuaptions  of  EQR  aodelling  we  have  obtained  exact 
solutions  for  the  paraaeters  using  the  H nornalization  constant. 
This  model  can  he  easily  used  and  provides  the  performance 
criteria  in  which  we  are  interested.  In  fact  it  provides  more 
flexibility  than  we  have  had  before.  Not  only  do  we  obtain 
individual  route  thruputs  and  response  times,  but  we  can  easily 
designate  individual  route  think  tines.  This  capability  adds 
another  degree  of  flexibility  to  our  network  aodelling. 

Ose  of  the  programs  developed  at  Appendices  E and  F provide 
a means  of  ready  model  solution  and  network  analysis. 
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Chapter  viii:  CONCLDSIONS 


General  Conclusions 


There  are  several  conclusions  to  be  drawn  from  this 
research  concerning  the  development  of  the  three  EQN  models 
presented  and  from  the  models  themselves.  First,  recall  that 
the  basic  objective  was  to  develop  a model  that  could  be  used 
for  trade-off  analysis  of  large  scale  computer  networks  and  not 
just  "local"  computer  systems.  Secondly,  we  wanted  a design  aid 
to  guide  in  the  solution  of  service  rates  for  the  processors  in 
the  network.  It  was  felt  that  these  objectives  could  be 
obtained  as  long  as  we  were  willing  to  accept  insights  and  not 
concentrate  on  details,  select  response  time  as  a major 
performance  criterion,  and  constrain  the  network  model  in  terms 
of  functionalized  processors.  Furthermore,  by  the  use  of 
functionalization,  one  is  able  to  impart  generality  to  the  model 
in  the  sense  that  the  network  may  be  viewed  as  being  composed  of 
several  defined  functions.  Since  the  model  is  constructed  of 
functionalized  processors,  it  is  then  the  responsibility  of  the 
designer/analyst  to  ensure  that  the  characterization  of 
functions  in  the  computer  network  is  an  accurate  network 
representation. 

Through  the  use  of  the  exponential  queuing  network 
approach,  we  have  seen  three  methods  of  solving  the  problem. 
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These  Methods  have  been  presented  to  accoDoodate  a large  nanber 
of  jobs  and  nodes  so  that  an  effective,  total  conputer  network 
say  be  Modelled. 

In  the  developnent  of  the  models  te  have  seen  that  the 
General  Purpose  EQN  solution  programs  such  as  the  versions  of 
ASQ  and  QNET4  can  work  well  with  small  networks,  few  job  classes 
and  few  total  jobs.  They  sees  to  be  best  suited  for  the  network 
representation  of  a local  computer  system.  When  we  wish  to 
model  on  a larger  scale  it  is  necessary  to  develop  more  specific 
program  models.  Cur  approach  through  the  use  of  the  computer 
program  nodules  of  Appendices  A - F permits  the  representation 
and  solution  of  such  models,  in  fact,  with  modification  and 
extension  of  these  modules,  a general  purpose  system  could  be 
developed.  Furthermore,  the  nodules  provide  the  job  route 
(type)  response  times  of  the  interactive  system  being  modelled. 
Of  the  other  known  solution  programs  discussed  in  Chapter  III, 
such  a capability  is  not  available.  In  fact,  the  HNOBH  model 
provides  the  response  times  directly  for  an  interactive  system 
with  multiple  job  types. 
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nodel  Pobastness 

With  the  use  of  ECN  theory,  we  have  restricted  ourselves  to 
the  assuBiption  of  exponential  service  distributions  with  FCFS 
service  disciplines  or  to  other  service  distributions  with  a 
rational  Laplace  transfora  and  service  disciplines  of  processor 
sharing  (PS)  or  last-coae-f irst-served  preemptive-resuae 
(LCFSPP)  as  discussed  by  Chandy  et  al  [1972],  This  restriction 
has  made  it  possible  to  obtain  solution  forms  for  the  EQN. 
However,  these  assumptions  are  not  in  the  least  way  restrictive. 
The  EQN  models  developed  are  extremely  robust,  particularly  when 
we  consider  that  we  are  not  searching  for  details  of  network 
performance  but  for  insights  and  trends.  Other  distributions 
may  in  fact  be  utilized.  To  give  a specific  example,  different 
service  distributions  (other  than  exponential)  were  selected  and 
used  in  the  GPSS  simulation.  Recall  that  selection  of  service 
distributions  is  one  of  the  advantages  of  the  GPSS.  The  PIP  and 
COnP  were  assigned  hypo-exponential  distributions  of  degree  two 
[Hellerman  and  Conroy  7*5],  expressed  as  a function  of  time, 
which  is: 

S(t)  = [1  ♦ 2 (1  / u(i))  t]  **  [-2  (1  / u(i))  t]  (VIII-1) 

where  S (t)  = pr[  no  job  copletes  by  time  t] 
and  1 - S(t)  = Pr[a  job  completes  by  time  t] 

for  i = F,  C.  Note,  however,  that  the  same  mean  service  rates 
(u(i))  are  to  take  on  the  same  values  used  for  the  EQN  models. 
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Such  a distribation  is  one  that  is  expected  for  the  FIP  and  COHP 
functions  in  an  interactive  environment.  The  DBH  and  PEP  were 
assigned  hyper-exponential  distributions  for  their  expected 
processing  function  characteristics.  Expressed  as  a function  of 
time  as  was  equation  VIIl-1: 

S(t)  = c ♦♦  [-2  c (1  / u(i))  t] 

♦ (1  - c)  *♦  r-2  (1  - c)  (1  / u(i))  t]  (VIII-2) 

where  c = .45  for  i = D,  P.  Again  note  that  the  u(i)  are  the 
mean  service  rate  values  used  in  the  EQN  models.  Table  VIII- 1 
compares  the  use  of  balanced  u(i}  with  exponential  distributions 
to  the  same  values  of  balanced  u(i)  with  hypo-  and  hyper- 
exponential  distributions.  The  6PSS  simulations  produce  nearly 
identical  results.  It  should  also  be  noted  that  the  service 
disciplines  have  all  been  FCFS.  Such  a combination  of  service 
disciplines  and  distributions  is  a clear  violation  of  the  EQN 
assumptions.  Therefore,  this  example  supports  the  conclusion 
that  the  solutions  obtained  using  the  EQN  models  will  still 
yield  excellent  comparative  results.  Such  comparisons  and 
analysis  are  further  enumerated  in  Buxen  [ 1975]  and  Giammo 
[1976].  In  our  specific  example,  the  performance  criterion  is 
response  tine,  which  is  known  to  provide  even  less  favorable 
comparisons  than  the  performance  criteria  of  utilization  and 
thrupnt. 
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Results  of  Hyper-,  Hypo-Exponential  Distributions 

GPSS  Sifflulation 


N 

Distribution  Proportion  of 

Types  p(k) 

Tr[N  (k)  ] 

(sec) 

Tr  (N) 
(sec) 

10 

Exp 

C .729, .150,. 053, .067] 

[ .891,10.7,21.3,13.7] 

4.32 

Hyper, 

Hypo 

[ .731, .147, .052,  .068  ] 

[ . 754,11.  1,2  1.7,  12.9  ] 

4.21 

% error 

= 2.6 

20 

Exp 

[ .767,.  129,  .044,  .059  ] 

[ 1.53,17.7,32.9,20.3  ] 

6.13 

Hyper, 

Hypo 

[ .763, .130, .044, .060  ] 

[ 1.23,16.6,31.3,18.9  ] 

5.65 

% error 

= 8.5 

30 

Exp 

[ .793,. 1 13, .038, .054] 

[ 3.  16,27.2,47.8,29.3  ] 

9.02 

Hyper, 

Hypo 

[.790, .117, .038, .053] 

[2.69,24.8,46.0,28.3 ] 

8.  32 

% error 

= 8.4 

Note:  The  mean  values  of  the  distribution  were  those  obtained 

for  a balanced  systen  (Chapter  IV)  . 
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Hodel  Conpacisons 

This  section  contains  the  results  of  each  of  the  lodels 
developed,  Basic  Hodel  (BH) , Iterative  Hodel  (IH)  and  the  Exact 
HNCRH  Hodel  (HH) , as  well  as  the  GPSS  siaulation  verifications 
(S) . Results  are  presented  for  two  different  sets  of  service 
rates.  The  balanced  rates  and  those  rates  referred  to  as  FEDSIH 
(used  in  the  validation  of  the  G Noraalization  constant  at 
Federal  Siaulations  as  noted  in  Chapter  VI) . Both  sets  of  rates 
are  given  for  conparative  purposes.  Also,  the  GPSS  run  using 
the  Bean  service  rates  for  the  balanced  systen  tut  with  hypo- 
and  hyper-exponential  distributions  is  presented.  The  aodel 
designations  are  shown  as  a function  of  E,  F and  H depending  on 
the  use  of  balanced  or  Fedsiii  cates  with  exponential 
distributions  or  balanced  cates  with  the  hyper-,  hypo- 
exponential  distributions,  in  all  cases  the  zb  = 15  sec  and  the 
job  load  distribution  was  taken  such  that  HI  / N = 0.6,  112  / N = 
0.2,  H3  / N = 0.1  and  Nft  / N = 0.1  as  was  established  in  the 
paraaeter  definition  of  the  Basic  Hodel  and  used  for  the 
Iterative  and  HNOFH  aodels. 
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Conparative  Results 
N=10,  N1=6,  N2=2,  N3=N4=1 


Model 

Proportion  of 
Types  p()c) 

Tr[N  (k)  ] 

(sec) 

Tr  (N) 
(sec) 

* 

err 

BH  (B) 

C“  » 

4.  50 

4.2 

IH(B) 

[ .873,.058,.027,.OU3  ] 

[ 1.00,16.7,25.6,15.9] 

3.  22 

25 

HM  (B) 

[ .731, .148, .053, .068  ] 

[ .907,  11  . 1,21.6,13.6] 

4.  38 

1.4 

S(B) 

[ .729,.  150, .053, .067  ] 

[.891,10.7,21.3,13.7] 

4.32 

— 

S(H) 

[ .731, .147, .052,  .068  ] 

[.754,11.1,21.7,12.9] 

4.  21 

2.5 

BH  (F) 

[ » » » ” ] 

* 9 **  9 

--  ] 

4.  23 

2.7 

TH(F) 

[ .861, .071, . 028, .040  ] 

[ 1.03,13 

.4,24.4, 

17.2] 

3.  21 

22 

HH  (?) 

[ .722,.  159,. 054, . 065  ] 

[ .915,9. 

01,20.6, 

14.5] 

4.14 

0.5 

S(F) 

[ .721  ,.  159,  .054,  .064  ] 

[ .913,8. 

92,20.2, 

14.5] 

4. 12 

-- 

Table  VIII-2 
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Coaparative  Results 


N=2C,  B1= 

12,  N2=4,  N3=N4=2 

Model 

Proportion  of 

Tr[N(k)  ] 

Tr  (N) 

% 

Types  p(k) 

(sec) 

(sec) 

err 

BH  (B) 

] 

[ " » " » “ • ] 

7.  64 

25 

in  (B) 

[.840, .088, .832, .039] 

[ 1.68,21. 1,39.7,28.1] 

5.  58 

9.0 

HH  (B) 

[ .766,.  130, .044, .060  ] 

[ 1.60, 17. t, 33. 2, 20. 6] 

6.20 

1.1 

S(B) 

[ .767,.  129, .044,.  059  ] 

[ 1,53,17.7,32.9,20.3] 

6.  13 

-- 

S(H) 

[ .763, .130, .044, .060  ] 

[ 1.23,16.  €,31.3,18.9] 

5.65 

7.8 

BH  (F) 

[ • » ~ ] 

[ — . 

7.  26 

22 

in(F) 

[ .819,.  116, .026, .039] 

[ 1.82, 

15.1,44.4,28.1] 

5.  49 

7.6 

HH  (P) 

[.753, .146, .045, .055] 

[ 1.66, 

13.6,31.5,22.8] 

5.  92 

0.3 

S(P) 

[ .752, .146, .045,  .055  ] 

[ 1.69, 

13.6,31.2,22.8] 

5.  94 

-- 

Table  VIII-3 
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Ccmparative  Results 
N = 3C,  N1  = 18,  »I2=6,  N3  = N4  = 3 


Model 

Proportion  of 

Tr[N  (k)  ] 

Tr  (N) 

Types  p(K) 

(sec) 

(sec) 

err 

Bfl  (B) 

c ] 

[ » f » ~ 3 

12.7 

41 

in  (B) 

[ .807,. 110, .031, .052] 

[ 3.51,27.8,59.4,32.7] 

9.  43 

4.5 

Hn  (B) 

[ .793,.  1 15, .038,  .053  ] 

[ 3. 22, 26. £,47.  9, 30. 1 ] 

9.  08 

0.7 

S(B) 

[ .793,.  1 13, .038, .054  ] 

[ 3.  16,27.2,47.8,29.3] 

9.  02 

— 

S(H) 

[ .790,.  1 17  ,.038, .053  ] 

[2. £9, 24. 8,46. 0,28. 3] 

8.  32 

7.  e 

BK  (F) 

12.6 

42 

IP!  (F) 

c .794,.  126, .040, .055  ] 

[ 3.80,23.8 

,53.4,37.5] 

7.68 

13 

HH  (F) 

[ .778,.  1 34, .040,.  049  ] 

[ 3.48,20.8 

,45.6,34.2] 

8.96 

1.  1 

S(P) 

[ .778, , 133, .038, .048  ] 

[ 3. 38,20.8 

,46.4,33.8] 

8.86 

-- 

Table  VIII-U 
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As  the  Bodels  were  developed,  a ccnparison  to  GPSS  for 
verification  was  Bade  and  a soall  discussion  concerning  the 
model’s  effectiveness  was  presented.  Here  we  will  sake  a 
comparison  of  the  nodels  and  their  effectiveness  amongst 
themselves.  As  already  noted,  there  is  a definite  region  of 
acceptable  effectiveness.  Figure  VIII-1  is  a diagram  indicating 
this  fact. 


Figure  VlTl-1 


This  figure  illustrates  that  the  Easic  Kodel  (Bfl)  will  produce 
effective  results  when  loaded  below  the  saturation  value,  N*. 
The  Iterative  Hodel  (IB)  yields  the  best  results  when  there  is  a 
large  job  lead.  Depending  on  the  network,  this  region  may  well 
include  job  loads  below  the  designed  saturation.  For  instance, 
a network  with  a high  W*  lay  be  modelled  by  the  in.  The  Exact  H 
nodel  (HH)  perforas  extreaely  well  across  the  entire  spectrum. 
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RecoaBended  nse  of  the  Models 

Based  on  our  previous  discussion  and  Figure  VIII-1  the 
following  points  are  noted: 

1.  The  Basic  Model  Focaulaticn  provides  the  basis  for  a 
balanced  design  criterion.  Fven  though  its  response 
tine  results  could  exceed  a designer’s  desired  error 
tolerance  level  when  more  than  N*  users  are  active,  the 
utilization  results  are  sufficiently  accurate  to 
provide  good  design  guidance. 

2.  When  a quick,  rough  analysis  is  desired,  the  Basic 

Model  is  the  best  approach.  This  assumes  that  a 
General  Purpose  Program  such  as  ASQ  is  available  or  one 
has  been  built  from  the  basic  program  modules 

(A  ppendices  A-F)  . 

3.  The  Iterative  Model  (IM)  can  be  used  to  analyze 

interactions  between  job  types  and  their  loading.  By 

itself,  it  is  not  a particularly  useful  model. 

However,  the  G Computation  Program  established  for  its 
use  is  quite  beneficial  for  the  development  of  other 
program  solution  models. 

4.  The  Bxact  HMOPM  Model  (HM)  provides  more  information 

and  is  more  accurate  in  all  regions  than  the  BH  but 
requires  mote  computer  execution  time  and  storage 
requirements  than  the  BH.  The  HM  is  also  more 
versatile  since  one  can  consider  individual  job  route 


1 

Chapter  VIII 

141 

think  tines  and  it  yields  job  route  mean  response  tines 
rather  than  just  the  composite  response  tine.  j 


Chapter  IX:  SUHKSBY  AND  SOGGESTICNS  POP  FDTOEE  WO 3K 


Based  on  the  techniques  of  Network  Queuing  Theory,  three 
Exponential  Queuing  Network  (EQN)  models  for  the  analysis  and 
design  of  projected  computer  networks  have  been  developed.  The 
projected  computer  network  must  be  parameterized  in  terms  of 
functions  in  order  to  be  modelled.  Such  functionalization  is  a 
natural  result  of  the  computer  network  structure  and 
organization  as  noted  in  the  examination  of  major  networks. 
Functionalization  such  as  computation,  data  base  management, 
communications  and  terminal/access  interface  have  already  been 
considered  in  major  military  command  and  control  systems.  The 
EQN  models  of  computer  networks  permit  trade-off  analysis  of 
various  network  structures.  Workload  description  is 
accommodated  by  permitting  route  definitions  or  circuit  paths 
through  the  network  which  are  composed  of  the  defined  processing 
functions.  As  a consequence,  jot  typing  or  classification  is 
implemented  in  the  models.  Additionally,  a basic  design 
technique  is  provided  for  the  synthesis/determination  of 
processing  "computing  power"  to  ensure  optimal  computer  network 
load  balance.  Based  on  the  major  performance  criterion  of 
response  time,  a thorough  evaluation  or  trade-off  analysis  of 
various  computer  network  organizations  and  effectiveness  may  be 
conducted.  The  networks  are  considered  to  include  all  possible 
functions  in  an  interactive  mode.  These  models  provide  trends 
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and  insights  into  the  perfornance  areas.  They  are  not  intended 
to  yield  detailed  inforiation  concerning  perfornance,  nor  is  it 
expected  that  the  inforaation  available  concerning  the  computer 
networks  to  be  aodelled  will  be  available  in  detail. 

A method  for  the  computer  network  definition  and 
description  is  presented  in  Chapter  III.  Each  of  the  three  EQS 
models  developed,  the  Basic  Hodel  (EH)  of  Chapter  IV,  the 
Iterative  Hodel  (IH)  of  Chapter  VI,  and  the  Exact  HNOBH  Hodel 
(HM)  of  Chapter  VII  employ  the  sate  methods  of  network 
definition  and  description.  In  Chapter  v a classical  event 
driven  simulation  programmed  in  GPSS  has  been  developed  in  order 
to  verify  the  models.  This  GPSS  verification  is  necessary  since 
we  have  developed  models  to  represent  projected  computer 
networks  which  may  be  structured  as  a network  of  functionalized 
processors.  Based  on  the  verification  using  GPSS,  the 
advantages,  disadvantages,  and  the  recommended  uses  for  the 
models  have  been  presented.  The  BH  can  provide  quick,  overall 
insights  yielding  its  best  results  when  the  network  is  operating 
below  saturation.  The  IH  provides  more  detailed  information  and 
insights  concerning  job  type  or  class  interaction,  but  this  is 
done  at  the  expense  of  increased  computation.  The  HH  provides 
information  on  all  job  types  and  is  accurate  throughout  the 
range  of  network  loading. 

In  the  development  of  the  EH  and  IH,  a general  purpose 
single  job  type  program  analyzer.  Cuke  ASQ,  has  been  used. 


Chapter  IX 


104 


Additionally,  computer  program  mo 
models  are  provided.  These  prog 
to  tailor  the  model  solution  to  h 
the  basis  of  a general  purpose  mu 


There  are  three  areas  in  vhi 
needed.  The  first  area  to  cons 
models  for  the  trade-off  analysis 
in  the  military  command  and  cont 
we  have  visualized  a networX  as 
(Figure  Ix-1),  where  the  FIP  func 
security  or  access  control  a 
unintelligent  terminals. 


dules  for  the  solutions  of  the 
ram  modules  enable  the  designer 
is  application  or  they  can  form 
Iti-job  type  program  analyzer. 

ch  we  feel  additional  work  is 
ider  is  the  employment  of  these 
of  various  network  structures 
rol  environment.  For  instance, 
shown  in  the  diagram  below 
tion  includes  the  capability  of 
nd  TEBH  represents  rather 


Computer  Processing  Network 
Figure  IX-1 

'^n  the  other  hand,  we  could  consider  a network  structure  such  as 
the  one  in  Figure  lX-2  where  the  security  control  capability  has 
been  dispersed  throughout  the  other  elements  of  the  network  and 
the  remaining  capability  of  the  FIF  has  been  delegated  to 
intelligent  terminals. 
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Figure  II-2 

h coaparative  trade-off  analysis  could  then  be  aade  to  deteraine 
which  structure  perforaed  best  and  at  what  cost.  Although  the 
nuaber  of  different  scenarios  for  trade-off  analysis  is 
▼irtually  endless,  «e  are  able  to  selectively  provide  a 
quantitative  basis  upon  which  the  decision-aaker  nay  resolve 
questions  concerning  network  perforaance  and,  in  the  final 
analysis,  its  costs.  This  capability  is  particularly  geraaine 
in  the  area  of  ailitary  coaaand  and  ccntrol  where  one  needs  to 
know  in  advance  the  expected  response  tine  of  a network  under  a 
given  load  and  the  cost  of  achieving  this  response  tine. 


A second  area  for  future  work  is  in  the  refineaent  of  the 
progran  aodules  of  Appendices  A - F.  These  aodules  can  be 
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APPENDIX  A:  PROGBAH  FOR  SOLUTION  OF  B EQUATIONS  IN  N UNKNOWNS 


VISSOL:  PROC  OPTIONS  (SAIN)  ; 

DCL  HNODES  FIXED  BIN; 

DCL  {P(BNODES, HNODES)  ,PT(HNODES,BNCDES)  ) FLOAT  BIN  CTL; 

DCL  PAOG  (HNODES, HNODES)  FLOAT  BIN  CTl; 

DCL  TEHPCOL (HNODES)  FLOAT  BIN  CTL; 

DCl  VSOL  (HNODES-1)  FLOAT  BIN  CTL; 

DCL  (V (HNODES) , HU (HNODES) ,X (HNODES) ) FLOAT  BIN  CTL; 

/♦  P IS  THE  TRANSISTION  PROBAEILITT  HATHIX. 

PT  IS  THE  TRANSPOSE  OF  P.  P IS  SET  OP  WITH 
THE  BALANCE  EQUATIONS  SUCH  THAT  THE  V (I) 

FORH  A LEFT  EIGENVECTOR.  THE  V(I)  WILL 
EORH  A RIGHT  EIGENVECTOR  FOR  PT.  THE  PAUG 
IS  FOBBED  SUCH  THAT  THE  RIGHT  HOST  COLOHN  IS 
THE  V VECTOR  AND  IS  NOBHALIZED  50  THAT  V(T) 

IS  SET  TO  ONE. 

HNODES  CORRESPONDS  TO  !H!  ♦/ 


GET  LIST  (HNODES)  ; 

ALLOCATE  P, PT , PAUG , TEHPCOL , VSOL, V, HO,X; 
P=0.O;PT=0.0;PAOG=0.0;VSOl=0.0;V=0.C;HD=0.C;X=0.0; 
GET  LIST  (P)  ; 

POT  SKIP  (3)  DATA (HNODES)  ; 

POT  SKIP  (3)  DATA(P)  ; 

DO  1=1  TO  HNODES; 

PT  (I,*)  =P  (♦,!)  ; END; 

PAOG=PT; 

DO  1=1  TO  HNODES; 

PAOG(I,I)=PT(I,I)-1.0;  END; 

TEHPCOL  (*)  =PAOG  (♦,  1)  ; 

DO  JC=2  TO  HNODES; 

PAUG (*,30-1) =PADG(*, JO) ; END; 

PAOG  (•,HNODES)=TEHECCL (♦)  ; 

POT  SKIP  (3)  DATA  (PAOG); 

FREE  PT,  TEHPCOL; 

CALL  HATSOL  (H NODES, HNODES , PAUG ,VSO I ) ; 

V(1)=1.0; 

DO  1=2  TO  HNODES; 
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V (I)  =VSOL  (1-1)  ; END; 

/*  DETZRHINE  THE  X VALUES  FCF  NORM  CCHPOTATIGN  */ 

/*  NORMALIZE  ON  X (1)  */ 

GET  LIST  (BO) ; 

X(1)=V(1)/Bn(1)  ; 

DO  1=2  TO  BNODES; 

X(I)=V(I)/!10(I)/X(1)  ; 

END; 

X (1) =1.0; 

POT  SKIP  (3)  DATA  (V)  ; 

POT  SKIP  DATA  (BO)  ; 

POT  SKIP  DATA  (X)  ; 

/*  PROCEDORE  FOR  THE  SOLOTICN  OF  B EQOATICNS, 

N ONKNCWNS  */ 

y :tim******m*m***m****  ***************  ^ ******  *********/ 

BA'^SOL:  PROC  (H,N,XHAT,  BSOI)  ; 

DECLARE  (B,N)  FIXED  BINARY; 

DCL  XBAT(*,*)  FLOAT  BIN; 

DCL  HSOL(*)  FLCAT  BIN; 

DYNABIX; 

BEGIN; 


DECLARE  (SIDE  (B)  , TCP  (N-1))  CHAR  (1); 

DCL  BALPEA  CHAR(12)  INIT ( ' A BCDEPGH I J RI • ) 1 
DCL  EALPHA  CHAR  (12)  INIT ( 'ZYXiVOTSECFC • ) ; 

DO  1=1  TC  N-1; 

TOP  (I)  =S0ESTB  (EALPHA,N-I,  1)  ; END; 

DO  1=1  TO  B; 

SIDE (I) =SOBSTR  (BALPHA, I,  1) ; END; 

POT  PAGE; 

POT  SKIP(3)  LIST  (‘THE  INITIAL  BATFIX  IS:'); 
CALL  BPRINT  (XBAT) ; 

CALI  SOLVE  (XHAT,HSC1); 


EXCHANGE:  PROCEDORE  (MATRIX, PIVROW , IIVCOL) ; 
DECLARE  MATRIX  (*,♦)  FLCAT  BIN, 
(PIVFON,PIVCOL)  FIXED  BIN, 

(PIVOT, TEMP)  FLOAT  BIN, 

(ROWS, COLS, I, J)  FIXED  BIN; 

DECLARE  TEBPA  CHAR  (1); 

/*  DETERMINE  NOBBER  OF  RCNS  AND  COLOBNS  •/ 
RCNS  = HBCOND  (BATRIX.I); 

COLS  = HBCOND  (MATRIX, 2); 
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/*  BEPEH  TO  THE  PIVOT  ILEHENT  AS  ‘PIVOT*  */ 

PIVOT  = HATBIX  (PIVBOH,PIVCOL) ; 

/♦  APPLl  THE  'HECTASGIE  PULE* 

DO  I = 1 TO  PIVBOH-1,  PIVBOH+1  TO  BOWS; 

TEMP  = HATHIX  (I,PIVCOL); 

DO  J = 1 TO  PIVCOL-1,  PIVCCL*-1  TC  COLS; 

HATRIX  (I,J)  = BATEIX  (I,J)  - 
TEMP  * HATPIX  (PIVROW,J)  / PIVOT; 

END; 

END; 

/*  CHANGE  THE  OLD  PIVOT  POW  */ 

DO  J = 1 TO  PIVCOL-1,  PIVCOH-1  TC  COLS; 

HATBIX  (PIVHOH,J)  = -BATRIX  |FIVROW,J)  / PIVOT; 
END; 

/*  CHANGE  THE  OLD  PIVOT  COLDBN  ♦/ 

DO  I = 1 TO  PIVBOH-1,  PIVROH+1  TO  ROiS; 

BATRIX  (I,PIVCOL)  = BATRIX  (I,PIVCOL)  / PIVOT; 

END; 

/♦  CHANGE  THE  PIVOT  ♦/ 

BATRIX  (PIVROH,PIVCOL)  = 1 / PIVOT; 

/-»  ALTER  THE  »TOP’  AND  ‘SIDE*  AEFAYS  APPROPRIATELY  */ 
TEBPA  = SIDE  (PIVROH) ; 

SIDE  (PIVROH)  = TOP  (PIVCCL) ; 

TOP  (PIVCOL)  = TEBPA; 

END  EXCHANGE; 


BPRINT:  PROCEDURE  (BATRIX) ; 

DCL  BATRIX  (♦,*)  FLOAT  BIN, 

(BOHS,COLS,I,J)  FIXED  BIN;  ' 

/*  DETERBIME  NDBBER  OF  ROWS  AND  COLOBNS  */ 

ROBS  = HBOUHD  (BATRIX, 1); 

COLS  = HBOUHD  (BATRIX, 2); 

/*  PRINT  THE  *TOP*  ARRAY  ♦/ 

PUT  SKIP  EDIT  (*  *)  (X(4),A); 

DO  1=1  TO  COLS  -1; 

PUT  EDIT  (TOP  (I),*  •)  (X  (8)  , A,X  (6)  ,A)  ; 

END; 

PUT  SKIP; 

/*  PRINT  THE  BATRIX  LINE  BY  LINE  ♦/ 

DO  1=1  TO  ROHS; 

/*  EACH  BOH  OF  THE  HATHIX  STARTS 
ON  A NEH  OUTPUT  LINE  */ 

PUT  SKIP  EDIT  (SIDE  (I))  (X  (4)  ,A)  ; 

DO  J=1  TC  CCLS; 

/*  OBIT  THE  HOHD  ‘SKIP*  - SC  HE  HILL 
PRINT  THE  ENTIRE  BOH  ON  A NEH  OUTPUT  LINE  */ 
POT  EDIT  (HATHIX  (I  ,J)  ) (F  ( 16,  6)  ) ; 

END; 

END; 

POT  SKIP; 
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1 


END  J1PNINT; 


SOLVE:  PPOC  (HAT, SOL)  ; 

DCL  HAT(*,^)  FLOAT  BIN, 

SOL  (*)  FLOAT  BIN, 

(SUH,ANS)  FLOAT  BIN, 

(R,C,MIN,K,KK,P,Q)  FIXED  EIN; 

POT  SKIP  LIST  (’ENTEB  THE  SOLVE  POOTINE*); 

R = HBODND  (HAT,  1)  ; 

C=HBOOND  (HAT, 2)  ; 

HIN  = C - 1; 

IF  R < HIN  THEN  HIN=E; 

BEGIN; 

DCL  (TEHP(R,C))  FLOAT  EIN; 

DCL  TEHPTOP  (N-1)  CHAR(1); 

TEHF=HAT; 

TEHFTOF=TOP; 

DO  K=1  TO  HIN; 

CALL  CHOOSE  (HAT,P,Q) ; 

IF  P=0  THEN  DC; 

PDT  SKIP  LIST 

('SOLVE  COHLD  NOT  DO  ALL  INTENDED  EXCHANGES'); 
GO  TO  NCHCRE; 

END; 

CALI  EXCHANGE  (HAT,P,Q); 

END; 


NOHORE: 

PDT  SKIP  EDIT  {'THE  SCLDTICN  FOUND  IS:  ' ) (A)  ; 
SOL=0; 

DO  K=1  TO  N-1; 

DO  KK=1  TO  H; 

IF  TEHPTOP{K)  = SIDE(KK)  THEN  DO; 

SOL  (K)  = HAT (KK,C) ; 

KK  = H; 

END; 

END; 

PUT  EDIT  (TEHPTOP  (K)  , ' = *,SOL(K)) 

(X(5)  ,A,A,F{15,6)) 

END; 


PDT  SKIP(2)  LISTCVERIFY  OOP  RESOLTS'); 

PDT  SKIP  EDIT  ('THERE  ARE  ',R,'  EQOATICNS  TO  SATISFY') 
(A,F  (4)  ,A)  ; 

PUT  SKIP; 
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DO  K=^  TO  R; 

SOM  = 0; 

DO  FK=1  TC  N-1; 

SDa=SOM*TEHP  {K,KK) *SOL  (KK)  ; 

END; 

ANS  = -TEnP{K,C)  ; 

POT  SKIP  EDIT  (SUM,*  SHOULD  = •,ANS) 
(P(16,6)  ,A,F(16,6))  ; 

END; 

END;  /*  OE  BEGIN  BLOCK  IN  SOLVE*/ 

END  SOLVE; 


CHOOSE:  FROC  (HAT,PH,PC)  ; 

DCL  HAT(*,*)  FLOAT  BIN, 

(PR, PC)  FIXED  BIN, 

MAX  FLOAT  BIN, 

(I,J)  FIXED  BIN; 

MAX  = IE-6; 

PR  = 0; 

DO  1=1  TO  H; 

IF  SIDE(l)  < THEN 

/♦THIS  IS  AN  ELGIBLE  HON  •/ 

IF  ABS  (HAT  (I,J) ) > MAX  THEN 

/*  THIS  IS  THE  BEST  PIVOT  SC  FAR*/ 
DO; 

HAX=AES(HAT(I,J)) ; 

PR=I; 

PC=J; 


END; 

END  CHOOSE; 


END; 


POT  SKIP  DATA  (MSOL) ; 


END  DTNAMIX; 

/*  END  OF  PROC  M EA,  N OKN  */ 

/**************«******  ♦**♦****«  ****«««**«  ***/ 


END  VISSOL; 
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APPENDIX  B:  EVENT  DRIVEN  SINDLAIION  OSING  GPSS 


SIMULATE 

* 

* INITIALIZATION  CF  MODEL  PARAHETEFS 

* 


TIMER  CONTROLS  AND  RN  SEEDS 
5MOLT  600,6CC,600 

INITIAL  XSTIMEF,2160C000/X$TINCR,72C00C0 


PROCESSOR  MEAN  SERVICE  TIMES 
INITIAL  XH$MEIP,600 

INITIAL  XHSMCCMP, 1375 

INITIAL  XH$HCEM,U125 

INITIAL  XHSMPEP,660C 

INITIAL  XHJMTT1, 15000 

INITIAL  XHSMTT2,150  00/XH?MTT3,  1 5000/XH f HTT4 , 15 000 

INITIAL  XHSNRT1,6/XH$NBT2,2/XHSNRT3,  l/XHSNRTU,  1 


INITIAL 


XHSNOSER, 10 


RTCNT 

* 

MATRIX 

MX, 2, 4 

* 

ENT 

MATRIX 

MX, 5, 5 

FUNCTION 

DEFINITIONS 

EXPD 

FUNCTION 

RN1 ,C24 

c,o/.i. 

. 104/. 2,. 

222/. 3, .355/ 

* 


MAFTIX  EOF  ROUTE  COUNTS 
ROB  1 - COUNT 

ROW  2 - NORMALIZED  AT  1000 


EIP  SERVICE  TIME  DIST 
.509/. 5,.  6 9/. 6,, 9 15/. 7, 1.2 


BTT 

EDBH 

BPEP 

BPIP 


BRANCH  CONDITIONS  CN  BOOLEAN  VARIABLES 


EVARIABLE 

EVARIABLS 

EVARIABLE 

EVARIABLE 


P3'GE*1+P2*E' 1 
(P2*E»2+P2»E' 3) * (P3' E'C) 
(P2*E*3*P3'E» 1) ♦ (P2»  E»  4»P3' E»0) 
(P2' NE'3*F3 'E* 1) ♦ (P2 • E • 3*P3 » E*  2) 


• VARIABLES 

7INIT  FVA’»IABLE 
VRUN  FVARIABLE 
SDMI  VARIABLE 
SnM2  VARIABLE 
SUM  VARIABLE 
PERI  FVARIABLE 
PER2  FVARIABLE 
PER3  FVARIABLE 


XJTIMER/36000CO 
X$TINCR/36000C0 
MX5PTCNT(1,  1)  -t-MIERTCNT  (1,2) 
MX$HTCNT( 1,3) ♦MXIRTCNT  (1,4) 
VSSDM1»V$SDM2 
(MXIRTCNT  (1,  1)/VSSDM)  *1000 
(MXIRTCNT  (1 ,2) /VISUM) *1000 
(MXIRTCNT  (1,3)/VSSUM) *1000 


PI;  BESPOHSl  TIHE 

P2;  FOOTE  NUnEEP  (1  ,2,3,4) 

P3;  IN  - 0 

COT  - 1 POP  RCOTES  2 AND  4 

3 PCP  ROUTE  3 
P4;  KEAN  THINK  TIKE 


ENR1  GENEPATl  , , , XH SNRT 1 , , 6 , ? 
ASSIGN  2,1 

TRANSFER  ,TERH 


ENR2  GENEPATl  , , , XH JNPT 2 , , 6 , F 

ASSIGN  2,2 

TRANSFER  ,TEHH 
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* 


ENR3  GENERATE 
ASSIGN 
TRANSFER 

* 


, , , XH$NRT3,  ,6  ,F 

2,3 

,TERR 


ENR4  GENERATE  , , , X H SNPT4 , , 6 , F 
ASSIGN  2,4 


* 

* 

♦ 

* RCDEL  3EGHENT  FCR  THE  COBPCSITE  NETWORK 


TERR 

TEST  E 

P2, 1 , OTHER 

ASSIGN 

4,XHINTT1 

TRANSFER 

, ENTER 

OTHER 

ASSIGN 

4,XH JNTT2 

ENTER 

QUEUE 

TERBC 

TERM  STAGE 

eSAVEVALUE 

RTCNT-»,  1,P2 

, 1,NX 

HSA7EVALDE 

ENT*  , 1 ,P2, 1 

,MX 

ENTER 

TERMS 

ADVANCE 

P4 

LEAVE 

TERMS 

DEPART 

TEPMQ 

ASSIGN 

3,^ 

DESIGNATE  AS 

NARK 

1 

EST  INTO  THE 

* 

* 


FIP  PROC 

ESSOF 

FIP 

QUEUE 

FIPQ  FTP  STAGE 

MSAVEVALUE  ENT-*- , 2 , P2  , 1 , MX 

SEIZE 

FIPF 

ADVANCE 

XHSHFIP,FN$EXFD 

RELEASE 

FIPF 

DEPART 

FIPQ 

TEST  E 

EVS3TT, 1,CCHF 

TEST  E 

P2, 1 ,EFT1 

ROD"'E  RE 

SPONSE  TIME  TABDIATTONS 

TABULATE 

RESR  1 

T'’ANSFER 

,TERM 

INWARD  TO 
NETWORK 


THE  NETWORK 
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BPT1 

* 

BPT2 

* 

BPT3 

* 

BPT4 

* 

* 

* 

* 

C01P 


* 

* 

* 

DBH 


♦ 

* 

« 

PEP 


TEST  E 

TABOtATE 

TRANSFER 

TEST  E 

TABULATE 

TRANSFER 

TEST  E 

TABULATE 

TRANSFER 

TRANSFER 


P2,2,EFT2 

FESR2 

,TEBS 

P2,3,EPT3 

HESR3 

,TEFH 

P2,4,BPT4 

BESR4 

,TERf! 

,TEEP 


CCaP  PROCESSOR 

QUEUE  conpc 

BSAVEVALUE  ENT+ , 3 , P2 , 1 , H X 

SEIZE  COflPF 

ADVANCE  XH$HCCBE,FNJEXFD 

RELEASE  COBPF 

DEPART  COflPC 

TEST  E EV$BtBR,0 ,DBa 

TEST  E BVSBPEP,0,PEP 

TEST  E BVSBFIP,0,FIP 


EBB  PROCESSOR 

QUEUE  DBBQ 

BSAVEVALUE  ENT*  , 4 , P2 , 1 , BX 

SEIZE  DBBF 

ADVANCE  XH$BDBB,PNSEXFD 

RELEASE  DBBF 

DEPART  DBBQ 

ASSIGN  3*,1  INCREBENT  tIRECTICN  COUNT 

TRANSFER  ,COBF 


PEP  PROCESSOR 


QUEUE 


PEPQ 


BSAVEVALUE  ENT* , 5 , P2 , 1 , BX 
SEIZE  PEPF 

ADVANCE  XHSBPEP,PN$EXPD 
RELEASE  PEPF 

DEPART  PEPQ 

ASSIGN  3*,1  INCFEBENT  DIRECTION  COUNT 

TRANSFER  ,COBE 


I 

I 

I 


1 
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CNTRL 


CUT  PL 


CPTRL 


CNTPL 


* 

* 

* 

♦ 

CNTPL 


MODPL  SEQUENT  ECP  TIHEP  PUN  CONTPCL 


GENERATE 

nSAVEVALOE 

fISAVEVALOE 

nSAVEVALOE 

NSAVEVALUE 

SAVEVALUE 

ESAVEVALUS 

nSAVEVALOE 

nSAVEVALOE 

nSAVEVALOE 

■SAVEVALUE 

SAVEVALOE 

SAVEVALUE 

TEPniNATE 


XSTinER 

PTCNT,2,1,VSFEP1,nX 
HTCNT,2,2,7$PEP2,nX 
RTCNT,2,3,V$PEB3,HX 
RTCNT,2,4 ,VSPEFa,nX 
CREST, VSnPEST, XI 
ENT,  1, 5,VSTCST,nX 
ENT,2,5,V$PCNT,nX 
E»IT,3,5,V$CCNT,nX 
ENT,4,5,V$DCNT,nX 
ENT,5,5,7SPCNT,nX 
TINIT,V$VINIT,X1 
TPON,V$VRnN,Xl 
1 


NCPnALIZED  VALOES  IN  ROW  2 


CONTROL 

NCXPEP 

LOAD 

START 

INITIAL 

INITIAL 

RESET 

GENERATE 

START 

INITIAL 

INITIAL 

INITIAL 

CLEAR 

STORAGE 

GENERATE 

START 

INITIAL 

INITIAL 

RESET 

GENERATE 

START 

INITIAL 

INITIAL 

INITIAL 

CLEAR 

STORAGE 

GENERATE 

START 

INITIAL 

INITIAL 


CARDS 


DAG06,DAG07 

1 

nxi  (1-2, 1-4)  ,0 
nX2(1-5,1-5) ,C 

XTTINCR 

1 

PECONFIGOEE  FOR  20  JOBS 

XHSNDSER,20 

XHSNRTI  ,12/XHSNRT2,4/XH$NRT3,2/XH$NRT4,2 
XSTinER,432C0000/X$TINCB, 1440CC00 
XH1-XH13, X1-X2 
SSTEPnS,20 
XSTinPR 
1 

nxi  (1-2, 1-4)  ,0 
nx2  (1-5, 1-5)  ,0 

XSTINCR 

1 

RECONFIGOEE  FOR  30  JOBS 

XHSNOSEP, 30 

XHSNRT1,18/XH$NET2,6/XHINRT3,3/XH$NET4,3 

XSTIHER,64800000/X$TINCE,21600COO 

XH1-XH13,X1-X2 

SSTERnS,30 

XSTIREP 

1 

nxi  (1 
nx2  (1- 
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PESET 

CNTRL  GENEFATE 


START 

REPORT 

EJECT 

TEXT 

TEXT 

TEXT 

TEXT 

18 

TEXT 

18 

TEXT 

18 

TEXT 

TEXT 

26 

TEXT 

26 

TEXT 

26 

TEXT 

26 

TEXT 

SPACE 

FBS 

TITLE 

SPACE 

TEXT 

TAB 

INCLUDE 

SPACE 

QUE 

TITLE 

QDZ 

INCLUDE 

SPACE 

FAC 

TITLE 

FAC 

INCLUDE 

SPACE 

FBS 

TITLE 

EJECT 

XSTINCF 

1 


NOHBER  CF  JOBS  IS  tX BSSDSER , 2/XXX # 
IHITIAIIZATICN  HCOPS  IRE  #XL$TINIT, 2/XXX . XX# 
DATA  COLLECTICN  BOORS  ARE  iXLJTRON, 2/XXX. IX# 
THINK  TIFIES  ARE  - A;  #XH$BTT1, 2/XXXXXX  # 

B:  #XHSBTT2, 2/XXXXXX# 

C:  #XH*nTT3,2/XXXXXXI 
D:  #XHSBTT4, 2/XXXXXX# 

SERVICE  TIHES(HSEC)  APE  - 
PIP:  EXP  fXHSHFIP, 2/XXXXXX# 

COBF;  EXP  #XH$RCOBP, 2/XXXXXX# 

DBB:  EXP  #XHSHtBB, 2/XXXXXX# 

PEP:  EXP  #XHSHPEP, 2/XXXXXX# 

1 

RTCNT,THE  JOB  CCONT  BT  ROUTE 
3 

COBPCSITE  BEAN  F TIBE  IS#XL$CR EST, 2/IX XXX . XXX  I 
T1-T4/1,2,3,«» 

3 

,QOEOE  STATISTICS 

Q1-Q5/1,3,7 

3 

.PROCESSOR  STATISTICS 

F1-P4/1,2,3,4 

3 

ENT, PROCESSOR  ENT  COUNT  BT  RCOTE 


* FOLLCNING  DENSITY  PLOTS  ARE  OPTIONAL 

* 


TAB 


30 


TAB 


30 


TITLE 

EJECT 

GRAPH 

ORIGIN 

X 

T 

STATEBEHT 

ENDGRAPH 

EJECT 

TITLE 

EJECT 

GRAPH 

ORIGIN 

X 

T 

STATEBENT 

ENDGRAPH 

EJECT 


RSSR1, ROUTE  1 RESPONSE  TIBE  FREQUENCY 

TF.RESRI 
53,  1C 

,1,0, 100, 1,100, NO 
0,10,50,1 

1,31, ROUTE  1 RESPONSE  TIBE  FREQUENCY 


RESR2, ROUTE  2 RESPONSE  TIBE  FREQUENCY 

TF, RESR2 
53, 1C 

,1,0,2000,1 , 100, NO 
0»1,50,1 

1,31, ROUTE  2 RESPONSE  TIBE  FREQUENCY 


1 
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TAB 

TITLE 

EJECT 

PESR3,FCnTE  3 RESPONSE  TIME  FEZQOENCY 

GPAPH 

TP,RESR3 

CP.IGIN 

53,  1C 

X 

,1,0,2000,1 ,100, NO 

Y 

0,1,50,1 

30 

STATEMENT 

ENDGEAPH 

EJECT 

1,31, FOOTE  3 FESPONSE  TIME  FFECOENCT 

TAB 

TITLE 

EJECT 

FESRU, FOOTE  U RESPONSE  TIME  FFECOENCT 

GRAPH 

TF,aESF4 

CPIGIN 

53,10 

X 

,1,0, 2000, 1,100, NO 

Y 

0,1,50,1 

30 

STATEMENT 

ENDGEAPH 

EJECT 

END 

1,31,POO'^E  4 FESPONSE  TIME  FFECOENCT 

I 


1 
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APPENDIX  C:  G COnEOTATION  PFCGHAH  FCE  COBPOTER  NETWORKS 
WITH  HDLTIPIE  JCB  TYPES 


NORH:  PRCC  OPTIONS  (BAIN); 

/♦  COBPDTES  THE  NOPB ALIZATICN  CONST  ;STS  FOP  H NODES  AND 
DP  TO  FGDP.  JOE  CLASSES*. 

A - 1ST  JOB  CLASS 
NA  - NDBBER  OF  A JOES 
B - 2ND  JCB  CLASS 
NB  - NOBEER  OF  E JOBS 
C - 3RD  JOB  CLASS 
NC  - NOBEER  OF  C JOBS 
D - 4TH  JOB  CLASS 
ND  - NDBBER  OF  t JOBS 

G(NA,NB,NC,ND)  IS  AN  ARRAY  FOR  THE  LAST  OP  HIGHEST  VALOE  OF 
B AND  IS  BAINTAINED  IN  STORAGE.  THE  B IS  IBPLIED. 

XA  (B)  ,XB  (B)  ,XC  {B)  ,XD  (H)  APE  THE  RESPECTIVE,  RELATIVE 
OTILIZATIONS  FCR  JOB  CLASSES  A,E,C,D.  */ 

/*  INPDT  IS  FREEFCRBAT: 

NO  OF  CLASSES,  NCLASS; 

B,  THE  NDBBER  CF  NODES; 

HO  CF  IHDIVIDOAL  CLASSES,  NA,NE,NC,ND; 

RELATIVE  OTILIZATIONS,  XA,XB,XC,XD. 

THESE  ARE  VECTORS  CF  LENGTH  B; 

SERVICE  HATES  AT  THE  NODES  (1,2, 3, 4); 

A FLAG  FOR  PRINTER  OOTPOT  OF  THE  G,  -Y  CP  N; 


A OOTPOT  FILE,  G400T,  BOST  BE  ALLOCATED  FOP  THE 
G BATRIX.  V 

DCL  (LA,LB, LC,1D)  FIXED  BIN; 

DCL  (B, NCLASS, NA,NB,NC,ND)  FIXED  BIN  ; 

DCL  (IA,IB,IC,ID)  FIXED  BIN; 

GET  LIST  (NCLASS, B)  ; 

POT  SKIP  DATA  (NCLASS, B)  ; 


/♦ 

DETERBIHE 

THE 

NOBEER 

OF  CLASSES 

IP 

NCLASS=1 

THEN 

GC 

TO 

FINISH; 

IF 

NCLASS=2 

THEN 

GC 

TO 

FINISH; 

IP 

NCLASS=3 

THEN 

GC 

TO 

FINISH; 

IF 

NCLASS=4 

THEN 

GC 

TO 

CLASS4; 

Appendix  C 


16C 


/*  ESTAELISHED  FOR  « CLASSES. 
OTHERS  ARE  SOECASES  CF  4 AND  ARE 
NOT  INCLUDEED  IN  THIS  LISTING.  */ 


CLASS4:  /*  CASE  CF  FOUR  JCE  CLASSES  */ 

GET  LIST  (NA,NB,NC,ND)  ; 

POT  SKIP  DATA  (NA,NE,NC,ND)  ; 

BEGIN; 

DCL  (XA  (H)  ,X3  (H)  ,XC  (M)  ,XD  (H)  ) FLOAT  EIN; 
DCL(G(0:NA,0;NE,0:NC,0:ND)  , GC  (0 ; NC,  0 : H)  , GD  (0:  ND,0:  M)  , 
GA  (0:NA,0:H)  , GB  (C:NB,0:B) ) FLOAT  EIN; 

DCL  (GAB  (0:NA,0:NB,0:n)  ,GAC  (0:NA,0;NC,0:H)  , 

GBC  (0:NB,0:NC,0:K)  ) FLOAT  BIN; 

DCL  (GAD (0: NA,0:ND,0:n)  ,GBD  (0  :NB,0:Nr,0:a) , 

GCD (0;NC,0;ND,0:n)  ) FLOAT  BIN; 

DCL  (GABC (0:NA,0: NE,0;NC,0: H)  , 

GABD  (0:NA,0:NE,0:ND,0:«)  , 

GACD(0:NA,0:NC,0:ND,0:B)  , 

GBCD (0:NB,0: NC,0:ND,0: F)  ) FLOAT  EIN; 

DCL  (0  (PI)  ,0A  (M)  ,OB  («)  ,OC{B)  ,DD  («)  ) FLOAT  EIN; 

DCL  (T  {?!)  ,TA  (B)  ,TB  (B)  ,TC(H)  ,TD  (f!)  ) FLOAT  EIN; 

DCL  BO(H)  FLOAT  BIN; 

DCL  (TSA,TSB,TSC,TSD)  FLOAT  BIN; 

DCL  GPRNT  CHAR  (1)  ; 

GEDIT:  FOEBAT  (A,F  (NA/IO+I)  , A,F  (NE/10  + 1)  ,A,F  (NC/10+ 1) , 

A,F(ND/10  + 1)  ,A,E(15,8)  ,SKIP)  ; 

TEDIT;  FORBAT  (A,SKIP,  (3)  (A,  F (2)  ,X  (3)  ) ) ; 

GET  LIST  (XA,XE,XC,XE)  ; 

POT  SKIP  DATA  (XA,XE,XC,XD)  ; 

GET  LIST  (BO)  ; 

POT  SKIP  DATA  (BO)  ; 

GET  LIST  (GPRNT)  ; 

IF  NA=0  THEN  POT  SKIP  LIST  ('CLASS  EEECR'); 

/*  INITIALIZE  THE  G */ 

G=0.0;  GA=0.0;  GB=0.0;  GC=0.0;  GD=0.C; 

GAB=0.0;  GAC=0.0;GEC=0. 0; 

GAD=C.C;  GBD=0.0;  GCD=0.0; 

GA3C=0.0;  GABD=0. 0 ;G ACD=0 . 0 ; GECD=0 . 0 ; 

CALL  ONEJOB (NA,XA,GA) ; 

CALL  ONEJOB  (NE,XB, GB) ; 

CALL  CNEJOB  (NC,XC,GC) ; 

CALL  ONEJOB  (ND,XD,GD)  ; 
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CALL  TWOJOB (NA,NB,XA,XE,GAB,GE,GA)  ; 
CALL  TWOJOB  (N A,«IC,XA,XC,GAC,GC,GA)  ; 
CALL  TB0.70B(NE,!JC,XE,XC,GBC,GC,GB)  ; 
CALL  TUOJOB  (N A,ND,XA,XD,GAD,GE,GA)  ; 
CALL  TSOJOB (NB,BD,XB,XD,GEC,GD,GB)  ; 
CALL  TBOJOB  (NC,ND,XC,XD,GCD,GD,GC)  ; 


CALL  THREJOB(NA,NB,HC,XA,XB,XC,GABC,GEC,GAC,GAB)  ; 

CALL  THeEJOB(NA,NB,ND,XA,XB,XC,GABB,GBD,GAD,GAB)  ; 

CALL  THREJOE(NA,NC,IIE,XA,XC,XD,GACr,GCD,GAD,GAC)  ; 

CALL  THEEJOE  {NE,NC,BD,XB,XC,XD,GBCD,GCC,GED, GEC)  ; 

/*  CCHPOTE  THE  G */ 

CALI  PODRJOB  (NA,NB,RC,ND,XA,XE,XC,XE, 

G,GBCD,GACD,GAEE,GABC) ; 

NT  = NA<-NE+NC  + NB; 

IF  GPHUT=*1*  TREK  CALI  PDTG  (R A , NB , NC , ND , G) ; 

ELSE  DO; 

PET  FILE{Gt»OOT)  EDIT  (NA,  NE,NC  , RE)  ( (4)  (F  (2)  , X (2)  ) , SKIP)  ; 
POT  PILE(G400T)  EDIT  (G)  ( (4)  (E  ( 1 5 ,8)  , X (3)  ) , SKIP)  ; 

END; 

POT  SKIP(2)  EDIT(*G  (»,NA-1, *, • ,NB, • , '^NC, • , • ,ND, •)=  », 

G (NA-1,NB,NC,ND))  (B(GEEIT)); 

POT  SKIP  (2)  EDIT{'G{*,HA,*,’»NE-1,',*,NC,*,*  ,ND,')= 

G (NA,NB-1,NC,ND))  (R  (GEEIT))  ; 

POT  SKIP(2)  EDIT  ('G  f*  ,HA,  ’ , SNE,  ’ , * ,NC-1,’ ,•  ,ND,  •)  = '# 
G(NA,NE,NC-1,HE))  (B  (GEEIT)) ; 

POT  SKIP  (2)  EDIT('G(*,NA,*,',NE,',',NC,*,',ND-1,*)= 

G {NA,NB,NC,ND-1) ) (F (GEDIT) ) ; 

POT  SKIP(2)  EDIT('G{*,NA,',',NB,*,',NC,*,',NB,')=  *, 

G (NA,NE,NC,ND) ) (B  (GEDIT) ) ; 

/*  CORPDTE  THE  OTIIIZATION  AND  THROPDT  FACTORS  */ 

DO  1=1  TO  B; 

DA  (I)  =XA  (I)  *G  (NA-1  ,HB,NC,HD)/G  (»  A , NB  , HC  , ND)  ; 

TA(I)=OA  (1)*B0  (I)  ; 

OB  (I)  =XB  (I)  *G  (NA,NB-1,NC,ND)  /G  (BA,NB,NC  ,HD)  ; 

TB  (1)  =0B  (I)*BO  (I)  ; 

OC(I)=XC  (I)  *G  (»A,NB,NC-1,ND)/G  (NA,NB,NC  ,ND)  ; 

TC(l)=OC  (I)*BO  (I)  ; 

OD(I)=XD  (I)  *G  (NA,NB,BC,HD-1)/G(KA,NB,NC,ND)  ; 

TD  (I)  =OD  |I)*BD  (1)  ; 

0(1)=  DA  (I)-»OB  (I)+OC  (I)  *00(1)  ; 

T(I)=TA(I)+TB  (I)+TC(I)  ♦TD(I)  ; 

POT  SKIP  (2)  EDITC'NODE  *,1,*;  0=  *,0(1), 

• 0A=  •,0A(I),»  0B=  ',08  (I),*  0C=  •,OC(I),*  0D= 
OD(r),»  T*  »,T{I),'  TA=  *,TA(I),'  TB=  '^TBd), 

• TC*  TD=  »,TD(I)) 

(A,P(2)  , A,  SKIP,  (2)  ((5)  (X(3),A,F(5,3))  ,SKIP))  ; 


END 
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/*  CCnFDTE  BRANCHING  RATICS  */ 

PA  = TA  {1)/T  (1)  ; 

PB=T3{1)  /T  (1)  ; 

PC  = TC  {1)/T  (1)  ; 

PD='ID  (1)/T  (1)  ; 

PUT  SKIP  (3)  EDIT { 'THRDPOT  BRANCHING  IS:', 

• PA=  ',PA,'  EE=  •,PB,'  EC=  *,PC,'  PD=  •,PD) 
(A, SKIP,  (4)  <X  (3)  ,A,F  (6,3)  ))  ; 

/♦  COHPDTE  SHORT  CIRCUIT  THBUPUT  ♦/ 

TSA=G  (N A-1,NB,NC, NE) /G (NA, NE, NC,ND)  ; 

TSB=G  (NA,NB-1,NC, NE) /G (NA , NE, NC,ND)  ; 
TSC=G(NA,NB,NC-1,NE) /G (NA, NB,  NC,ND)  ; 
TSD=G(NA,NB,NC,ND-1)/G (NA , NE, NC,ND)  ; 

PUT  SKIP  (3)  EDITCSHORT  CIRCUIT  THRUFUT  FOLLOWS:  • 

• TSA=  »,TSA,'  TSB=  •,TSE,'  TSC=  * ,TSC, 

• TSD=  *,TSD) 

(A, SKIP,  (4)  (X  (3)  ,A,F(6,3)))  ; 


/♦COMPUTE  THE  THRUPUT  AS  A FUNCTION  CF  THE 
NUHBER  OF  JOBS  IN  EACH  JCE  CLASS  WITH  THE 
OTHER  JOES  HELD  FIXED.*/ 

PUT  SKIP  (3)  EDITCTHRUPDT  VS.  AJOBS', 

♦NB=  •,NB,'NC=  •,NC,'ND=  ' ,ND)  (R  (TEDIT) ) ; 
DO  IA=1  TO  NA; 

TSA=G  (IA-1  ,NB,NC,ND) /G  (IA,NB, NC ,ND)  ; 

OUT  SKIP  EDIT('NA=  •,IA,'TSA=  ' ,TSA) 
(A,F(2),X(3),A,F(7,3)); 

END; 

POT  SKIP(3)  EDIT ('THRUPUT  VS.  EJOBS*, 

'NA=  ',NA,'NC=  ',NC,'ND=  ', ND)  (E  (TEDIT)  ) ; 
DO  IB=1  TO  NB; 

TSB=G  (NA,IB-1 , NC, ND) /G (N A , IB , NC , NO)  ; 

POT  SKIP  EDIT('NB=  ',IE,'TSB=  ',TSB) 
(A,F(2)  ,X(3),A,P(7,3))  ; 

END; 

POT  SKIP  (3)  EDIT ('THRUPUT  VS.  CJOBS', 

'NA=  ',NA,'NB=  ',NB,'RD=  ',ND)  (B  (TEDIT)  ) ; 
DO  IC=1  TO  NC; 

TSC=G(NA,NE,IC-1,ND)/G  (NA,NB,  IC,ND)  ; 

PUT  SKIP  EDIT(*NC=  •,IC,»TSC=  ' ,TSC) 
(A,F(2)  ,X(3)  ,A,F(7,3))  ; 

END; 

PUT  SKIP  (3)  EDITCTHRUPUT  VS.  EJOBS', 

'NA=  ',NA,'NB=  ',NB,«NC=  ', NC)  (R  (TEDIT) ) ; 
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DO  ID=1  TO  ND; 

TSD=G  (NA,NB,NC,  ID-1) /G  (S A , NB,  NC  , ID)  ; 

POT  SKIP  EDIT(«ND=  *,IE,*TSD=  ’ ,TSD) 

(A,F(2)  ,X  (3),  A,P(7,3))  ; 

END; 

END;  /*  BEGIN  BLOCK  C1ASS3  •/ 

GO  TO  FINISH; 

ONSJOB:  PB0C(N1,X1,G1) ; 

DCL  N1  FIXED  BIN; 

DCL  XI  (*)  FLOAT  BIN; 

DCL  G1  {*,*)  FLCAT  EIN; 

G1  (C,*)  = 1.0; 

DO  Ln=1  TO  M; 

DO  LA=1  TO  N1  ; 

G1  (LA,LH)=G1  (IA,LR-1)  ♦ X1(Ln)  ♦ G 1 {LA- 1 ,LH)  ; 

END; 

END; 

END  CNEJCB; 


TWOJOB:  PROC(N1,N2,X1,X2,G2,G221,G2Z2) ; 

DCL  (N1,N2)  FIXED  EIN; 

DCL  (XI  (♦)  ,X2  (•) ) FLCAT  BIN; 

DCL  (G2  (*,♦,♦)  ,G2Z1  (*,*)  ,G2Z2  (*,*) ) FLOAT  BIN; 

DO  Ltl=1  TO  H; 

G2 (*,0,Ln)=G2Z2  {♦,!«)  ; 

G2 (0,*,LH)=G2Z1  (*,LB)  ; 

DC  LB=1  TC  N2; 

DO  LA=1  TC  N1; 

G2(IA,LB,Ln)  = G2  (LA,lE,in-1)  ♦ 

XI  (LN)  *Q2  (LA-1, LB, IB) 

♦ X2  (LB)  *02  (LA, LB-1, IB)  ; 

END; 

END; 

END; 

END  TWOJOB; 

THREJOB;  PRCC  (N1 , N2 ,N3 ,X1 , X2 , X3,G 3 , G3Z 1 , G3Z 2 ,G3Z3) ; 

DCL  (N1,N2,N3)  FIXED  BIN; 

DCL  {XI  (•)  ,X2(*),X3(*))  FLCAT  EIN; 

DCL  (G3(*,*,*,*),G3Z1(*,*,*), 

G322(*,*,*)  ,G323  {*,♦  *' ) FLOAT  EIN; 

DO  LB=1  TO  B; 

G3(0,*,*,LB)=G3Z1(*,*,IR)  ; 

G3  (*,0,*,  LB)  =G3Z2  (•,*,!>!)  ; 
G3(*,*,0,LH)=G3Z3(*,*,1H)  ; 


Appendix  C 


16U 


DO  LC=1  TC  N3; 

DO  LB=1  TC  N2; 

DO  1A=1  TO  N1; 

G3 (LA,L3,LC,1H)=G3 (LA,LB,LC,LH-1) 

♦ XI  (IS)  *G3  (LA-1  ,LE,1C,LI1) 

♦ X2  (L«)  *G3  (LA,LB-  1,1C,L?)) 

♦ X3  (L?1)  *G3  (LA,LB,  IC-1,LB)  ; 

END; 

END; 

END; 

END; 

END  THBEJOE; 

FODPJQB:  PHOC (N1 ,N2,N3,N4,X1 ,X2,X3,X4, 

Gt:,G4Z  1,G4Z2,G4Z3,G4Z4)  ; 

DCL  (N1,N2,N3,N4)  FIXED  BIN; 

DCL  (XI  (*) ,X2  {*) ,X3  ♦),X4(*))  FLOAT  BIN; 

DCL  (G4 (*,*,*,*)  ,G4Z1  (*,♦,*,•)  ,G4Z2  (♦,*,*,*), 

G423 (*,*,*,*)  ,G4Z4  (*,*,*,*) ) FLOAT  BIN; 

DO  Ln=1  TO  n; 

G4  (0,*,*,*)=G471  (*,*,*,!(!)  ; 

GU(*,0,*,*)=G4Z2  (*,*,*,!«)  ; 

G4  (*,*,0,*)=G4Z3{’*,*,*,L-)  ; 

G4  ( *,*,*,  0)=G4Z4  (*,*,*,!»)  ; 

DO  LD=1  TC  N4; 

DO  IC=1  TC  N3; 

DO  LB=1  TC  N2; 

DC  LA=1  TC  N1; 

G4 (LA,LE,IC,LD)= 

G4  (LA,LE,LC, IE) 

♦ XI  (LM)  *04  (LA-  1,LB,LC,  LD) 

♦ X2  (LH)  *04  (LA,LE-1,LC,LD) 

♦ X3  (Lf!)  *G4  (LA,LE,LC-1,  LD) 

♦ X4  (LS)  *G4  (L  A,LE,LC,LD-1)  ; 

END; 

END  ; 

END; 

END; 

END; 

END  FOORJOB; 


I 
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POTG:  PROC (N1 , N2,N3, Na,GF)  ; 

DCL  GP {*,*,*,*)  FLCAT  BIN; 

DCL  {N1,N2,N3,N«,NT)  FIXEC  BIN; 

POT  SKIP(3)(*G  FOLLOWS:*); 

POT  SKIP  (2)  ; 

NT=NUN2  + N3  + N4; 

LCOUNT=3; 

DO  IA=0  TO  N1; 

DO  IB=0  TO  N2; 

DO  IC=0  TO  N3; 

DO  1D=0  TO  N4; 

IF  LCCONT=0  THEN  DO; 

POT  SKIP; 

LCOONT=3;  END; 

POT  EDITC  G {•,IA,',',IB,',  ',10, ,ID, 

*)  = ’ ,GP  (IA,IB,IC,ID)  ) 

(A,F(N1/10»1)  ,A,F(N2/10+1),A,P(N3/10*1)  ,A, 
F(N4/10*1)  ,A,E  (15,7))  ; 

LCCUNT=LCOONT-  1 ; 

END; 

END; 

END; 

END; 


END  POTG; 
FINISH:  END  NOPB; 
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APPENDIX  C:  THFOPUT  COtiPUI ATIC N USING  G VALUES 


TPUT:  PRCC  OPTIONS (KAIN)  ; 

DCL  (NA,NB,NC,ND)  FIXED  BIN; 

DCL  in  (4)  FLOAT  BIN; 

DCL  NCASE  FIXED  BIN; 

DCL  (ZA,23,ZC,2D)  FICAT  BIN; 

/♦  AN  INPUT  FILE,  GVALOES,  BUST  BE  ALLOCATED.  NOTE 
THAT  IT  CONTAINS  THE  PARAKETEFS  OF  THE  NETWORK 
(SEE  FOLLOWING  STATEMENT)  AS  BELL  AS 
G COMPUTED  FROM  NOFH,  •/ 

GET  FILE  (GVALUES)  LIST (N A , NE , NC, N D , f CAS E , Z A, ZB ,ZC , ZD , HU ) ; 
PUT  SKTP(3)  EDIT('VAIUES  OF  G MATRIX  ARE;',  •NA=  ',NA, 
'NB=  •,NE,'NC=  •,NC,»ND=  *,ND) 

(A,SKIP,  (4)  (A,F(2)  ,X(3)))  ; 

PUT  SKIP  (3)  EDIT  (• NETWORK  CHARACTERISTICS  ARE;', 'CASE  ', 
NCASE, •ZA=  ',ZA,'ZB=  •,ZB,'ZC=  •,ZC,'ZD=  ’,20, 
'HD(1)=  •,HO(1)  ,'10(2)=  •,HD(2), 

•HU(3)=  • ,HU  (3)  ,'1U(4)  = *,HU(4)) 

( A.  SKIP,  A,  F (2)  ,X  (3)  , (4)  (A,F  (6 ,2)  ,X  (3))  ,SKIP, 

(4)  |A,E(15,7)  ,X  (2)  ) ) ; 

BEGIN; 

DCL  (TSA,TSB,TSC,TSD,NV,TTCT)  FLOAT  BIN; 

DCL  G(C:NA,0;NE,0:NC,0;ND)  FLOAT  BIN; 

DCL  (VA, V3,VC,VD)  FIXED  BIN; 

DCL  (LISTF,LOADP)  CHAR(1); 

DCL  RFLAG  FIXED  BIN; 

DCL  (PA, F3, PC, PD)  FLOAT  BIN; 

TEDIT;  FORMAT  ( A , X (3)  , (3)  (A  , F (2)  , X (3)  ) , SKIP)  ; 

GET  FILE  (GVALOES)  LIST(G); 

POT  SKIP{3)  LIST(’G  LISTING  DESIRED?  - I OR  N ;'); 

GET  LIST(LISTF); 

IF  LISTF='Y'  THEN  CALL  PUTG (N A , NB , NC , N D , G)  ; 

/*  COMPUTE  THE  DESIRED  THROFUTS  */ 

/•  COHPOTE  SHORT  CIRCUIT  THPDPOT  */ 

RESTART: 

PUT  SKIP  (3)  LIST(*TCTAL  JOB  VECTOR;'); 

GET  LIST  (N A,SE,NC,NE)  ; 

RFL AG=0; 

PUT  SKIP(3)  ECITCJCE  LOADING  IS;', 

'NA=  ',NA,'NB=  ',NE,'NC=  ',NC,'Nr= 

(A, SKIP,  (4)  (A,F(2)  ,X  (2)))  ; 

TSA  = G (NA-1 ,NB,NC,ND)/G (NA,NE, NC,ND)  ; 

TS3  = G (NA,NB-1 ,NC,NE) /G ( N A , N B , NC, ND)  ; 

TSC=G (NA ,NB,MC-1, NC)/G(NA,NB, NC,ND) ; 


•,ND) 
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TSD=G  (NA,NE,NC,ND- 1)/G{NA,NE,NC,ND)  ; 

TTOT='rSA+TSB*TSC  + TSD; 

POT  SKIP (3)  EDIT('SHCBT  CIPCOIT  THPDPOT  FOLLCHS;*, 

• TTOT=  •,TTOT,'  TSA=  ‘^TSA,*  TSE=  ',TSE,'  TSC=  '/TSC, 

• TSD=  *,150) 

(A, SKIP,  (5)  (X  (3)  ,A,F(6.3)))  ; 

/*  COMPUTE  ERANCHIKG  PATICS  */ 

PA=TSA/TTOT; 

PB=TSB/TTOT; 

PC=TSC/TTOT; 

PD=TSD/TTOT; 

POT  SKIP(3)  EDIT ('THPOPOT  EPASCHING  IS:', 

• PA=  ',PA,'  FE=  ',08,'  EC=  PD=  *,00) 

(A, SKIP,  (4)  (X(3),A,F(6,3)))  ; 

/♦COMPOTE  THE  THPOPOT  AS  A FONCTION  CF  THE 
NOHBEP  OF  JOBS  IN  EACH  JOE  CLASS  PITH  THE 
OTHER  JOBS  HELD  FIXED.*/ 

PFLAG=1 ; 

NEPLOAD: 

IP  RFLAG=0  THEN  DO; 

POT  SKIP  (2)  LIST 

(•INPUT  THE  ACTIVE  VALOES  CF  JOB  CLASSES;'); 

POT  SKIP; 

GET  LIST (VA,VE,VC,VD) ; 

END; 

ELSE  DO;  VA=N A ; VB= NB ; VC=NC; VD= F C ; END ; 

POT  SKIP  (3)  EDIT  ('THPOPOT  VS.  AJOBS', 

•NB=  •,VB,'NC=  •,VC,'NC=  • , VD)  (E  (TEDIT) ) ; 

DO  IA=1  TO  VA; 

TSA=G  (IA-1,NB,NC,ND) /G  (IA,NB, NC,ND)  ; 

POT  SKIP  EDIT(»NA=  •,IA,'TSA=  * ,TSA) 

(A,F(2)  ,X(3),A,F(7,3))  ; 

END; 

POT  SKIP{3)  EDIT ('THPOPOT  VS.  EJOBS', 

•NA=  •,VA,'NC=  •,VC,*ND=  * ,VD)  (F  (TEDIT)  ) ; 

DO  IB=1  TO  VB; 

TSB=G (NA,IB-1 ,NC, FD)  /G (N A ,IB , NC ,ND)  ; 

POT  SKIP  ECIT('NB=  •,IE,'TSB=  • ,TSB) 

(A,F(2)  ,X(3)  ,A,F  (7,3))  ; 

END; 

POT  SKIP  (3)  EDIT ('THPOPOT  VS.  CJOBS', 

•NA=  •,VA,'NB=  ',VB,*NC=  ', VD )( P (TEDIT)  ) ; 

DO  IC=1  TO  VC; 

TSC=G  (NA,NE,IC-1,ND)/G  (NA,NB,  IC,ND)  ; 

POT  SKIP  EDIT('NC=  ',IC,'TSC=  ' ,TSC) 
(A,F(2),X(3),A,F(7,3)); 

END ; 
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?U?  SKIP|3)  EDIT {'THPOPDT  VS.  CJOBS', 

»NA=  •,VA,»NB=  •,VB,*NC=  ' , VC ) ( E (T ED IT)  ) ; 

DO  IB=1  TO  VD; 

TSC=G  (NA.SB.NC,  ID-1)  /G  (NA,NB,NC,ID)  ; 

PBT  SKIP  EDIT{’15D=  ’,IC»’TSD=  * ,TSD) 

{A,F  (2)  (3)  ,A,P  (7,  3))  ; 

END; 

HPLAG=0 ; 

POT  PAGE; 

PUT  SKIP  (6)  LIST 

{'NEW  LOADING  DESIPED?  - Y,N  OF  P FOE  RESTART;'); 

GET  LIST  (LCADF)  ; 

IF  LOADF='T'  THEN  GO  TO  NEWICAD; 

IF  LOADF='R'  THEN  GO  TO  RESTART; 

END;  /*  BEGIN  BLOCK  V 

PDTG;  PFOC(N1,N2,N3,N4,GP) ; 

DCL  r,P  {*  ,*,*,*)  FLOAT  BIN; 

DCL  (N1,N2,N3,N4,NT)  FIXED  BIN; 

PUT  SKIP{3)  ('G  FOLICWS:'); 

POT  SKIP  (2)  ; 

?1T=N1  + N2+N3+N4; 

LCOONT=3 ; 

DO  TA=0  TO  N1; 

DO  IB=0  TC  N2; 

DO  IC=0  TC  N3; 

DO  ID=0  TO  N4; 

IF  LCCUNT=0  THEN  DO; 

PUT  SKIP; 

ICODNT=3;  END; 

PUT  EDITC  G{',IA,',»,IE,',  '.IC,  •,',ID, 

•)  = ',GP(IA,IB,IC,IC)) 

{A,F(N1/10+1),A,F(N:/10+1),A,F(N3/10>1),A, 
F(N4/10  + 1)  ,A,E  (15,7))  ; 

LCCUNT=LCOnNT-1 ; 

END; 

END; 

END ; 

END; 

END  POTG; 


END  TPOT; 
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APPENDIX  E:  EXACT  HNOEM  CCHEDTATION  FCB  COBPOTER  NETWORKS 

WITH  nOLTTPlE  JOB  TYPES 


HNORHU:  PROC  OPTIONS  (WAIN); 

/♦  COMPOTES  THE  NORM ALIZATICS  CONSTANTS  FOR  M NODES  AND 
OP  TO  POOR  JOE  CLASSES: 

A - 1ST  JOB  CLASS 
NA  - NOMBER  OF  A JOBS 
B - 2ND  JOE  CLASS 
NE  - NOMBER  OF  E JOBS 
C - 3Pr  JOB  CLASS 
NC  - NOMEER  OF  C JOBS 
D - 4TH  JOE  CLASS 
ND  - NOMEER  OF  C JOBS 

H (NA,NB,NC,ND)  IS  AN  ARRAY  FOR  THE  LAST  OB  HIGHEST  VALOE  OF 
M AND  IS  MAINTAINED  IN  STORAGE.  THE  M IS  IMPLIED. 

XA  (H)  ,XB  (M)  ,XC  (M)  ,XD  (M)  ARE  THE  RESPECTIVE,  RELATIVE 
OTILIZATICNS  FCB  JOB  CLASSES  A,E,C,C.  */ 

/♦  INPOT  IS  FREEPCRHAT: 

NO  OP  CLASSES,  NCIASS; 

M,  THE  NOMEER  CF  NODES; 

NC  CF  INCIVIDOAL  CLASSES,  NA,NE,NC,ND; 

RELATIVE  OTILIZATICNS,  XA,XB,XC,XD. 

THESE  ARE  VECTORS  CF  LENGTH  M; 

SERVICE  RATES  AT  THE  NODES  (1,2, 3, 4); 

A FLAG  FCB  PRINTER  OOTPOT  OF  TEE  H,  -Y  CR  N; 


A OOTPOT  FILE,  HDATA,  HOST  BE  ALLOCATED  FCR  THE 
G MATRIX.  V 

DCL  (LA,LB,  LC,LD)  FIXED  BIN; 

DCL  (M, NCLASS,NA,NB,HC,ND)  FIXED  BIN  ; 

DCL  (IA,IB,IC,ID)  FIXED  BIN; 

GET  FILE  (HDATA)  LIST (NCLASS , H)  ; 

POT  SKIP  DATA  (NCLASS, H)  ; 

DETERMINE  THE  NOMEER  OF  CLASSES  */ 


IF 

NCLASS=1 

THEN 

GC 

TO 

FINISH; 

IF 

NCLASS=2 

THEN 

GO 

TO 

FINISH; 

IF 

NCLASS=3 

THEN 

GC 

TC 

FINISH; 

IF 

NCLASS=4 

TEEN 

GC 

TO 

CLASS4 ; 
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/*  ESTABLISHED  EOP  4 CLASSES. 
OTHERS  ARE  SUBCASES  OF  4 AND  APE 
NOT  INCLUDEED  IN  THIS  LISTING.  */ 


CLASS4;  /♦  CASE  CF  FOUR  JCE  CLASSES  */ 

GET  FILE(HDATA)  LI  ST  {N  A , NB  , NC  , ND)  ; 

POT  SKIP  DATA  (NA,NB,NC, ND)  : 

BEGIN; 

DCL  (XA  (H)  ,XB  (f1)  ,XC  (P!)  ,XP  (f1)  ) FLOAT  EIN; 

DCL  (H(0:NA,0;NB,0:NC,0;ND)  , HC  {0:  NC  ,0;  H)  , HE  (0:  ND,0  ;H)  , 
HA  (0:NA,0:M)  , HB  {0:NB, 0: B) ) FLOAT  EIN; 

DCL  (HAB  (0:NA,0:NB,0:n),HAC(0:NA,0:NC,0:B), 

HBC  (0;NB,0:NC,0:n)  ) FLOAT  BIN; 

DCL  (HAT  (0:NA,0:ND,0:H)  ,HBB(0:NB,0:NC,0:H)  , 
HCD('):NC,0:NB,0:B))  FLOAT  BIN; 

DCL  (HABC(0;NA,0:NE,C:NC,0:H)  , 

HABD(0:NA,0:NE,0;ND,0:B)  , 

HACD(0;NA,0:NC,0;ND,0:B) , 

H’3CD(0;NB,0;NC,C:ND,0:H)  ) FLOAT  EIN; 

DCL  (U  {B)  ,UA  (B)  ,OB  (fl)  ,UC  (B)  ,OE  (B)  ) FLOAT  BIN; 

DCL  (T  (B)  ,TA(H)  ,TB  (B)  ,TC(H)  ,TD  (B)  ) FLOAT  BIN; 

DCL  BU(B)  FLOAT  BIN; 

DCL  (EA,EB,7C,ZD)  FLOAT  BIN; 

DCL  (TRA,TRB,TRC,TEB,TR)  FLOAT  BIN; 

DCL  (PA, PB, PC, PC)  FLCAT  BIN; 

DCL  {BNA,BNB, BNC,HND)  FLOAT  EIN; 

DCL  (VA,V3,VC  ,VD)  FLCAT  BIN; 

DCL  (BTPA,BTPE,BTPC,BTPD,HTPTCT)  FLCAT  BIN; 

DCL  PFLAG  CHAR (1) ; 

HEDIT:  FOFBAT  (A,F  (NA/10^1)  ,A,F(NE/10+1) ,A,F  (NC/10+1) , 

A,F  (ND/IO+T)  , A,E  (15,8)  ,SKIP)  ; 

BED  IT:  FOP  BAT  (A,  SKIP,  (4)  (X  ( 3)  , A , F (6  , 3)  ) ) ; 

GET  FILE{HDATA)  L I ST  ( X A , XB  , XC  , XD)  ; 

POT  SKIP  DATA {XA,XE,XC, XD) ; 

GET  FTLF{HDATA)  LIST  (BO)  ; 

POT  SKIP  DATA  (BO)  ; 

GET  FILE  (HDATA)  LI  ST  (Z  A ,Z  B , ZC  , ZD)  ; 

POT  SKIP  DATA  (ZA,ZE,ZC,ZD)  ; 

IF  NA=0  THEN  POT  SKIP  LIST(*CIASS  ERROR*) ; 

/•  INITIALIZE  THE  H */ 

H=1.');  HA=1.0;  HB=1.0;  HC=1.0;  HD=1.C; 

HAB=1.0;  HAC=1 .C;HEC=1. 0; 

HAD=1.0;  HBD^I.O;  HCD=1.0; 

HABC=1.0;  HABD=1. 0 ;HACD=1 .0;HECD=1.C; 


Appendix  E 


171 


i 


CALL  CNEJOB  (NA,](A,BA,ZA)  ; 

CALL  ONEJOE {NE,XB,aB,ZB)  ; 

CALL  ONEJOB (HC,XC,HC,2C) ; 

CALL  ONEJOB (ND,XD, EC,ZD)  ; 

CALL  TiOJOB  (BA,NB,XA,XE,HAB,HB,HA,ZA,ZB); 

CALL  THOJOB  (?IA,NC,XA,XC,HAC,HC,HA,ZA,ZC)  ; 

CALL  TWOJOB  (NB,NC,XB,XC,HEC,iiC,HE,ZB,ZC)  ; 

CALL  TWO JOB  (N A , MD, X A ,XD , HAD , HD ,HA , Z A ,ZD)  ; 

CALL  TWO JOE  (H E, N D , XB ,XD , EBD , HD , HB, Z E ,Z D)  ; 

CALL  TM0J03  (NC,»ID,XC,XD,HCD,HD,HC,ZC,ZD)  ; 

CALL  THREJOB  {NA,f)B,NC,XA,  XB,XC,HABC,HBC,  HAC,HAB,ZA,ZB,ZC)  ; 
CALL  THREJ0E(NA,NB,SD,XA,XB,XD,HABC,HBD,HAD,HAB,ZA,2B,ZD)  ,* 
CALL  THREJOE  (HA,HC,ND,XA,XC,XD,HACE,HCD,HAD,HAC,ZA,ZC,ZD) ; 
CALL  THREJOB  (NB, NC, ND, XB, XC , X D , HBC C , HCD , HBD , BBC , Z E,Z C, ZD) ; 

/*  CCHPDTE  THE  H */ 

CALL  FOOPJOB  (RA,NB,RC,ND,XA,XE,XC,XC, 

H,HECD,HACD, HABC,HABC,ZA,ZB,ZC,2C)  ; 

POT  SKIP  (2)  LIST  ('PRINTCDT  OF  H DESIRED?  - Y , N,  F OR  B •); 
GET  LIST  (PPLAG)  ; 

IF  {PFLAG='T'  PFLAG='B’)  THEH  CALL  POTH  (NA , NB, NC,RD, H)  ; 

IF  (PFLAG=»B'  PF1AG=»F')  THEN  DO; 

POT  FILE{H400T)  LI  ST  (N  A , NB  , SC  , N E)  ; 

POT  SKIP{2)  LIST('SHAT  IS  THE  CASE?*); 

GET  LIST  (NCASE)  ; 

POT  FILE(H40DT)  LIST  (NCASE); 

POT  FILS(H400T)  LIST  (Z A, ZE,ZC,ZC)  ; 

POT  FILE(H400T)  LIST  (HO)  ; 

POT  FILE(H40DT)  EDIT  (H) 

((4)  (E(15.8)  ,X(3))  ,SKIP)  ; 

END;  /*  THEN  DO  */ 

PDTSKIP(2)  EDIT  (•  H (' ,NA-1,  * , • ,NB,  • , * ,NC,  • , • ,ND,  ' )=  *, 

H (NA-1,NB,NC,ND))  (R  (HEEIT))  ; 

POT  SKIP(2)  EDITC  H (•  ,NA,  * , * ,NE-1  , * , * ,NC,»  ,*  ,ND,  •)=  *, 

H (NA,HB- 1,NC,ND))  (R(HEDIT)); 

POT  SKIP{2)  EDIT  («H  (•,NA,  * , •,NE,  *,  * ,NC-1  ,*  ,*  ,ND,  • )= 
H(NA,RB,HC-1,ND))  (R(HEDn)); 

POT  SKIP(2)  EDIT  (•  H (*  ,NA,  * , •,NB,  * , • ,NC,*  ,*  ,NC-1,  * )=  *, 

H (NA,NB,NC,ND-1) ) (R(HEEIT)); 

POT  SKIP(2)  EDIT  {'H  (*,NA,  * , '^NE,  « , * ,NC,*  , * ,ND,*)  = *, 

H (NA,NB,NC,ND) ) (P{HECTT)); 


/*  CCBPUTE  HEAR  THPOPDT  */ 


7A=NA;  VE=NB;  7C=NC;  VD=ND; 

nTPA  = H (NA-1,MB,NC,ND)/H  (N A , NB , NC, N D) ♦VA/ZA; 
HTPB=H (SA,NB-1,NC,SD)/H (NA,NB,NC,ND) *76/28; 
HTPC=H  (NA,NE,NC-1,BD)/H  (N A , NB , NC, N D)  ♦7C/2C; 
I1TPD=H  (NA,NB,NC,ND-1)/H  (NA,NB,NC,NC)  *7D/ZD; 
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PUT  SKIP(3)  EDITCSHCPT  CIPCOIT  THRUEO'"  FOLLCWS:*, 
’ MTPA=  *,!1TPA,»  1TPE=  ’,rtTPB,*  MTPC=  •,HTPC, 
• nTPD=  ' , STPC)  (F  (BEDTT)  1 


/♦  JOB  PPOBABILITT  FOLLOWS  */ 

HTPTOT=HTPA*f1TPE*MTPC  + *1TPD; 

PA=HTPA/MTPTOT ; 

PB=mPD/MTPTOT  ; 

PC=HTPC/MTFTOT; 

PD=«TPD/MTFTOT; 

POT  SKIP{3)  EDITCJCB  PROBABILITIES  APE:*,'PA=  • , PA, 
’PB=  •,PE,'PC=  •,PC,»PD=  * ,PD)  (F  IHEDIT)  ) ; 


/•  MEAN  QUEUE  LENGTHS  FOLLOW  */ 


HNOPH=H  (NA,NB,NC, NE)  ; 

aNA=VA* (1.0- (H (NA-1 ,NB,NC,NE)  /HNOHM)  ) ; 
f1NB  = VB*  (1. 0-  (H  (NA,  NE-1  ,NC,ND)  /HNORH)  ) ; 
f1NC  = VC*  (1. 0-  (H  (NA,NB,NC-1  ,ND)  /HNORH)  ) ; 

HND=VD* (1. 0- (H (NA, NE,NC,ND- 1) /HNORH) ) ; 

POT  SKIP(3)  EDIT(»eEAN  QUEUE  LENGTHS  ARE:’, 

'HNA=  •,MNA,’HNB=  •,HNE,'KNC=  •,HNC,'HNC=  •,HND) 
(P  (HEDIT) ) ; 

/*  RESPONSE  TIMES  */ 


TP A=HNA/HTPA; 

TEB=HNB/HTPE; 

TRC=HNC/HTPC; 

TPD=nND/HTPD; 

TP  = TPA*PA  ♦ TRB*PE  ♦ TRC*FC  ♦ TPE^PE; 

PUT  SKIP(3)  EDIT(*PEAN  RESPONSE  TIHES  ARE:', 

•THA^  •,TPA,''rRE=  *,THB,’TPC=  ',TRC,’TEC=  •,TBD) 
(P (HEDIT) ) ; 

PUT  SKIP  EDIT('TR=  ’ ,TR  ) (I  (3)  , A,  F (6 , 3) ) ; 


END;  /♦  BEGIN  BLOCK  CLASS3  */ 
GO  TO  FINISH; 


ONEJOB:  PR0C(N1,X1,H1,2;1)  ; 
DCL  N1  FIXED  BIN; 

DCL  XI (*)  FLOAT  BIN; 

DCL  21  FLOAT  BIN; 

DCL  HI  (*,*)  FLOAT  BIN; 

HI  (*,0) =1.0; 

HI (C,*) =1.0; 

HI (C,0) =1.0; 


i 

i 


Appendix  E 


173 


DO  lfl=1  TO  B; 

DO  IA=1  TO  N1; 

HI <LA,Ln)=Hl (LA,LB-1)  ♦ XI  (IB)  ♦LA /Z 1*H  1 (L A- 1 , LB)  ; 

END; 

END; 

END  ONEJOB; 


TWOJOB;  PPOC(N1,N2,X1,X2,H2,H2Z1,H2Z2,21,Z2) ; 

DCL  (N1,N2)  FIXED  EIN; 

DCL  (XI  (*)  ,X2  (*) ) FLOAT  BIN; 

DCL  (Z1,Z2)  FLOAT  BIN; 

DCL  (H2  (*,*,♦) ,H2Z  1 (*,♦) ,H2Z2  (♦,*) ) FLOAT  BIN; 

DO  LH*1  TO  B; 

H2 (*,0,LH) =H2Z2  (*,IB)  ; 

H2 (0,*,LH)=H271  (♦,LB)  ; 

DO  1B=1  TO  N2; 

DO  LA=1  TO  N1; 

H2(LA,LB,LH)  = H2  (L  A , LB , LH-1)  ♦ 

X1  (LB)  *B2  (LA-1,LB,1H)  *LA/Z  1 
♦ X2  (LB)  *H2  (LA,LB-1,IB)  *LE/Z2; 

END; 

END; 

END; 

END  THOJCB; 

THREJOB:  PROC  (N1,N2,N3,X1 ,X2,X3,H3, B321,H3Z2,H3Z3 ,Z1 ,Z2,Z3)  ; 

DCL  (N1,N2,N3)  FIXED  BIN; 

DCL  (XI  (♦)  ,X2  (♦)  ,X3  (*) ) FLOAT  EIN; 

DCL  (H3(*,^,^,^),H321(*,*.*), 

fl3Z2(*,*,*)  ,H3Z3  (*,♦,♦))  FLOAT  EIN; 

DCL  (Z1,Z2,Z3)  FLOAT  BIN; 

DO  LH=1  TO  fl; 

H3  (0,*,*,LB)=H3Z1  (*,*,LB)  ; 

H3(^,0,*,LB) =H3Z2(*,*,IR)  ; 

H3(*,*,0,LB)=H3Z3(*,^,LH)  ; 

DO  LC=1  TO  N3; 

DO  LB=1  TO  N2; 

DO  1A=1  TO  N1; 

H3  (1A,LB,IC,LH)=H3  (L A ,LB,LC,LH-1) 

♦ XI  (LB) ♦H3  (LA-1,LE, IC,LH) *LA/Z1 

♦ X2  (LB)  ♦H3  (LA,LB-  1,LC,LH)  *LE/Z2 

♦ X3  (LB)  *H3  (LA,LB,  IC-1,LH)  ♦LC/23; 

END; 

END; 

END; 

END; 

END  THREJOB; 
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FOnRJOE;  PROC  (N1,II2,II3,N4,X1,X2,X3,X4, 

H4,H4Z  1,H4Z2,H473,  H4Z4,Z1,Z2,Z3,Z4); 

DCL  (N1,N2,N3,tl4)  FIXED  BIN; 

DCI  (XI  {«) ,X2 (*) ,X3 <*) ,X4  (*) ) FLOAT  BIN; 

DCL  (H4  (♦,*,*,*)  ,h4Z1  (*,*,*,♦)  ,H4Z2 

H4Z3  (*,*,*,*)  ,H4Z4  (♦,*,*,*) ) FLOAT  BIN; 

DCL  (21,Z2,Z3,Z4)  FLOAT  BIN; 

DO  LB=1  TO  *!; 

H4  (0,*,*,*)=H4Z1  (*,*,*,!«)  ; 

(♦,0,*,*)=H4Z2(*,*,*,LH) ; 
H4{*,*,0,*)=H4Z3(*,*,*,LK) ; 

Htt (*,*,*,0)=H4Z4(*,*,*,LM) ; 

DC  ID=1  TC  N4; 

DC  LC=1  TC  N3; 

DO  LE=1  TO  N2; 

DO  LA=1  TC  N1; 

H4 (LA,LB,IC,LD) = 

H4 (LA,LE,LC, ID) 

♦XI  (LB)  ♦H4  (LA-1,LB,LC,  ID)  *LA/Z  1 

♦ X2  {L«)  *H4  (LA,1E-1,IC,  ID)  *LB/Z2 

♦ X3  (LM)  ♦H4  (LA,LE,LC-1,  LD)  *I,C/Z3 

, ♦X4  (Id)  *H4  (LA,LB,LC,LD-1)  *LD/Z4; 

END; 

END; 

END; 

END; 


END; 

END  FODEJOB; 

POTH:  PROO (N1, N2,N3,N4,HP) ; 

DCL  HP{*, *,♦,♦)  FLOAT  BIN; 

DCL  (N1 ,N2,N3,N4,NT)  FIXED  EIN; 

POT  SKIP  (3)  ('H  FOLICNS:'); 

POT  SFIP  (2)  ; 

NT=N  UN2  + N3  + N4 ; 

LCOONT=3; 

DO  TA=C  TO  N1; 

DO  IB=0  TC  N2; 

DC  IC=0  TC  N3; 

CO  ID=0  TO  N4; 

IF  LCCONT=0  THEN  DO; 
POT  SNIP; 
ICOONT=3;  END; 
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POT  IDIT(» 

•)  = '/HP 

(A,F(m/10+1),R,F(N2/104-1), 
F(Na/10>1)  ,A,E{13,6)  ,A) 
LC0DNT=LCCDNT-1; 

END; 

END ; 

END; 

END; 

END  POTH; 

FINISH;  END  HNOEBU; 


» ’,IC, ,ID, 

’) 

A,F(N3/10*1)  ,A, 


I 

I 


J 
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APPENDIX  F:  HNCPfl  HCDFI  FESPONSE  TIME  CORPOTATION 


HP.TinE:  PPCC  CPTICNS  (FAIN)  ; 

DCL  H (0:NA,0: NE,0: NC,0: ND)  FICAT  BIK  CTL; 

DCL  NCASE  FIXED  BIN; 

DCL  {2A,2B,ZC,ZD)  FLOAT  BIN; 

DCL  (NA,NB,NC,ND)  FIXED  BIN; 

DCL  HNOFfl  FLOAT  BIN; 

DCL  (TPA,TEB,TFC,TFD,TR)  FICAT  BIN; 

DCL  {HNB,nHC,PND,aNA)  FLOAT  BIN; 

DCL  (i*TPA,HTPE,!«TPC,  J1TPD,BTPTCT)  FLOAT  BIN; 

DCL  NO  (4)  FLOAT  BIN; 

HEDIT;  FORMAT  (A , P (NA/ 1C ♦ 1 ) ,A,F(NE/10  + 1),A,F (NC/IO+I) , 

A,F (ND/IO  + I)  , A,E (15,8)  ,SFIP) ; 

MEDIT;  FORMAT  (A, SKIP,  (4)  (X  (3)  ,A,F  (6,  5)  ) ) ; 

DCL  PFLAG  CHAP  (1)  ; 

/*  MOST  ALLOCATE  AN  INPUT  FILE,  RVALUES, 

CONTAINING  THE  NETWORK  FAPACETERS  ANt  H.  */ 

GET  FILE  (RVALUES)  L 1ST (N A , N B , NC, ND , I C AS E , Z A , ZB ,ZC ,ZD , MU)  ; 

ALLOCATE  H; 

GET  FILE  (RVALUES)  LIST(H); 

POT  SKIP(3)  EDITCVALOES  OF  H MATRIX  ARE:’,  ’ NA=  •,NA, 
•NB=  ',NE,'NC=  •,NC,’NC=  ',ND) 

(A, SKIP,  (4)  (A,F  (2)  ,X  (3))  ) ; 

PUT  SKIP(3)  EDIT  (•  NETWORK  C H A R ACTEP  I STICS  APE;', ’CASE  ', 
NCASE,’ZA=  •,ZA,'ZB=  ',28, ’20=  •,ZC,’ZD=  ’,20, 
'MU(1)=  • ,M0  (1)  ,’MU  (2)  = ',Mn(2), 

'MU(3)=  • ,MU  (3)  , 'MU  (4)  = ’,00(4)) 

(A,SKIP,A,F(2)  ,X  (3)  , (4 ) ( A , F (6 , 2 ) , X (3))  , SKIP, 

(4)  (A,E(15,7)  ,X  (2)  ))  ; 

POT  SKIP{3)  LISTCINPOT  THE  LOAD  VECTOR;'); 

GST  LIST (NA,NB,NC, ND) ; 

POT  SKIP  (2)  LIST  ('PRINTCOT  OF  H DESIRED?  - Y OR  N;')  ; 

GET  LIST  (PFLAG)  ; 

IF  (PPLAG='Y')  THEN  CALL  POTH  (N A , N E , NC, ND , H)  ; 

RESTART; 

POT  SKTP(2)  EDIT  (*H  (’ ,N  A-1,  * , ’ ,NB,  ’ , * ,NC,  • , • ,ND,  • )=  ’, 
H(NA-1,NB,NC,ND)) (P(HECIT)) ; 

POT  SKIP(2)  EDIT('H  (•,NA, ', •,NB-1, • ,',NC,',’  ,ND, •)=  ', 

H {NA,NE-1,NC, ND) ) (R(HEDIT)); 

POT  SKrP(2)  BDIT(’H{',NA, ',  •,NB,',  •,NC-1,',' ,ND,  •)=  ', 

H (NA,NB, NC-1 , ND) ) (R(HEEIT)); 

PUT  SKIF(2)  EDIT  (’H  (’ ,NA, ',  •,NE,', • ,NC, ’,ND-1, •)=  ', 
H(NA,N3,NC,ND-1))  (R(HEEIT)); 
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PUTSKIP(2)  EDIT(*H(',NA,*,*,NE,*,’,NC,','»NE,')= 
H (NA, NE, NC,ND) ) (P(HEEIT))  ; 

/*  COnPOTE  flEAN  THBOFOT  */ 


flTPA  = H (NA-1  ,NB,NC,  KE)/H  (NA,NE,NC,NE)  *NA/ZA; 
flTPB=H  (NA,NB-  1 ,NC,  ND)/H  (NA,NB,NC,NE)  •KE/ZE; 
J1TPC=H(NA,NB,NC-1,ND)/H  ( N A , NB  , NC , N D)  *NC/ZC; 

«TPD=H (N A,NB,NC,ND-1)/H  (N A , NB , NC, N E) *ND/ZD; 

POT  SKIP  (3)  EDITCSHCPT  CIPCOIT  THPUFDT  FCILCSS:  « , 

• HTPA=  •,flTPA,'  HTPB=  *,HTPB,'  RTPC=  ‘^HTPC, 

* HIPD^  • ,HTPE)  (P  (HEDIT)  ) ; 


/♦  JOB  PPOBABILITT  FOLLOWS  V 

BTPTOT=nTPA*nTPB+HTPC+f!TPC; 

PA=nTPA/HTFTOT; 

PB=HTPB/nTFTOT; 

PC=BTPC/HTPTOT ; 

PD=HTPD/HTFTOT; 

PUT  SKIP(3)  EEITCJCE  PROBA  EIIITIES  ARE:',»PA=  •,PA, 
•PB=  •,PB,'PC=  •,PC,'PE=  ',PD1  (f  (AEDIT)  ) ; 


/*  BEAN  QOEOE  LENGTHS  FOLLOW  ♦/ 


HNORB=H (NA,NB,NC, NE) ; 

HNA=NA*  (1.0-  (H  {NA-1,NE,NC,NB)/HNOHB) ) ; 

HNB  = NB*  (1.0-  (H  (NA,NB-1,NC,NE)/HNORB)  ) ; 

HHC=NC*  (1.0-  (H  (NA,NE,NC-1,HD)/HNOEH)  ) ; 

BND=ND*  (1. 0-  (H  (NA,NE,NC,ND-1)/'HHORB)  ) ; 

POT  SKIF{3)  EEIT(’HEAN  OOEOE  LENGTHS  ARE:*, 

*HNA=  •,HNA,*BNE=  »,HNB,'HNC=  •,BNC,*BNE=  '^BND) 
(P (BEEIT) ) ; 


/*  RESPONSE  TIBES  */ 

TRA=BNA/HTPA; 

TRB=BNB/BTPB; 

TRC=BNC/BTPC; 

TPD=HND/HTPD; 

TR  = TRA*PA  ♦ TRB*FB  ♦ TRC^PC  ♦ TRD*PE; 

POT  SKIP(3)  EDIT(*CEAN  PESPCNSE  TIBES  APE:', 

•TRA=  •,TRA,'TRB=  •,TRE,*TBC=  •,TPC,'TBE=  •,TRD) 
(P (BEDIT) ) ; 

POT  SKIP  EDIT(*TR=  * ,TR  ) (X  ( 3)  , A , F ( 6 , 3)  ) ; 


POT  PAGE; 

POT  SKIP(6)  LIST(’NEW  LOAD  VECTOP  IS:'); 
GET  LIST  (NA,NB,NC,NE)  ; 

IP  NA<100  THEN  GO  TO  RESTART; 
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POTH:  PRCC(N1,N2,N3,Na,HP)  ; 

DCL  HP(*, *,♦,*)  FLCAT  BIN; 

DCL  (N1,N2,N3,Na,NT)  FIXEC  EIN; 

POT  SNIP(3)(*H  FOLICMS:’); 

PUT  SKIP  (2)  ; 

NT=N1-»-N2+N3  + NU; 

LCOUNT=3; 

DO  IA=0  TO  N1; 

DO  IB=0  TO  N2; 

DO  IC=0  TC  N3; 

CO  ID=0  TO  N4; 

IF  ICCDNT=0  THEN  DO; 

POT  SKIP; 
lCOONT=3;  END; 

POT  EDTT(*  H {»,IA,  ' , IB, ',10,  • ,ID, 

•)  = • ,HP  {lA, IB, 1C, ID)  ,•;  •) 

(A,F  (N1/10*1)  ,A,F{N5/10+1)  ,A,F(N3/10  + 1)  ,A, 
F(NU/10  + 1)  ,A,E  (13,6)  ,A)  ; 
LC0DNT=LCCDNT-1; 

END; 

END; 

END; 

END; 

END  POTH; 

END  HFTIBE; 
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Addendum 

to 

EQN  Models  for  the  Analysis  & Design 
of  a 

Computer  Network  of  Functionalized  Processors 


Appendix  A; 


1 . In  procedure  SOLVE  the  second  EDIT  statement  after  NOHORE 
must  have  a semicolon  inserted  (;)  to  read: 

F {15,  6)); 

2.  In  subroutine  CHOOSE  the  following  lines  of  code  need 
to  be  inserted  after  "/*  THIS  IS  AN  ELIGIBLE  ROW*/": 

DO  J=1  TO  N-1 ; 

IF  TOP(J)  > 'M'  THEN 

/*  THIS  IS  THE  BEST  PIVOT  SO  FAR*/ 


p.  75:  Eq.  IV-43  should  read  Tr(l)  = (|M|-l)m. 


